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Scattering of X-Rays from Powdered Crystals 


G. E. M. JAUNCEY* AND ForD PENNELL, Department of Physics, Washington University 
(Received February 9, 1933) 


When monochromatic x-rays are scattered by a powdered 
crystal, the scattered rays consist partly of those rays 
which are scattered in special directions and which give 
rise to the Debye-Scherrer circles in a photograph and 
partly of those rays which are diffusely scattered and which 
give rise to the background between the circles. If the 
scattered rays from a powdered crystal enter a wide 
window of an ionization chamber, the intensity of the rays 
per unit solid angle is an average for several Debye- 
Scherrer circles and the diffuse background between the 
circles. Coven has experimentally examined this ‘‘average 
scattering’ for the case of monochromatic x-rays, but his 
results are inconclusive. In our experiments we have used 
not only a wide window (the extreme angular width of 
the rays entering the window was 6°) but also a wide band 
of wave-lengths obtained by passing the x-rays from a 
tungsten target x-ray tube operated at 54.9 kv through 
3.25 mm of aluminum or its equivalent. The spectral 
distribution of the intensity of these rays was obtained by 
reflection from rocksalt. The scattered intensities at aver- 
age angles varying from 5° to 90° from powdered crystal 
brickettes of KCl, CaS, NaF and MgO were compared 
with the scattered intensity at 90° from paraffin. From 


Jauncey and Harvey’s theory of the diffuse scattering 
from crystals together with the theory of the intensity 
of the Debye-Scherrer circles, a formula has been obtained 
for the ionization current produced by the rays entering 
a narrow window and consisting of a wide band of wave- 
lengths of known intensity distribution. This formula was 
then integrated graphically so as to obtain theoretical 
values for the case of a wide window. The calculation of 
the theoretical values involves f and f’’ values. James and 
Brindley’s f values as obtained from wave mechanics were 
used. Values of f’’ were calculated from data given in 
James and Brindley’s paper and are tabulated for F~, Nat, 
Cl-, K*, neon and argon. Excellent agreement between the 
theoretical and experimental values of the scattered in- 
tensity is shown over the whole range of average angles of 
scattering for each of the crystals. It is shown that, in 
experiments on the diffuse scattering of x-rays by crystals, 
an average wave-length as determined by absorption in 
aluminum may to a fairly close approximation be treated 
as if it were a single wave-length and so the S curves 
obtained for KCl by Harvey and for NaF by Jauncey and 
Williams are reasonably accurate. 





I. INTRODUCTION 


CCORDING to the theory as it has 
been developed by Jauncey, Harvey and 
Woo!: ?: 3. 4. 5 the intensity of the x-rays diffusely 


* The senior author was aided in part by a grant from the 
Rockefeller Foundation to Washington University for 
research in science. 

1G. E. M. Jauncey, Phys. Rev. 37, 1193 (1931). 

*G. E. M. Jauncey and G, G. Harvey, Phys. Rev. 37, 
1203 (1931). 

*G. E. M. Jauncey, Phys. Rev. 42, 453 (1932). 

*Y. H. Woo, Phys. Rev. 38, 6 (1931). 

*Y. H. Woo, Phys. Rev. 41, 21 (1932). 


scattered from a crystal consisting of atoms of 
one kind is given by 


S=Si+S:2/(1+a vers ¢)*, (1) 
Si= (f?— F*)/Z (2) 
= 1—f'""/Z*. (3) 


where 


and 


Furthermore, the intensity of the x-rays scat- 
tered from a monatomic gas is given by Eqs. (1) 
and (3) together with 


Si=f?/Z. (4) 
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For the definitions of the various quantities in 
Eqs. (1) to (4) the reader is referred to Jauncey’s 
paper.*® Jauncey and Harvey‘ have further shown 
that the intensity of the x-rays scattered by an 
inert gas is related to the intensities of the 
x-rays diffusely scattered and regularly reflected 
from the corresponding crystal by the formula 


Sgas= (S+ F?/Z)erystal- (5) 


The experimental S and F values for the right 
side of Eq. (5) have been obtained in different 
experiments and in different laboratories. It 
would be an advantage to devise an experiment 
in which at one and the same time the values of 
the sum of S and F*/Z could be obtained. In a 
private discussion with A. H. Compton it de- 
veloped that this might be accomplished by 
scattering the x-rays from a powdered crystal and 
by using a wide window in front of the ionization 
chamber so as to collect arcs of several Debye- 
Scherrer circles together with the diffuse scat- 
tering in between the circles. Following Compton’s 
suggestion Coven’ has scattered monochromatic 
x-rays of wave-length 0.71A from powdered 
crystals of KCl, NaF and MgO. Using Coven’s 
notation we shall represent the sum of the 
intensities of the diffuse scattering and of the 
Debye-Scherrer diffraction circles by S, and 
shall call this the total scattering. Coven found 
that the relation 


S: = ;. (6) 


holds very well for KCl and argon but does not 
hold nearly so well for NaF or MgO and neon, 
the S,,, values for argon and neon being those 
obtained by Wollan.$ 

Since the conclusions to be drawn from Coven’s 
experimental findings are uncertain and since the 
theory given in his paper is somewhat inade- 
quate, we have again attacked the problem of the 
total scattering from powdered crystals from 
both the experimental and theoretical sides. 
Furthermore, the essential point of Compton’s 
suggestion is that x-rays for a band of values of 
(sin 3¢)/X enter the window of the ionization 
chamber for a given setting of the spectrometer. 


®G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 38, 
1071 (1931). 

7A. W. Coven, Phys. Rev. 41, 422 (1932). 

8 E. O. Wollan, Phys. Rev. 37, 862 (1931). 
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This may be accomplished by using (1) a wide 
chamber window with a single wave-length, (2) a 
narrow window with a wide band of wave- 
lengths, or (3) a wide window together with a 
wide band of wave-lengths. Coven used the first 
method, while in the research to be described we 
have used the third method. We have measured 
the total scattering from powdered crystals of 


CaS, KCl, MgO and NaF. 


II. EXPERIMENTAL METHOD 


A specimen of the substance to be examined 
was finely powdered and pressed into a brickette 
in the hollow of a brass ring. The ring containing 
the brickette was then mounted on the axis of an 
x-ray spectrometer. In the case of both CaS and 
MgO, waterproof celophane was sealed over the 
two faces of the brass ring to prevent CO, and 
moisture from coming in contact with the sample 
and decomposing it. The general radiation from 
a tungsten tube was used to illuminate the 
brickette. The scattered rays were then measured 
by means of an ionization chamber and elec- 
trometer. 

An oil immersed x-ray outfit employing the 
balanced circuit as described by Bennett® was 
used to give a constant source of radiation. For 
reasons of convenience, the x-ray tube was 
mounted directly in front of the oil tank in a lead 
box. The maximum excitation potential as de- 
termined experimentally by the Duane-Hunt” 
relation was 54,900 volts and the current through 
the tube was 8.1 milliamperes. The ionization 
chamber was 35.2 cm long and was filled with 
methyl bromide at 74 cm pressure at a tempera- 
ture of 22°C. The diameter of the chamber was 
such that the x-rays travelled the whole length of 
the chamber without hitting the sides. The 
window was made of thin glass. 

The readings at an angle @ were made as 
follows: The paraffin was set at an angle of 45° to 
the main beam between the x-ray tube and the 
powdered crystal brickette which was mounted 
on the spectrometer axis in the Crowther 
position. The scattered intensity was measured 
at an angle ¢. The brickette and the paraffin were 


*R. D. Bennett, N. S. Singrich and W. C. Pierce, Rev. 
Sci. Inst. 2, 226 (1931). 

10 W. Duane and F. L. Hunt, Phys. Rev. 6, 166 (1915). 

J. A. Crowther, Proc. Roy. Soc. A86, 478 (1912). 
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then interchanged. The ionization chamber was 
moved to 90° and the scattered intensity from the 
paraffin measured. The reason for this procedure 
is that the intensity of the primary beam 
penetrating the substance in the scattering posi- 
tion is the same irrespective of whether the 
paraffin or the crystal brickette is in that 
position. 

In a preliminary experiment, the crystal 
brickette was removed and the paraffin slab set 
so as to scatter x-rays into the ionization chamber 
at an angle of 30°. Sheets of aluminum were 
placed between the paraffin and the x-ray tube 
and the logarithm of the intensity of the x-rays 
entering the chamber was plotted against the 
thickness of the aluminum. For thicknesses of 
aluminum greater than 3.25 mm the graph as 
shown in Fig. 1 becomes a straight line. In the 
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Fic. 1. Absorption in aluminum. 


actual scattering experiments, a thickness of 
aluminum was placed in the main beam such 
that the absorption in the aluminum plus that in 
the crystal brickette and that in the paraffin slab 
was equal to the absorption in 3.25 mm of alumi- 
num. The mass absorption coefficient in alumi- 
num for the straight portion of the graph in Fig. 
1 is 0.952 per g/cm*. Using Compton's tables,'” 


® A. H. Compton, X-Rays and Electrons, p. 184. 


this corresponds to a wave-length of 0.39A. 
However, we have not used this average wave- 
length in the way described by Harvey" and by 
Jauncey and Williams" but have instead obtained 
the distribution of the energy in the spectrum of 
the x-rays which have passed through 3.25 mm of 
aluminum. 

The energy distribution with respect to wave- 
length was obtained in a separate experiment. 
We reflected the x-rays from a (1, 0, 0) face of a 
single crystal of rocksalt into the same ionization 
chamber as was used in the scattering experiment. 
The integrated intensity'® of x-rays of wave- 
length \ as reflected by a mosaic crystal is given 
by 

Ww And'F? e* 1+ cos? 20 
— n? ep RINE ’ (7) 
I 4u m’c* sin 20 





where @is the glancing angle for the reflected 
rays. Eq. (7) holds for monochromatic rays when 
the crystal is rotated at an angular speed w. W 
is the total energy reflected into the ionization 
chamber while the crystal is being turned through 
a reflection position for the wave-length \. The 
purpose of the rotation of the crystal is to bring 
all the minute crystals of the mosaic crystal into 
the reflection position. However, when the con- 
tinuous spectrum is being analyzed the difference 
in the orientations of the minute crystals of the 
mosaic causes wave-lengths of a range \ to 
A+dd to be reflected and consequently each 
minute crystal is able to find some wave-length 
to reflect when the face of the large crystal is set 
at a glancing angle @ with the primary beam. 
Hence, at least approximately, 


W =const. X (J\* F?/)(1+ cos? 26)/sin 20, (7a) 


where W is the intensity of the x-rays in the 
wave-length range \ to A+dd reflected into the 
ionization chamber, u is the absorption coefficient 
of rays of wave-length \ in the crystal and F is 
the atomic structure factor for the crystal. Fis a 
function of (sin @)/A. We have used the F values 
for rocksalt as given by James and Firth.'* The 


3G, G. Harvey, Phys. Rev. 38, 593 (1931). 

4G. E. M. Jauncey and P. S. Williams, Phys. Rev. 41, 
127 (1932). 

6 See reference 12, p. 128. 

‘6 R. W. James and E. M. Firth, Proc. Roy. Soc. A117, 
62 (1928). 
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values of » used were those for the ordinary 
absorption coefficient in rocksalt. We made no 
attempt to correct for extinction. After obtaining 
values of W, Eq. (7a) was solved for J, the 
intensity of the x-rays in the range \ to A\+d\ 
existing in the primary beam. Values of J were 
then plotted against \. Next, it was necessary to 
correct for second and third order reflection of 
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Fic. 2. Distribution of the intensity in the continuous 
spectrum. For ease in calculation the curve is drawn to 
meet the wave-length axis at \=0.68A. 


x-rays. The shortest wave-length present was 
0.225A, so that second order reflection began at 
an apparent wave-length of 0.450A. For re- 
flection from a given set of planes in a given 
crystal—for instance, the (1, 0, 0) planes—F is a 
function of the order of reflection. Hence, 
knowing the value of J as calculated from the 
first order reflection of, say, \=0.30A, the values 
of W for the second order reflection of \=0.30A 
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can be obtained from Eq. (7a). This value is 
subtracted from the experimental value of W 
found for a value of @ corresponding to first order 
reflection of \=0.60A. The remaining value of W 
is the true value of W for \=0.60A. Also at 
A= 3X0.225=0.675A third order reflection be- 
gins and a further correction must be made. The 
final corrected curve for the spectral distribution 
is shown in Fig. 2. It should be noted that this 
curve is for the spectral distribution of the x-rays 
as absorbed in the ionization chamber. The mass 
of the methyl bromide in the chamber is not 
sufficient to absorb completely x-rays of all the 
wave-lengths in the spectrum. However, the 
same ionization chamber was used as in the 
scattering experiments. The scattered rays of 
various wave-lengths are absorbed to the same 
extent in the methyl bromide as the reflected 
rays of the same wave-lengths in the spectral 
distribution experiment. We shall see later that 
the spectral distribution as measured by ab- 
sorption in the chamber is all we need to 
know. 

As previously mentioned, Harvey" and Jauncey 
and Williams" in their experiments measured the 
absorption coefficient of the x-rays in aluminum 
and from this coefficient determined the average 
wave-length and then used this average wave- 
length as if it were the wave-length of mono- 
chromatic x-rays. From the spectral distribution 
curve of Fig. 2, the average wave-length was 
found to be 0.41A, while the absorption coeffi- 
cient obtained from the straight portion of the 
curve of Fig. 1 leads to an average wave-length 
of 0.39A. We submit that the agreement of these 
two average wave-lengths is reasonably good. 
We shall return later to the question of the use of 
an average wave-length as if it were a mono- 
chromatic wave-length. 


III. ForRMULAS FOR CALCULATION 


Referring to the paper by Jauncey and 
Williams," we see that the energy of ionization 
produced per second is not equal to the rate of 
flow of x-ray energy into the ionization chamber 
unless all of the x-rays are absorbed in the gas of 
the chamber. When the absorption is incomplete, 
the situation is expressed by Eq. (20) of the above 
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paper. Then, eliminating s,; and sz by means of Eqs. (21) and (22) of that paper, we obtain 


ACt 





C.- ; ; 
' R?cosi¢ W mic 


where C, and C are the ionization currents 
produced respectively by the scattered and the 
transmitted rays, S; and S: are respectively the 
coherent and incoherent scattering factors, Ko 
and K, are respectively the fractions of the 
coherent and incoherent rays absorbed in the 
chamber, 7; is a quantity introduced by Jauncey 
and DeFoe" into the Crowther" formula for 
scattering by a slab in order to take account of 
the extra absorption of the incoherent scattered 
rays in the slab itself, A is the area of the 
jonization chamber window, and R the distance 
of the window from the scattering slab. 

Eq. (8) is for a given wave-length. However, in 





A NZpt 1+cos?¢@ e! 
° R? W = 2cos 3 mec4 





D 


For the wave-lengths used the x-rays scattered by 
paraffin at 90° were almost entirely incoherent, 
as pointed out by Harvey"™ and Coven.’ Conse- 
quently for scattering from paraffin at 90° we 
obtain an equation similar to Eq. (9), but with 
i=0 and S:=1. We shall prime (’) those 
quantities in this equation which specifically 
refer to paraffin. Then, dividing Eq. (9) for the 
crystal by the similar equation for paraffin, we 

obtain 
Ds, ZptW’ (1+ cos? ¢) cos 45° B 


, ne (10) 
cos 5¢ B 


D' y¢° Z' p't'W 








where B is an integral for the crystal and B’ is the 
corresponding integral for the paraffin. 
We shall now consider B’, where 


K'90°T' 90° 
B= { —— 1a. 
K,(1+a)§ 


According to Jauncey and Williams," 


(11) 


7G. E. M. Jauncey and O. K. DeFoe, Phil. Mag. 1, 
711 (1926). 


NZp e* 1+cos’?¢ {s SoK4T¢ 





(8) 
Ko(1+a vers ¢)? 


our experiments we used the band of wave- 
lengths shown in Fig. 2. The intensity in the 
range dd at d is J(A)dd where J(A) is proportional 
to the ordinate of the graph and hence we replace 
C in the right side of Eq. (8) by Z(A)dd. Also we 
shall write the left side in the form Cydd because 
C;dd is the ionization current produced by the 
scattered rays which are due to the primary 
rays in the range d) at \. To obtain the ionization 
current D, when the whole spectrum of wave- 
lengths is scattered, we integrate over all wave- 
lengths present in the primary beam. Several of 
the quantities on the right side of Eq. (8) are not 
functions of \ and so we obtain 








SoK 4T¢ 
f {Si |1an (9) 
Ko(1+ a vers ¢)' 
T=(1—e-*)/g, (12) 


where 


g=(u2—)t/cos 3¢ (13) 


and ye and yw; are the linear absorption coeffi- 
cients of the incoherent and coherent scattered 
rays, respectively, in the scattering slab. The 
value of p’t’ for the paraffin slab was 0.439 g/cm. 
From Compton’s absorption tables” we have 
calculated the absorption coefficient of paraffin 
for various wave-lengths. Putting values in Eqs. 
(12) and (13) we find that 7’90° varies from 0.990 
to 1.000 over the range of wave-lengths 0.7A to 
0.2A and so may be regarded as practically a 
constant whose value is 0.995. If mw; is the 
absorption coefficient of wave-length \ and pe 
that of wave-length \+ 6A in methyl bromide, the 
fractions absorbed are 


Kyo=1-—e™! (14) 
and 
Koo? = 1—e-", (15) 


where ¢ is the length of the gas in the ionization 
chamber. The change of wave-length 6A is 
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0.0242A, the Compton change at 90°. Consider 
the quantity 


y = Koo"/Ko(1+a)*. 


The graph of y plotted against \ is shown in Fig. 3 
and it is seen that y has the practically constant 
value of 0.884 over the range A=0.25A to 
\=0.65A, which is the range of the wave-lengths 
of appreciable intensity in the spectrum, as 
shown in Fig. 2. Hence the factor of J in the 
integrand of the right side of Eq. (11) is constant 
with respect to the wave-length and so may be 


ffs 


after putting in numerical values. For con- 
venience we shall denote the value of the right 
side of Eq. (17) when experimental values are 
inserted by S.xp. It is interesting to note that this 
same numerical formula for S.x, is obtained if 
either of the average wave-lengths mentioned in 
the last paragraph of §2 had been used. 


(16) 





aK.T. 





{ran /f'1an=0.314% . — 
Ko(1+e vers ¢)? D'yo° 1+cos? @ Zot 
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Fic. 3. Graph of Kg/Ko(1+a)* against wave-length. 


taken outside of the integral sign and we have 
B’=0.880 f Idd. Solving Eq. (10) for B/ fJIdh, 
we obtain 

Ds 


cos3}¢ W 





(17) 





IV. EXPERIMENTAL RESULTS 


Values of the experimental ratio Ds/D'99° for 
the respective powdered crystals were obtained 
and substituted in the right side of Eq. (17) and 
the corresponding values of Six, calculated. 
These values of S.x, are shown in Tables I and 
II. The values of the mass per unit area of the 


TABLE I. Total scattering from powdered crystals of KCl and CaS. Values of Dg/D'90° and S.zp. 












































KCl CaS i KCl ~ CaS 
. @ D, D : , 
> sin 3 ve Sexp ve Ss xp ? sin s - Sexp “ee Sexp 
5.0° 0.044 2.64 3.97 18° 0.156 3.27 5.10 3.22 4.36 
6.0° .052 5.02 7.59 20° .174 2.87 3.92 
7.0° .061 6.51 9.85 FS 191 2.78 3.83 
7.5° .065 6.00 7.95 24° .208 2.29 3.67 2.41 3.36 
8.0° .070 6.72 10.19 6.36 8.39 30° .259 1.80 2.99 1.99 2.87 
8.5° .074 6.89 10.46 6.60 8.70 36° .309 1.44 2.50 1.59 2.39 
9.0° .078 6.67 10.15 6.80 8.96 42° .358 1.13 2.05 1.32 2.07 
9.5° .083 6.56 9.96 6.60 8.70 48° 407 0.948 1.80 
10.0° .087 6.05 9.20 6.49 8.57 54° 454 0.773 1.54 0.876 1.52 
10.5° .092 5.71 8.71 60° .500 0.656 1.37 
11.0° .096 5.79 8.77 5.95 7.88 66° .545 0.632 1.31 0.707 1.33 
12.0° .105 3.15 7.85 5.37 7.13 a .588 0.505 1.12 
14.0° .122 4.45 5.94 78° .629 0.510 1.14 0.599 1.16 
15.0° 131 3.94 6.09 84° .669 0.500 1.12 
16.0° .139 3.68 4.95 90° .707 0.470 1.00 0.597 1.10 
TABLE II. Total scattering from powdered crystals of NaF and MgO. Values of Dg/D'y0° and Sezp. 
NaF MgO o NaF - MgO 
‘ D D.  & Py 
1) sin S De Sexp ve. Sexp 7) sin 9 ve Sexp D’ soo Sexp 
9° 0.078 9.21 4.55 30° 0.259 2.63 1.42 1.06 1.63 
10° .087 10.17 5.04 2.71 3.81 36° .309 2.19 1.25 0.837 1.34 
ty .096 10.25 5.08 2.96 4.17 42° .358 1.98 1.17 .738 1.23 
* .104 9.84 4.90 3.14 4.44 54° 454 1.54 1.11 .608 1.12 
15° .130 6.97 3.50 2.58 3.67 66° 545 1.49 1.05 .529 1.06 
18° .156 5.34 2.71 1.90 2.75 78° .629 1.37 1.00 497 1.02 
24° .208 3.54 1.85 1.32 1.95 90° .107 1.40 0.97 486 0.96 
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various powdered crystal brickettes were as 
follows: NaF, 0.678; MgO, 0.229; KCI, 0.218; 
CaS, 0.242 g/cm*. The width of the window of 
the ionization chamber subtended an angle of 4° 
at the axis of the spectrometer. The widths of the 
defining slits in the primary beam were such that 
the total variation in ¢ for a given setting of the 
ionization chamber was about 6° so long as the 
average angle was less than 30°. In addition to 
the defining slits, screening slits were added so as 
to prevent rays scattered by the defining slits 
from entering the ionization chamber. These 
screening slits are very important at small 
angles. 


V. COMPARISON WITH THEORY 


The left side of Eq. (17) contains S; and S, and 
these quantities* are dependent upon f and f” as 
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well as upon F. The atom form factors f and f” 
can be calculated from wave mechanics. Tables 
of f values so calculated have been published by 
James and Brindley.'* Furthermore, James and 
Brindley give tables of fy values. The symbol fo 
refers to the same quantity as does the symbol E 
used by Jauncey':* and Woo,' and hence values 
of f” have been calculated from James and 
Brindley’s tables by means of the formula*® 


f"=ZDEP. (18) 


Because calculated values of f’’ may be of use to 
others working in this field, these values are 
given in Table III. The values for neon and for 
argon in Table III are the averages of the 
corresponding values for F~ and Na* and for Cl- 
and Kt, respectively. We have used Froman’s 


TABLE III. Atom form factors. Values of f and f”. 




















sin }¢ F- Nat Neon a Kt Argon 
r f } f - f i f : ¥ Da f ca 
0.00 10.00 10.00 10.00 10.00 10.00 10.00 18.00 18.00 18.00 18.00 18.00 18.00 
a 8.73 8.76 9.55 9.58 9.14 9.17 15.20 15.50 16.40 16.50 15.80 16.00 
Pr 6.69 6.89 8.18 8.23 7.44 7.56 11.50 13.10 13.36 14.05 12.43 13.57 
a 4.77 5.38 6.62 6.82 5.70 6.10 9.34 12.00 10.78 12.53 10.06 12.26 
4 3.52 4.63 $5.20 5.71 4.39 5.17 7.98 10.92 8.80 11.50 8.39 11.21 
5 2.74 4.29 4.08 4.88 3.41 4.59 7.16 9.97 7.74 10.67 7.45 10.32 
6 2.23 4.10 3.22 4.43 2.73 4.27 6.47 9.01 7.04 9.82 6.75 9.42 
R 1.82 3.92 2.62 4.15 2.22 4.04 5.86 8.13 6.44 8.95 6.15 8.54 
8 1.60 3.76 2.22 4.00 1.91 3.88 5.12 7.36 5.91 8.18 5.51 7.77 
9 1.53 3.59 1.94 3.84 1.73 3.72 4.41 6.75 5.32 7.56 4.86 7.16 
1.0 133 342 1.73 3.68 1.63 3.55 3.77 6.25 4.72 6.88 4.25 6.56 
1.5 0.96 2.14 2.23 3.86 
2.0 0.59 1.32 1.32 2.35 
3.0 0.28 0.62 0.60 1.14 























method" for extrapolating the f values beyond 
(sin 3@)/A= 1.0. For the extrapolation of the f” 
values we have applied Froman’s method to 
James and Brindley’s E (or fo) values for each 
type of electron and have then used Eq. (18). 
For a given value of ¢, S; and S: are functions 
of A alone. Since in the case of the diffuse 
scattering from a single crystal S,; and S» are 
given by Eqs. (2) and (3), respectively, it is seen 
that S,; and S:; may be calculated from wave 
mechanics. In our experiments we have not used 
single large crystals but have scattered x-rays 
from brickettes of powdered crystals and so it is 





*R. W. James and G. W. Brindley, Phil. Mag. 12, 
81 (1931). 





necessary to find formulas for S; and S2 for our 
experimental case. The formula for S2 remains as 
given in Eq. (3), but the formula for S; is 
changed. We now proceed to the derivation of the 
formula for S; for the case of a powdered 
crystal. 

Referring to Eq. (5.23) on p. 131 of A. H. 
Compton’s X-Rays and Electrons, we see that for 
a powdered crystal brickette placed in the 
Crowther position 

P, N*p*\'F* et 1+cos?¢ pltp’ 


P W? mic* 2sin®?@ 4rRp 
where P, is the intensity of the Debye-Sherrer 
19D. K. Froman, Phys. Rev. 36, 1339 (1930). 





» (19) 


| 
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circle diffracted into an ionization chamber 
window whose height is /, P is the intensity of the 
primary beam penetrating the brickette, p is the 
density of a little crystal in the brickette and p’ 
the overall density of the brickette. Now the 
intensity J, of the coherent part of diffusely 
scattered x-rays entering the chamber window is 
given by 
e* 1+ cos? ¢ p’tS; 
Giecenens , (20) 
P RF W mic 2 cos $@ 
where S; is given by Eq. (2). If both a Debye- 
Scherrer circle and the coherent part of the 
diffuse scattering on either side of the circle 
enter the ionization chamber window, Eqs. (19) 
and (20) combine to give 
Iyt+P, A NZ e 1+¢cos?¢ p’t 


W mic4 2 














P R cos @ 
Npxr* F? pLR 
| Sit fae , (21) 
W Z 167A sin? 6 cos 6 


where 0= 3¢. Hence, if we represent the coherent 
scattering factor for a powdered crystal by Sy’, 
we have 


Npd;? F? PAR 
W 2Z 167A sin? 0, cos 6, 











S$,’ =S\,+ (22) 


Because the second member of Eq. (22) depends 
upon the particular Debye-Scherrer circle which 
enters the chamber, various quantities are sub- 
scripted with r. For the r circle we have 


A, = 2d, sin 0, (23) 
where, for a simple cubic crystal, 
d,=d,/(r)'. (24) 
The counter 7 is given by 
r=P+P+7’, (25) 


where a, 6 and y are integers or zero but not 
more than two of them may be zero. The counter 
r according to Eq. (25) may be any integer 
excepting 7, 15, 23, 28, 31 and so on. But for a 
simple cubic crystal the principle grating space is 
given by 


» — di=(W/Np)}. (26) 


Now the area A of the ionization chamber 
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window may be replaced by /w where w is the 
width of the window. Furthermore, the width w 
may be replaced by Ré¢, where R is the distance 
of the window from the axis of the spectrometer 
and 4¢ is the angle subtended by the width of the 
window. Hence, after eliminating \, from Eq. 
(22) by means of Eqs. (23), (24) and (26), we 
obtain 





Z (n)} 2nbb 27) 


This quantity S;’ replaces S, in the integrand in 
the numerator of the left side of Eq. (17). 
However, since, for a given angle ¢, F, only has 
values for particular values of \, the integral as 
applied to the second term on the right of Eq. 
(27) becomes a summation. In this summation 
the Z(A)d\ which appears under the integral is 
replaced by J(A,)6A,. The quantity 6d, is due to 
the angular width 6¢ of the window. Differ- 
entiating Eq. (23) and then dividing the equation 
so obtained by Eq. (23), we obtain after re- 
arrangement 


5\,= 2A, cot 0,56. (28) 


The left side of Eq. (17) in virtue of Eqs. (2), (3) 
and (28) becomes 


nnd ~ 8+ L, (29) 


where 


2 (1—f"/Z)KeT , 
a= f F Macs “tran / | Idd (30) 
Z Ko(ite vers ¢)* . 


and 





I= [ (/zytax / [tan (31) 


jl _ F2 padlQr.) . 
io me, piel see Idv. (32) 
lar * Z  (r)3 


The numerical value of the right side of Eq. (29) 
is determined by the wave-mechanics values of f 
and f’’ and hence we shall designate this nu- 
merical value by Sm. If our theory is correct, 
Sexp as given by the right side of Eq. (17) should 
equal Sy. We have not used the theoretical 
values of F in the calculation of Sj, because of 
the uncertainty of the Debye-Waller temperature 
factor, but have used the values of F as experi- 


and 
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Fic. 4. Curves I, F?/Z averaged over the wave-length 
band shown in Fig. 2 and an ionization chamber window 
width of 6°; curves II, L (see Eq. (32)) averaged over the 
wave-length band; curves III, L averaged over the wave- 
length band and the window width. 


mentally obtained for KCl by James and Brind- 
ley,”? for NaF by Havighurst,”! and for MgO by 
Froman.” : 

Values of /7 as defined by Eq. (31) have been 
calculated for a given crystal for various values 
of ¢ and have been plotted against ¢. Each of 
these values of J/ is for a small range of angles 
about ¢. However, in our experiments, the range 
of angles was intentionally large, being about 6° 
for a given setting of the ionization chamber. 


**R. W. James and G. W. Brindley, Proc. Roy. Soc. 
Al21, 155 (1928). 

* R. J. Havighurst, Phys. Rev. 29, 1 (1927). 

* D. K. Froman, Phys. Rev. 36, 1330 (1930). 


We have therefore obtained the value of H as 
averaged over a range of 6° for each average 
value of the angle. The graphs for these average 
values of ZZ are designated by the Roman 
numeral I in Fig. 4. Likewise, L as defined by Eq. 
(32) was plotted against @ (curve II) and from 
this curve average values of L were obtained. 
The graphs for these average values of L are 
designated by the Roman numeral III in Fig. 4. 
It is seen that for a given crystal there is an 
angle above which (L—#2) is nearly zero. For 
angles above this particular value, therefore. 
Sm=G approximately. But G, in virtue of Eq. 
(30), is the value of S,, for the corresponding 
inert gas. Hence the relation as expressed by 
Eq. (6) should approximately hold for angles 
greater than the particular angle for each crystal. 
As the angles are made smaller than the par- 
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Fic. 5. Scattering of x-rays by powdered crystals. 
Curves I, theoretical S values for argon or neon; curves 
II, theoretical S values for powdered crystals (see Eq. 
(29)); black circles, experimental points. 
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ticular angle, the value of (L—#J) first increases, 
then decreases, and finally becomes a large 
negative number. At zero angle, (L—H) very 
nearly equals — Z. 

Using the values of f and f” for neon and argon 
shown in Table III, we have calculated the 
values of G as defined by Eq. (30). As for the 
graphs of the average values of HJ and L, the 
values of G were averaged over an angular width 
of 6°. These average values of G are plotted 
against sin }¢ and the resulting graphs are 
designated by the Roman numeral I in Fig. 5. 
Graphs of Sj, averaged over an angular width of 
6° are now obtained and are designated by the 
Roman numeral II in Fig. 5. Since G very 
nearly equals +Z at zero angle, Sj, in virtue of 
Eq. (29) becomes Z—Z or zero at zero angle. 

Values of S.x, are shown as black circles in Fig. 
5. It is seen how closely these circles fall upon the 
respective S,, curves. Particularly is this the 
case for KCl where the experimental points were 
taken down to angles as low as 5°. The experi- 
mental points follow the theoretical curve up to 
the maximum and then follow it by rapidly 
dropping as the angle approaches zero. An St, 
curve for CaS is not shown because F values for 
this crystal were not available. 


VI. CoNCLUDING REMARKS 


We feel that the agreement between our 
theory and the experimental results as shown by 
the curves and experimental points in Fig. 5 is 
excellent. Such discrepancy as appears is proba- 
bly due to numerical approximations in the 
theory and to experimental errors. 

With regard to Coven’s experiment’ with 
monochromatic rays, we believe that the relation 
expressed by Eq. (6) is only approximately true 
for particular angles of scattering, but is not 
generally true. The argument as developed in 
Eqs. (19) to (27) leads to an equation similar to 
Eq. (27), but with a summation sign in front of 
the second term on the right side of this equation. 
For monochromatic rays and a wide window, we 
arrive at 

F? FP pr 
S:=Sgaa—— +> J P . 
Z * Z (r)i 2rd 


tan 6, 











(33) 


It is seen that the discrepancy between S,; and 
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Sgas is equal to the algebraic sum of the second 
and third terms on the right side of Eq. (33), 
Values of the second term with a positive sign 
have been plotted against the angle for the case 
of KCI and \=0.71A as shown by the curve in 
Fig. 6. Values of the third or summation term in 
Eq. (33) have been calculated for particular angles 
and are shown in Fig. 6 as black dots and white 
circles. The black dots refer to the case for a 10° 
window and the white dots to that for a 7° win- 
dow. The curve in Fig. 6 depends also upon the 
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Fic. 6. Scattering of \=0.71A from a powdered crystal 
of KCI. Curve, values of F?/Z averaged over a window width 
of 6°; circles, values of third term on right side of Eq. (33) 
for different window widths. 


width of the window but only slightly. The curve 
shown is for a 6° window. If a 7° window is used, 
the discrepancy between S, and S,,; for a given 
angular setting of the ionization chamber is 
equal to the difference between the ordinate of a 
white circle and that of the curve for this setting. 
For a 10° window, the discrepancy is equal to the 
difference between the ordinate of a black circle 
and that of the curve. It is seen that, generally 
speaking S; is not equal to Sgas. 

In previous experiments performed in this 
laboratory on the diffuse scattering of x-rays 
from crystals an average wave-length has been 





—__—— 

















SCATTERING OF X-RAYS 515 


determined from an absorption curve such as 
shown in Fig. 1. This average wave-length has 
then been treated as a single wave-length and 
values of S plotted against (sin $¢)/A. There has 
been some objection to this procedure. However, 
we have pointed out at the end of the last 
paragraph of §2 that the average wave-length as 
obtained from the spectral distribution is very 
nearly that obtained from absorption measure- 
ments. Also, at the end of §3 it is pointed out that 
the numerical formula for S.x» is very nearly 
independent of the precise wave-length assumed 
if the x-rays are supposed to consist of a single 
wave-length. Hence, in experiments such as those 
of Harvey" and of Jauncey and Williams," the 
values of S are reasonably accurate. The trouble, 
however, arises in the values of (sin 3¢)/A to 
which the respective values of S should be as- 
signed. As a result of our present investigation it 
seems that, excepting in those portions of the 
graphs of S versus (sin }¢)/ where the slope is 
large, the values of (sin }¢)/A may be calculated 
from the angle and that value of \ which is found 


from the absorption curve. It seems, then, that a 
beam of x-rays which has passed through a 
sufficient thickness of aluminum acts for scatter- 
ing purposes very nearly as if it were made up of 
monochromatic rays. The homogeneity of x-rays 
which had passed through a sufficient thickness 
of aluminum was noted several years ago by C. 
G. Barkla. The S curves shown in the papers by 
Harvey and by Jauncey and Williams are 
therefore reasonably accurate, excepting at small 
angles. As a result of this present investigation, 
we now believe that the proper procedure is to 
obtain the distribution of energy in the spectrum 
of the x-rays used and to calculate from theory 
the scattering to be expected at each angle. 
These theoretical values are then compared with 
the experimental values. For convenience the 
curve of the theoretical values and the points 
representing the experimental values are plotted 
against (sin 3¢). 

In conclusion, we wish to thank Dr. J. A. Van 
den Akker for much help and advice in regard to 
the assembling of the all-in-oil x-ray outfit. 
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The (1-1-2) and (1-0-0) planes of a tungsten crystal 
were bombarded at normal incidence with primary elec- 
trons and the intensity of the full-velocity secondary 
beams measured as a function of azimuth, co-latitude (6), 
and primary voltage. A new magnetic deflection method 
of analyzing the secondaries permitted observations at co- 
latitudes down to zero. The crystal was outgassed 1550 
hours at temperatures up to 1600°C at pressures of 10-7 
mm of Hg or less. Strong sharp beams were observed in 
the AA’ azimuth (Fig. 1) of the (1-1-2) plane at all wave- 
lengths which were used, and were found to be governed in 
every case by the volume equation m.\=d/6)+-2(d/6!) 
sin (30°—@) and in no case by the surface equation. This 
is the first time that such a pronounced deviation from the 
usual theory has been observed. The final average experi- 
mental value for W, was 5.52 volts which, combined with 
photoelectric and thermionic data yields a value for W; of 
about 1 volt. The values obtained for W. were unusually 
consistent. A few beams of the usual type obeying both 
interference equations were observed in the B and C 


azimuths of the (1-1-2) plane and also in an azimuth 30° 
from B. They were all quite broad, however, and con- 
sequently no attempt was made to estimate W, from them. 
Beams were found in the case of the (1-0-0) plane which 
could be detected for various wave-lengths within the 
immediate neighborhood of their predicted location and 
which showed a tendency to vary in latitude in obedience 
with the volume equation, but became too weak to observe 
if the wave-length chosen was very far from the predicted 
value. Slight variations in the angle of incidence were 
found to have a very great effect on the observations. 
Some of those portions of the observations which differ 
from typical results obtained in similar investigations with 
other metals can be correlated qualitatively with the 
peculiarities of the tungsten crystal lattice. It seems likely 
that the explanation of the peculiarly governed beams 
observed in the AA’ azimuth of the (1-1-2) plane and of 
Kikuchi’s ‘‘N-pattern’”’ for mica is somehow contained in 
the relation of atomic plane population to interplanar dis- 
tance. 





I. INTRODUCTION 


HE present work was undertaken partly for 

the purpose of trying out a new type of 
apparatus! to see how it would perform in actual 
service, and partly to see whether a change from 
the low melting point face-centered cubic crystals 
which have so far been the usual subject of low- 
speed single crystal work to a high melting point 
body-centered cubic crystal would produce any 
noteworthy difference in the results. Having 
decided upon a new type of apparatus and a 
different type of crystal, conventional procedure 
was further violated by selecting as the first plane 
for study the (1-1-2) plane which has several 
unique characteristics as pointed out in section 
IV. The (1—0—0) plane was studied subsequently. 
The apparatus is unusual in two respects. (1) It 
enables one to analyze diffracted electron beams 
emitted at angles with the normal as small as one 
pleases, down to zero, even though the primaries 
are incident normally. (2) It eliminates the 


1 Sproull, R. S. I. 4, 193 (1933). 


necessity of allowing for contact potentials 
between the filament and other parts of the 
electron gun when determining the velocity of 
the incident primary electrons. 


Il. THe CRYSTAL 


The crystal was cut from a bar of tungsten 
about 20 mm long by 7 mm wide by 5 mm thick 
for which we are indebted to Professor P. I. 
Wold of Union College. The material reached 
Professor Wold rather indirectly from Germany, 
and formed part of a larger bar with which he 
proposes to carry out Hall-effect measurements. 

The bar contained perhaps a dozen crystals 
ranging in size from about 1X3X1 mm up toa 
large one which was approximately 7 X7 X4 mm. 
A section of the bar containing this crystal along 
with two or three other much smaller ones was 
carefully sawed out using a thin sheet of copper 
armored with carborundum powder. The bar was 
moved back and forth along the edge and held in 
the fingers so as not to damage it by clamps. 

The positions of the various atomic planes 
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were determined before the crystal was sawed 
from the bar, taking advantage of the well-known 
fact that chemicals which attack a metal crystal 
dissolve away the surface in such a way as to 
expose facets which coincide with various of the 
more prominent atomic planes. The bar was 
dropped into a boiling 3 percent solution of H2O2 
and etched for about five minutes and then 
mounted on a goniometer. The angular positions 
of the reflected beams of light from an incident 
beam of parallel rays were determined with 
respect to one edge of the bar and the normal to 
the surface. By locating a dozen or so of these 
beams, it is possible to determine from which 
plane each individual beam is reflected, and thus 
to determine accurately the orientation of the 
unit atomic cubes in the bar. It was found that 
the large 7X7 mm face of the bar lacked only 6° 
40’ of coinciding with the (1-1-2) plane. The 
(1-1-2) plane is an important one in a body- 
centered cubic, and is the plane which tends to 
develop most prominently in tungsten when it is 
etched, according to Smithels.2 From this, we 
suspected that it might also tend to become more 
prominent under the action of outgassing in a 
vacuum, 

Next, the crystal was sawed from the bar. A 
steel wedge about an inch square with its faces at 
an angle of 6° 40’ with each other was machined, 
and a hole about 1 cm in diameter was bored 
through it. The wedge was then bolted to a 
steel plate an inch square, the bolt heads being 
countersunk so as to leave the upper face of the 
wedge clear. The upper and lower faces of the 
tungsten segment containing the crystal were 
parallel. Consequently, when it was placed in the 
hole in the wedge with its bottom face resting on 
the steel plate, and then turned to the proper 
azimuth, the part which protruded above the 
upper face of the wedge was just the part which 
should be removed in order to expose a (1-1-2) 
plane. Wood’s metal was poured around the 
crystal, and when it froze, the crystal was ready 
for grinding. It was ground with fine carborundum 
powder and water on plate glass until flush with 
the upper face of the wedge. It was then polished 
to mirror brilliance on a high-speed linen wheel 
with emery flour. After removal from the wedge, 


*See Smithels’ book Tungsten, Chapman and Hall, 
London, 


the crystal was re-etched and mounted in an 
x-ray camera, and from the positions of the 
Bragg x-ray lines it was determined that the 
(1-1-2) plane really coincided with the geo- 
metrical surface within less than a degree. 

For the electron diffraction experiments, the 
crystal was inserted in a molybdenum socket 
having jaws which were too small to be struck by 
any appreciable number of primary electrons, the 
socket being rigidly fastened to a long heavy 
quartz tube mounted in bearings so that the 
azimuth could be varied by rotating the quartz 
tube. The tube was rotated and the molybdenum 
socket bent until no perceptible wobble of the 
(1-1-2) face remained. In this manner the 
(1-1-2) face was aligned with its normal within 
10 or 20 minutes of the axis of the quartz tube, 
for when the error was 1°, the wobble was very 
noticeable. The quartz tube was then aligned 
with its axis parallel to the direction of motion of 
the primary electrons striking the crystal, this 
adjustment being made by heating and bending 
the outer Pyrex walls of the experimental tube, 
and being accurate to }°. 


III. OuTGAssING TREATMENT AND REPRODUCI- 
BILITY OF RESULTS 


A vacuum of 10-7 mm of mercury or better was 
maintained throughout the work, and the crystal 
was outgassed continuously except when a set of 
readings was being taken, the temperature at 
the start being about 1000°C, and being increased 
gradually to about 1600°C after 1550 hours at the 
close of the experiments. After the first 100 hours 
of outgassing, the contact potential between the 
tungsten crystal and the surrounding molyb- 
denum chamber became constant and remained 2 
volts during the rest of the work. This early 
attainment of stability is in agreement with the 
observations of Dowling*® and Glasoe.* After the 
first 300 hours of outgassing, it was found that 
diffraction data could be duplicated almost 
exactly if a given set of observations were 
repeated, even after a lapse of several weeks. 

Langmuir has found that a monatomic layer of 
oxygen clings to tungsten surfaces even at very 
high temperatures, half of it being driven off in 15 


* Dowling, Phys. Rev. 25, 812 (1925). 
4 Glasoe, Phys. Rev. 38, 1490 (1931). 
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seconds at 1800°C. Boas and Rupp’ in their work 
with polycrystalline tungsten found ‘additional 
beams” fitting the usual diffraction theory if one 
postulates a lattice spacing of 4.46A which they 
attributed to this oxygen layer. Such beams were 
not found in the present work, which might 
indicate that the oxygen layer and other im- 
purities were drawn off as positive ions by the 
electric field while the crystal was being bom- 
barded at 2000 volts although they might not 
have been removed by simply heating the crystal 
to the same, or even a higher temperature. 


IV. PECULIARITIES OF THE (1-1-2) PLANE 


According to Davey,® the length of an edge of 
the unit cube in tungsten is d= 3.155A. Fig. lisa 
plan of the (1-1-2) plane. If the circles be taken 
as atoms at the corners of the unit cubes, the 
squares will represent “‘body-centered atoms,”’ 
and vice versa. Those numbered 1 are in the 
surface layer; those numbered 2 are in the first 
layer beneath the surface, those numbered 3 in 
the second, etc. It is to be noted that there is line 
symmetry about the line BC in Fig. 1, but not 


‘ 
;OHOUNOHO 
| — 








-n® GY © © BH o-« 


ro) 





OHOWORO 
1 


Fic. 1. The tungsten lattice viewed along the perpendicular 
to the (1-1-2) planes. 


about the line AA’. Hence electron beams in the 
azimuths marked A and A’ should be alike, but 
beams in the B azimuth should be different from 
those in the C azimuth. Atoms in the sixth layer 
beneath the surface lie directly beneath surface 
atoms, contrasting with the situation in the 


5 Boas and Rupp, Ann. d. Physik [5], 7, 983 (1930). 
* Davey, Phys. Rev. 26, 736 (1925). 
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(1-0-0) plane, where atoms directly beneath the 
surface atoms are found in the second layer 
beneath the surface, and in the case of the 
(1-1-1) plane, the third, these planes being cited 
for comparison because they are the ones other 
experimenters have chosen for study. The atoms 
in the (1-1-2) planes are staggered to an unusual 
extent with respect to the normal to the planes, 
and the perpendicular distance between the 
(1-1-2) planes is only d/6! or 1.29A, very much 
shorter than the corresponding distance for the 
crystal faces usually bombarded by other experi- 
menters with single crystals. Another unusual 
feature of the (1-1-2) plane is that the distance 
measured in the AA’ azimuth between rows of 
atoms parallel to the BC azimuth is 43A, much 
longer than the “grating spacings” that have 
been used in most of the single crystal electron 
diffraction experiments. Consequently, the ratio 
of the grating spacing to the interplanar distance 
is extraordinarily high in this azimuth. Another 
fact possibly of significance is that the rows of 
atoms parallel to the BC azimuth happen to be 
the most densely packed rows in the crystal. 
That is, the atoms in these rows (which are the 
diagonals of the lattice cubes) are 3'd/2 apart, 
and this is the distance of closest approach for 
tungsten. These peculiarities are described at 
length because they may help to account for the 
unusual behavior of the secondary beams 
diffracted in the AA’ azimuth. 


V. THEORY 


Throughout the paper, we are dealing ex- 
clusively with the ‘‘full velocity secondaries,” 
that is, those electrons which have undergone 
elastic collision with atoms of the crystal and 
have recoiled with a velocity equal to that of the 
incident primary electrons. 

In Fig. 1, imagine a plane perpendicular to the 
paper (i.e., perpendicular to the (1-1-2) plane) 
and lying in the azimuth marked AA’. The 
arrangement of the atoms in this section is 
represented in the inset of Fig. 2. By examining 
this inset it is easy to see that phase waves of 
length \ incident normally upon the crystal and 
diffracted by the atoms Z and X at a co-latitude 
§@ as shown will be in phase after diffraction 
provided 

n\\= 2!d sin 6, (1) 
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Fic. 2. Graphs of the interference equations for the (1-1-2) plane in the azimuths A or A’ of 


Fig. 1. Inset is a vertical section of the (1-1-2) planes in the AA’ azimuth. 


where 7; is any integer and d is the edge of a unit 
cube= 3.155A. This is the condition for con- 
structive interference between all the atoms of 
any one layer, by an obvious generalization. 

If, in addition, we are to have constructive 
interference between the rays diffracted from all 
the layers, the waves diffracted from such atoms 
as X and Y must be in phase. That is, NY+ YM 
must be an integral number of wave-lengths, m2, 
whence we have 


nor\= d/6'+2(d/6') sin (30°—@) (2) 


as the second interference condition. Beams are 
usually expected only when \ and @ have such 
values as to satisfy both (1) and (2) simultane- 
ously. Eq. (1) is often called the “surface 
condition” and (2) the ‘‘volume condition”’ of 
interference. Following usual practice, these two 
equations are graphed (main part of Fig. 2) with 
sin @ as ordinates and X as abscissae, the inter- 
sections of the two sets of curves predicting the 
co-latitude @ and the wave-length \ (and hence 
the primary voltage V, since \= h/mv= (150/ V)!) 


at which the beams should occur. The straight 
lines represent the surface condition for m,=1, 2, 
etc., and the curved lines represent the volume 
condition for m2=1, 2, etc. 


OD) — ated) 
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x x Ka 
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x 
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Q-1-) (-@0) (1-1-1) 


Fic. 3. The tungsten lattice viewed along the perpendicular 
to the (1-0-0) planes. 


If we do the same thing for the (1-0-0) plane 
in the azimuths marked (1-0-0) (Fig. 3), we 
obtain the surface condition 


n,\=d sin 0 
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and the volume condition 
n2r\= $d+d/2! sin (45°—8@) 


illustrated and plotted in Fig. 4. 
Corresponding equations for some of the other 
cases studied are: (1-1-2) plane, B azimuth 


mA= 3'd/2 sin 0, (surface) 

n2= (2/3)'d+3'd/2 sin (70° 30’—8@); (volume) 
(1-1-2) plane, C azimuth 

n\= 3'd/2 sin 6, (surface) 

nov= d/6'+3d sin (54° 45’—@). (volume) 
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Fic. 4. Graphs of the interference equations for the + 


(1-0-0) plane in the (1-0-0) azimuths of Fig. 3. Inset is 
a — section of the (1-0-0) planes in the (1-0-0) 
azimuth. 


These sets of equations differ because there is not 
line symmetry about the line AA’ (see Section 
IV). For the case of diffraction in the azimuths 
marked X due to the rows of atoms indicated by 
the dotted lines in Fig. 5, one obtains 


m\= (6/11)'d sin 6 (surface) 

and 

n2d= d/6'+ (2/11)'d sin (73° 10’—@); (volume) 

(1-0-0) plane, (1-1-1) azimuths (see Fig. 3) 
n,\\= 21d sin 6, (surface) 
n2d\= 3d+3'd/2 sin (35° 16’—@). (volume) 





Fic. 5. Plan of the (1-1-2) plane similar to Fig. 1, but 
showing only one layer. 





SPROULL 


In the present apparatus, the collector can be 
moved into such a position as to catch secondary 
electrons leaving the crystal at normal exodence, 
so that the case 6= 0 can be studied. One should 
note that this represents the zeroth order of 
surface interference, that is, ™,=0 for all the 
surface equations, which degenerate to 0=0. In 
Figs. 2 and 4 this means that the surface equation 
is represented merely by the \ axis which should 
therefore be regarded as one of the system of 
straight lines in the graphs. At @=0, therefore, 
one should expect to find maxima in the collector 
current at certain specific values of \ (and hence 
V) determined by the intersections of the curved 
lines with the \ axis. When @= 0, all of the volume 
equations for the various azimuths of any one 
plane become identical, as one would expect. 


VI. THE OBSERVED BEAMS 


Each set of observed beams will be described 
upon the basis of the following aspects: 

A, Intensity of the beams compared to the 
background of random scattering upon which 
they are superimposed. 

B, Sharpness in co-latitude. That is, main- 
taining the bombarding potential constant at the 
value at which a typical beam of the set is fully 
developed, does the electron current at various 
colatitudes in the neighborhood of the maximum 
decrease rapidly or slowly for a given azimuth? 

C, Sharpness in azimuth. Maintaining the 
bombarding potential constant at the value at 
which a typical beam of the set is fully developed, 
does the electron current at various azimuths in 
the neighborhood of the one being studied 
decrease rapidly or slowly for a given co-latitude? 

D, Rate of growth and decay. If one explores a 
typical beam of the set by varying the co- 
latitude with constant azimuth, not only at the 
critical bombarding potential at which the beam 
is fully developed, but also at various other 
values of V in the neighborhood, does the beam 
fade out or “‘decay”’ rapidly as V recedes from 
this critical value? 

E, Mode of growth and decay. If one maintains 
the bombarding potential at a constant value 
differing slightly from the critical value at which 
a typical beam of the set is fully developed, is the 
location of the beam in co-latitude such as to 
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agree closely with the surface interference con- 
dition or with the volume interference condition, 
or with neither? (A beam should be fully 
developed in a given azimuth only when both @ 
and V are chosen so as to simultaneously satisfy 
both conditions; when V is varied slightly, both 
conditions can not be satisfied.) 


1. Sets of beams in the B and C azimuths of the 
(1-1-2) plane. (Fig. 1) 

A, Very weak and difficult to detect. B, 
Moderately sharp in co-latitude. C, Moderately 
sharp in azimuth. D, Very rapid decay. E, Mode 
of decay indeterminate, because decay is too 
rapid. 


2. The sets of beams in the azimuths of the 
(1-1-2) plane marked X in Fig. 5 


Exactly similar to the sets just described. 


3. Beams in the (1-1-1) azimuths of the (1-0-0) 
plane (Fig. 3) 

A, Weak. B, Very broad in co-latitude. C, 
Broad in azimuth. D, Fairly rapid decay. E, 
Tendency to satisfy the volume condition rather 
than the surface condition. 


4. Beams in the (1-0-0) azimuths of the (1-0-0) 
plane (Fig. 3) 

A, Rather weak. B, Very broad in co-latitude. 
The ,;=1=m2 beam is shown plotted in co- 
latitude (@) in Fig. 6A where the radii represent 
collector current for the various values of 6. The 
arrow indicates the peak which is approximately 
at its maximum development at this voltage. 
The 2,;=1; m2=2 beam is shown similarly in 
Fig. 6C. C, Broad in azimuth. The »,;=1=n,2 
beam is shown plotted in azimuth in Fig. 7, the 
radii again representing collector current. D, 
Unusually slow growth and decay. The n,;= 1= m2 























Fic. 6. Co-latitude curves for the (1-0-0) plane. A, 
The nm,=1=nz beam at full development. B, Two curves 
showing the same beam at a bombarding potential 10 
volts less than A.C, The m;=1; mg=2 beam. 























Fic. 7. An azimuth curve of the m,=1=m, beam for 
the (1-0-0) plane taken at a co-latitude of 30°; bombarding 
potential V =23.8 volts. 


beam could be detected at any value of V 
between 25 and 50 volts. Fig. 6B shows co- 
latitude curves taken at a value of V ten volts 
less than that at which the beam is fully de- 
veloped, the peaks being indicated by the 
arrows. These curves were both taken in (1-0-0) 
azimuths, but are 90° apart in azimuth (see 
Fig. 3) and were taken about two weeks apart. 
The curves of Figs. 6A, B, C were all corrected for 
change of solid angle subtended by the collector 
opening as discussed in the article in R. S. I.! E, 
Close agreement with the volume condition. 


5. Beams in the A and A’ azimuths of the (1-1-2) 
plane (Fig. 1) 

A, Very strong. B, Exceedingly sharp in co- 
latitude; see Fig. 8. C, Very sharp in azimuth; see 
Fig. 9, note curve marked 22° (co-latitude). D, 
No decay whatever; the various co-latitude 
curves of Fig. 8 were taken at the various 
voltages indicated, and it is seen that the beam 
is equally well developed at every voltage which 
was tried. Only the central portion of each peak 
was plotted because it was desired to record the 
whole set in one afternoon under exactly the 
same conditions. These peaks are so narrow that 
it was not necessary to make correction for 
change of solid angle subtended by the collector 
opening. E, After allowance for refraction (to be 
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Fic. 8. Co-latitude curves showing the unusual series of m2: =1 maxima found in the AA’ azimuth 
of the (1-1-2) plane at ten different values of V as indicated. 























Fic. 9. Azimuth curves of the 2=1 maxima in the AA’ 
azimuth of the (1-1-2) plane at a bombarding potential 
of 44.2 volts at co-latitudes of 20, 22 and 25 degrees as 
indicated 


discussed in Section VIII), the positions of the 
beams agree almost perfectly with the volume 
condition (2), and since there is no voltage of 
maximum development, it seems that these beams 
are governed exclusively by Eq. (2) and not at all 
by Eq. (1). 


VII. CoRRELATION OF THE OBSERVED BEAMS 
WITH THE THEORY 


The sets of beams described in the preceding 
section were enumerated in order of their agree- 
ment with the theory of Section V. The first 
beams described behaved most nearly in accord 
with the theory, while the ones described last 
departed markedly from it. The behavior of the 
beams observed in the (1-0-0) azimuth of the 





(1-0-0) plane at voltages in the vicinity of the 
critical value of full development can be illus- 
trated graphically in Fig. 4 by saying that they 
could be observed throughout the regions X and 
Y indicated by the dotted lines near the inter- 
sections of the curves. The tendency to obey the 
volume condition rather than the surface con- 
dition is indicated by the fact that these regions 
lie along the curved (volume) graphs and not 
along the straight (surface) graphs. Other experi- 
menters’ studying the process of growth and 
decay of the electron beams have found a 
tendency to obey the surface condition more 
closely than the volume condition. The circles 
numbered 1, 2 and 3 in Fig. 4 mark the values of 6 
and \ at which the curves of Fig. 6A, B and C 
were taken respectively. 

The beams found in the AA’ azimuth of the 
(1-1-2) plane and shown in Figs. 8, A-J are 
located at values of @ and X indicated by the 
circles in Fig 2, where one sees immediately that 
they are falling along the »2=1 volume curve.® 
They do not fall exactly on the curve because we 
have so far neglected to correct for refraction. 
This correction will be made in the next section. 


7See Davisson and Germer, Phys. Rev. 30, last para- 
graph, p. 734 (1927) or Farnsworth, Phys. Rev. 34, top 
of p. 693 (1929). 

8 The striking behavior of these beams as regards growth 
and decay can be illustrated graphically by comparing 
Fig. 2 with Fig. 3 obtained by Farnsworth for copper, 
Phys. Rev. 34, p. 685 (1929), or comparing Fig. 8 with 
Fig. 10 obtained by Davisson and Germer for nickel, 
Phys. Rev. 30, 716 (1927). 
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The cross in Fig. 2 indicates a beam found after 
only 28 hours of outgassing which evidently falls 
in line with the circled beams which were dis- 
covered later after 440 hours of outgassing at 
about 1200°C. The beams were followed right up 
to d=0. 


VIII. CoRRECTION FOR REFRACTION AND CAL- 
CULATION OF THE SURFACE WORK 
FuNcTIONS W, AND W; 


A. Normal exodence 


In the case of the experiments at 6=0, one 
does not need to make any correction for bending 
at the surface. The only effect of refraction is to 
make the length of the phase waves within the 
crystal less than it is outside. Therefore one can 
merely subtract the experimental value of the 
bombarding potential V at which a maximum 
occurs from that at which it is predicted and say 
that the difference represents the surface work 
function W,. If one sets @= 0 in Eq. (2) and solves 
for \ and then for V using A= (150/ V)', it turns 
out that maxima in the collector current should 
be expected at 6=0 for the (1-1-2) plane for the 
values of V given in the second column of Table 
I, corresponding to the choices of m2 in the first 
column. Proceeding similarly for the (1—0—0) 


TABLE I. (1-1-2) plane. 


























Veal Vobs. Wa( = Veale. —_ Vobe.) 
? (volts) (volts) (volts) 
1 22.5 16.0 (218 hrs. outgassing) 6.5 
2 90.0 84.6 (259 hrs. outgassing) 5.4 
1 22.5 19.25 (383 hrs. outgassing) 3.25 
2 90.0 84.6 (383 hrs. outgassing) 5.4 
TABLE II. (1-0-0) plane. 
n Veale. obs. W.( rad Veale. _ obs.) 
. (volts) (volts) (volts) 
1 15.1 7.9 7.2 
2 60.2 56.8 3.4 








plane, we obtain Table II. The average of the six 
values of W, in Tables I and II is 5.19 volts. The 
variation in the six values may be due to small 
Variations in the angle of incidence as pointed out 
by Farnsworth.® The 7.9 volt maximum of Table 





* Farnsworth, Phys. Rev. 40, 684 (1932). 
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Fic. 10. The m;=0; m2=1 (normal exodence) maximum in 
the (1-0-0) azimuth of the (1-0-0) plane."® 


II is shown plotted" in Fig. 10 which shows some 
indication of an accompanying satellite of the 
type found for copper and silver by Farnsworth.® 
The positions of these 2= 0 maxima are indicated 
in Figs. 2 and 4 by triangles. 


B. Oblique exodence 

The extreme sharpness and great intensity of 
the beams found in the A and A’ azimuths of the 
(1-1-2) plane made it seem advisable to base the 
estimates of W, exclusively on these beams, for 
all of the other beams were so broad in com- 
parison that estimates of comparable accuracy 
would be impossible. In Fig. 2, it is supposed that 
the points lie near the volume curves instead of 
on them because no correction has been made for 
refraction. If the curves are now corrected for 
refraction until they are shifted over to coincide 
with the points, one can estimate from the 
amount of the correction necessary what must 
have been the value of the surface work function 
W, to which the refraction is due. We now 


proceed to do this. 


10 Tn Fig. 10, the abscissae give the apparent bombarding 
potential as read by a voltmeter (connected from filament 
to grid) and it is seen that the voltage of the maximum 
read in this way is about 15 }, although the true bombarding 
potential was 7.9 volts as marked, the discrepancy being 
due to contact potentials in the electron gun. These contact 
potentials which are not only large, but quite variable 
from day to day are of no importance in the present work, 
but Fig. 10 was plotted so as to emphasize them because 
they must be considered in using apparatus of the usual 
type, and their wide fluctuation may account for the 
irregular variations of W, reported by some observers. 
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It must be remembered that due to refraction, 
the length of the phase waves within the crystal 
is not \, but \’, and that the angle @ at which the 
beams are observed is not identical with the 
angle @’ at which the diffraction is really oc- 
curring. Consequently, Eq. (2) should be cor- 
rected by replacing \ by ’ and @ by @ using the 
relation 

\/d\’ =y=sin 6/sin 6’ (3) 


mw being the index of refraction given by yu 
=(1+W./V)' where W, is the surface work 
function. If Eq. (2) is thus corrected, and if at the 
same time one substitutes for d its value (3.155A), 
one obtains for the first order equation (m2= 1) 


\’=1.29(1+2 sin (30°—90’)) (4) 


or substituting for \’ and 6’ their values in terms 
of X, 0, and w from Eq. (3), one obtains 


d sin? @\ } sin 6 
= =1.29(1+(1- ) ~1.732 ). (5) 
im we m 


Solving this equation for yu, one obtains 
(A/1.29)?+2.685A sin 6+4 sin? 0 
~——-9(A/1.29+1.732 sin 8) 








(6) 





In Table III, the first column gives the co- 
latitude at which the beam occurred as read 
from Fig. 8. The second column gives the 
primary voltage V at which the beam was 


TABLE III. A—A’ azimuth (1-1-2) plane. 











6 V r Wa, 
(degrees) (volts) (Angstroms) B (volts) 

0 17.62 2.918 1.1324 5.0 
12 27.55 2.334 1.0949 5.5 
19 37.75 1.994 1.0773 6.1 
22 44.2 1.843 1.0743 6.8 
23 47.45 1.779 1.0659 6.5 
243 53.05 1.682 1.0544 5.9 
253 $7.75 1.612 1.0448 . 
27 63.35 1.539 1.0421 5.4 
28 68.1 1.485 1.0396 a5 
31 83.0 1.535 1.0318 5.5 
35 118.5 1.125 1.0198 4.7 

Mean: 5.65 








actually observed, this being computed from the 
magnetic field as explained in the article in R. S. 
I.! The wave-length \ corresponding to this 
primary voltage as computed from the equation 
A=(150/V)! is given in the third column. 
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Substituting the values of @ and \ from columns 
1 and 3 into Eq. (6), one computes the values of 
mw given in the fourth column. Then using the 
equation »=(1+W./V)! one computes the 
values of W, recorded in the fifth column. These 
values are seen to be fairly consistent and their 
average is 5.65 volts. If one includes in the 
average the values for W, given in Tables I and 
II, one obtains a final average for W, of 5.52 
volts, which the writer regards as his best 
experimental value. 


IX. EFFECT OF VARIATION OF ANGLE 
OF INCIDENCE 


In a recent paper,? Farnsworth emphasizes the 
importance of slight variations in the angle of 
incidence. An azimuth curve was taken with the 
normal to the (1—0—0) plane inclined at an angle 
of 2 or 3 degrees to the direction of motion of the 
incident primary electrons at an angle @ and 
voltage indicated by the star in Fig. 4. The curve 
obtained is the solid one in Fig. 11. The tube was 























Fic. 11. Azimuth curves showing the effect of slight 


variations in the angle of incidence. 


then cut down, the crystal re-aligned, followed by 
reassembly of the tube, the usual two weeks of 
baking, and 100 more hours outgassing of the 
crystal. The azimuth curve was then repeated at 
the same latitude and voltage, the result being 
shown by the dotted curve in Fig. 11. The solid 
curve exhibits peaks very strongly in the (1-0-0) 
azimuth at 90° and 180°, and weakly at 270° and 
0°, but the dotted curve exhibits no peaks at all. 
The importance of small variations in the angle of 
incidence is thus again demonstrated. 
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X. SUBSEQUENT EXAMINATION OF THE CRYSTAL 


Microscopic examination of the crystal after 
the work was completed revealed slight pitting 
of the surface by the bombardment, but the 
crystal facets flashed up in unison over the 
entire face even more brightly than they had 
just after the original etching, indicating that if 
any recrystallization had occurred, it was ultra- 
microscopic. 


XI. Discussion 


A. General 


In the following discussion, the observed 
peculiarities of electron diffraction in tungsten 
are enumerated, followed in each case by a brief 
attempt to partly correlate them with the known 
structure of the tungsten lattice. 

(1) Except in the case of the AA’ azimuth of 
the (1-1-2) plane, the beams on the whole are 
less intense, less sharp in co-latitude and 
azimuth, and have a slower rate of growth and 
decay (these expressions were explained in Sec- 
tion VI) than those observed by others for 
nickel, copper, silver, etc. 

The (1-0-0) and (1-1-2) planes of tungsten 
are somewhat less densely populated than the 
planes examined in these other metals. For 
example, the populations of the tungsten (1—1-—2) 
and (1-0-0) planes, the (1-1-1) plane of nickel, 
and the (1-0-0) planes of silver and copper are 
approximately in the ratio 81 : 99 : 106 : 121 : 155 
respectively. This circumstance might enable the 
primary electrons to penetrate more deeply into 
the crystal, and the secondaries originating deep 
within the metal might suffer a second diffraction 
or undergo irregular scattering on their way back 
to the surface. Such successive diffraction might 
help to account for the generally lessened sharp- 
ness and rate of growth and decay and for the 
increased diffusely scattered background. 

(2) In general, the volume interference con- 
dition seems to play a more important rdle in 
tungsten than in other metals. This also could be 
explained on the basis of deeper penetration, 
which would tend to increase the intensity of the 
volume maxima and also to increase their 
sharpness due to improved resolving power. 

(3) A unique behavior is observed in the case of 
the beams diffracted in the AA’ azimuth when 
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the primary beam is incident normally upon the 
(1-1-2) plane. Here the beams were sharper in 
co-latitude and azimuth than are usually ob- 
served for other materials, and were found at 
every voltage tried, always at a value of @ 
agreeing with the volume Eq. (2). In this 
azimuth of this plane, the surface interference 
Eq. (1) has negligible effect, or none at all. 

The sparse population of the (1-1-2) planes 
coupled with their small interplanar distance and 
the unusual staggering of the atoms about the 
normal pointed out in Section IV might allow 
penetration of an unusual number of planes by 
the primaries. From Fig. 1, one can see that the 
secondaries diffracted at oblique exodence in the 
B or C azimuths will encounter many obstructing 
atoms in their journey back to the surface with 
the results already mentioned in Part 1 of the 
discussion, but secondaries diffracted at oblique 
exodence in the A or A’ azimuth will find a 
comparatively free path back to the surface if @ 
is less than about 35°. This makes the extreme 
sharpness, great intensity, and the accentuation 
of the volume interference somewhat more 
plausible, but fails to account for the absence of 
destructive interference when Eq. (1) is not 
satisfied. 

The unique behavior is observed with a 
sparsely populated closely spaced set of planes. 
The opposite extreme in behavior was observed 
by Kikuchi with the opposite extreme in struc- 
ture, though with primary velocities from 10,000 
to 85,000 electron-volts. Bombarding films of 
mica about 10-> cm thick normally to the 
cleavage planes which are spaced at the relatively 
great distance of 10A"™ and recording the trans- 
mitted secondaries by a photographic plate 
behind the film, he obtained” a so-called ‘‘N- 
pattern” which fits perfectly the surface inter- 
ference condition applied to the cleavage plane, 
while the volume condition had no inftuence 
which could be even detected. The same phe- 
nomenon was noted in mica with x-rays by W. 
Linnik.* 


It seems worth noting that while with tungsten 


the breakdown of the surface condition for 


1 Siegbahn, Phys. Rev. 8, 320 (1916). 

1 Kikuchi, Proc. Imp. Jap. Acad. 4, 271, 275, 354, 471 
(1928); Japanese Journal of Physics 5, 83 (1928). 

3 Linnik, Nature 123, 604 (1929). 
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constructive interference is associated with inci- 
dence on relatively empty but closely spaced 
planes, Kikuchi’s observed breakdown of the 
volume condition for constructive interference is 
associated with relatively full planes which are 
widely spaced. This suggests that the explanation 
of both anomalies is somehow contained in the 
relation of population to interplanar distance. 

This problem has already aroused the interest 
of others. W. L. Bragg'* suggested that the 
phenomenon could be “readily explained by the 
familiar laws of diffraction by a three dimensional 
grating, assuming a slight random warping of the 
mica planes,” and performed an experiment to 
demonstrate this possibility. It seems doubtful 
whether the present results can be explained on 
the basis of an opposite type of deformation of the 
tungsten crystal. 

The sharpness of the beams in the AA’ azimuth 
of the (1-1-2) plane enables one to compute that 
the phase waves assumed to be associated with 
the electrons have an effective lateral extent at 
right angles to the direction of propagation of at 
least 30A if one is permitted to apply the usual 
formulae for resolving power to the phase waves. 
A similar conclusion was reached by G. P. 
Thomson.'§ 


B. Values for W, 


According to Sommerfeld’s theory of metals,'® 
we have 


o= W.-W; (7) 


4 Bragg, Nature 124, 125 (1929). 

16 See Thomson’s book, Wave Mechanics of Free Electrons, 
p. 72. 

16 Sommerfeld, Zeits. f. Physik 47, 31 (1928). 
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where ¢ is the usual photoelectric or thermionic 
work function (4.55 volts for tungsten), W, is the 
height of the potential wall at the metal surface, 
and W; is the maximum internal energy of the 
free electrons in the metal. Therefore, if we take 
the final value of 5.52 volts for W, yielded by the 
present experiments, we get W;=about 1 volt, 
whereas an application of Sommerfeld’s theory to 
tungsten yields a theoretical value of 5.7 volts, 

The only other determination of W, for 
tungsten is that of Boas and Rupp® who obtained 
values of 0 and 10 volts for polycrystalline 
material. As regards consistency, the present 
results are considerably more satisfactory than 
those recently given by Farnsworth’ for a single 
copper crystal in which case the values of W, 
range from 5 to 31 volts. The fact that higher 
outgassing temperatures could be used with 
tungsten may account for the more consistent 
values obtained, and this consistency leads one to 
attach some possible importance to the difference 
between the resulting value of W; and that 
computed from the Sommerfeld theory. This 
difference may be due to the limitations of the 
Sommerfeld theory resulting from the fact that 
the variations in potential inside the crystal due 
to the atoms is not considered; or perhaps one is 
not justified in treating the diffraction from the 
crystal planes in such a highly idealized way, as 
though the problem were quite analogous to 
diffraction of light in optics; or it may be that 
unsuspected systematic experimental errors are 
responsible. 

In conclusion, the writer wishes to express his 
sincere appreciation of the invaluable help so 
willingly offered by Professor C. E. Mendenhall, 
who directed the work. 

















APRIL 1, 1933 


PHYSICAL REVIEW 


VOLUME 43 


The Fluorescence Yield from the Z;;, Level of Uranium 


REGINALD J. STEPHENSON, Ryerson Physical Laboratory, University of Chicago 
(Received February 3, 1933) 


The fluorescence yield, defined as the ratio of the number 
of fluorescence quanta emitted from the Li11 level to the 
total number of incident quanta absorbed in the Lr level 
has been measured for uranium; and the value of 0.67 
found. The relative intensities of the a to the #-lines in the 
K fluorescence spectra of molybdenum have been measured 
and found to agree with the value from the characteristic 
spectrum. The relative intensities of the a to the @-lines 


originating from the Li11 level have been measured in 
the fluorescence spectrum of uranium, and found to agree 
with the value from the characteristic spectrum. The 
atomic absorption coefficients for uranium for wave-lengths 
including the Zi11 discontinuity have been measured, 
and a value for 6=2.27 found, where 6 is the ratio of the 
absorption coefficients on the short and long wave-length 
sides of the Lr11 discontinuity. 





CCORDING to the theory of the production 

of fluorescence x-rays, when an electron in 

the K energy level absorbs a quantum of x-rays, 
it is ejected as a photoelectron. An electron from 
the L or outer levels will then fall into the vacant 
K level with the emission of K fluorescence radia- 
tion. Thus for each quantum of energy absorbed 
in the K level there should be emitted one 
quantum of fluorescence K radiation. However, 
it has been apparent for some years that this 
simple picture does not give a complete account 
of the phenomena associated with the reorganiza- 
tion of the atom. In some early experiments on 
the efficiency of production of fluorescence 
x-rays, excited by K characteristic radiation, for 
atoms belonging to the iron group Barkla and 
Sadler! found that only some 30 percent of these 
atoms ionized in the K level emitted K fluores- 
cence radiation. Barkla? and Beatty’ also found 
that the excitation of the K series was accom- 
panied by a sudden increase in the emission of 
the photoelectrons. Auger,‘ using a Wilson 
cloud-chamber with different inert gases, was 
able to photograph the tracks of the photoelec- 
trons and measure their energy. He found that in 
addition to the photoelectron whose energy is 
given by Einstein’s equation Eyin=hv—hyx for 
absorption by a K electron; there might be 


1 Barkla, Bakerian Lecture, Phil. Trans. 217, 317 (1917). 

? Barkla and Philpot, Phil. Mag. 25, 489 (1923). 

* Beatty, Proc. Roy. Soc. A85, 329 (1911). 

‘ Auger, C. R. 177, 169 (1923); C. R. 180, 65 (1925); 
Ann. de Physique 6, 183 (1926). 


emitted from the same atom an electron of 
energy approximately E=hvx—2hv, and others 
of a similar type. This simultaneous ejection of 
several 8-particles from the same atom has been 
called the “compound photoelectric effect.”” By 
knowing the energy and number of the £- 
particles, Auger was able to count the number of 
quanta absorbed in the K level, and the number 
of fluorescence K quanta emitted. He defined the 
fluorescence yield for the K level as the ratio of 
the number of fluorescence K quanta emitted to 
the total number of incident quanta absorbed in 
the K level. A similar definition may be applied 
to the L levels. Using the symbols wx and wy, to 
represent the fluorescence yield from the K and L 
levels, respectively, Auger found the following 
results: 


Auger’s Results on Fluorescence Yield in K and L Series. 


Atomic No. Gas WK WL 
19 argon 0.07 
36 krypton 0.51 0.1 
54 xenon 0.71 0.25 


Kossel® has suggested that the low efficiency of 
the K series emission may be explained on the as- 
sumption that the majority of such atoms un- 
dergo radiationless reorganizations of the Klein- 
Rosseland type resulting in the appearance of 
high-speed photoelectrons. Bothe,® on the other 
hand, has suggested that the results of Barkla 


5 Kossel, Zeits. f. Physik 19, 333 (1923). 
* Bothe, Phys. Zeits. 26, 410 (1925). 
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and Auger can be correlated on the assumption 
that they are manifestations of the phenomena 
of “internal absorption.’’ According to this pic- 
ture K radiation is always emitted in the re- 
organization of a K ionized atom, but the prob- 
ability that this radiation shall be absorbed in the 
L, M, N, etc., levels of the same atom is very 
much greater than that it shall be absorbed in the 
same level of some other atom. In other words, 
the absorption coefficient for the K series radia- 
tion in the L, M, N, etc., levels of the ‘‘parent”’ 
atom is very much greater than for radiation of 
the same wave-length emitted from some other 
atom. There appears, however, to be no method 
of deciding between the two theories from an ex- 
perimental standpoint. The excitation of the K 
series radiations in a solid or gas will be accom- 
panied by a sudden increase in the photoelectric 
emission as a result of the ‘internal absorption”’ 
of the fluorescence radiations in the outer levels 
of the atoms. This was observed by Barkla? and 
Beatty’ in the ionization produced relative to air, 
in methyl bromide, methyl iodide, and selenium 
hydride, when the wave-length of the radiation 
was slightly shorter than the K limit of the gas in 
question. Martin’ has calculated the value for 
wx from these data for bromine, iodine and 
selenium, and found it to give reasonable results. 
The importance of this effect in the measurement 
of relative intensities of x-ray lines with an 
ionization chamber has been well shown by 
Allison and Andrew.® 

Measurements on the fluorescence yield in the 
K series have been made by Compton’ and 
Martin,’ where it was shown that wx increases 
with the atomic number of the element, but does 
not vary with the wave-length of the incident 
x-rays. The latter conclusion involves the as- 
sumption, in the calculation of the energy of the 
x-rays, that the energy required to produce a pair 
of ions in the gas of the ionization chamber is 
constant and independent of the wave-length of 
the x-rays. This assumption appears to be fully 
justified by quite independent experiments.'° 

Since there are three L levels it was decided to 
make a determination of the fluorescence yield 

7 Martin, Proc. Roy. Soc. A115, 420 (1927). 

8 Allison and Andrew, Phys. Rev. 38, 441 (1931). 


* A. H. Compton, Phil. Mag. 8, 961 (1929). 
10 Crowther and Bond, Phil. Mag. 8, 961 (1929). 
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for a single one of these; and accordingly the Ly 
level of uranium was chosen. The method used in 
this experiment is essentially the same as that 
adopted by Compton.’ 
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Fic. 1. Diagram of apparatus used for measuring the 
fluorescence yield from the Lin level of uranium. 


The wave-lengths of the three L levels of 
uranium are given by Siegbahn" as: L;= 0.5680A; 
Ly=0.5913A; Lyy17=0.7208A; and the wave- 
length of the molybdenum Ka,=0.7078A; Mo 
KB,=0.6309A;; so that if this radiation from Mo 
is incident on a block of U only the lines arising 
from the Ly; level will be excited. Primary x-rays 
from a rhodium target fell on a block of molyb- 
denum at R; producing nearly homogeneous K 
radiation. A rhodium target was used in the x-ray 
tube because the wave-lengths of the Rh Ka; and 
RhK&8, are just a little shorter than the K absorp- 
tion limit of Mo, and consequently the fluores- 
cence lines from Mo will be excited quite strongly. 
Fluorescence radiation was used in preference to 
characteristic, because the former is entirely line 
radiation. 

The fluorescence Mo radiation after passing 
through collimating slits fell on a piece of U at 
R2, and the resulting radiation measured by an 
ionization chamber placed in position B. The 
intensity of the direct beam was measured by the 
ionization chamber in the position A. 


1 Siegbahn, Spektroskopie der Réntgenstrahlen, Second 
Edition, p. 274. 
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CALCULATION OF THE FLUORESCENCE YIELD 


Let pw’ be the total absorption coefficient of the 
primary Mo x-rays in the U radiator, then 


p’=r'+o’ 


where 7’ is the photoelectric absorption, and o’ 
the absorption due to scattering. In this case the 
greater part of the photoelectric absorption will 
be due to the electrons in the Ly level, i.e., 71’. 
For the sake of convenience the subscript III 
will be dropped, and it is to be understood that 
7, and w_, refer only to the Lyy; level of U and not 
to all of the levels. The Mo fluorescence beam 
was taken as composed of four lines; a, a, 81, Be, 
and their relative intensities taken from Meyer’s” 
data, namely: a; : ae : Bi : Be=100 : 50.6 : 23.25 
: 3.48. Each of these lines will make its own con- 
tribution to the fluorescence beam from the 
uranium. Consider first the contribution made by 
the a;-line; its intensity is 100/177.3 of the total 
intensity from molybdenum. The fluorescence 
beam coming from the uranium is composed of 
seven lines: a1, a2, Be, Bs, Bs, Bz, 1, and each one of 
these will make its separate contribution to the 
ionization current derived from the uranium. 
Allison’s values for the relative intensities in the 
Lin characteristic spectrum from U _ give 
a, : a2: Bo: Bs :B6 : Bz :1=100: 11:28 :6.4:1.6 
:0.4 : 2.4. Let us therefore first obtain an expres- 
sion for the power entering the ionization cham- 
ber due only to the U La; and excited by the 
Mo Ka. The wave-lengths and relative intensi- 
ties of the Mo and U lines were taken from the 
characteristic spectrum, but as is shown later 
this appears to be quite justified. 

Suppose a beam of x-rays from the Mo of in- 
tensity 7’ and cross section A’ traverses a thick- 





pe A” tie v’” 100 
ieee TL 


100 





where 100/149.8 is the fractional intensity of the 
U La; to the total U beam )’ and X” are the 
wave-lengths of the Mo Ka; and U La, re- 


® Meyer, Wissenschaftliche Verdffentlichungen aus dem 
Siemens-Konzen VII. 2 (1929). 
Also, Siegbahn, Spektroskopie der Réntgenstrahlen, 
Second Edition. 


[conde = Wy — 
v’ 149.8 177.3 YJ Agr? yp’ +p” 149.8 177.3 dX” 
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ness dx of the U radiator. The number of quanta 
absorbed by the L level and hence the number of 
photoelectrons ejected from the L level will be 


dn= (I'A'dx/hv’)(100/177.3)r1’, (1) 


where »’ is the frequency of the Mo Ka;. Hence 
the number of fluorescence L quanta escaping 
from the U atoms, due to the Mo Ka, will be 
dn,=wy,dn. The power in the fluorescence beam 
due to the U La, will be hv’’dn;, where v” is the 
frequency of the U La;. Barkla and Sadler have 
shown that the intensity of the fluorescence beam 
is the same in all directions. Thus the power in 
the U fluorescence beam, due to the Mo Ka; and 
U Lay, uncorrected for absorption; which enters 
the ionization chamber is 


a”? 


dP" =hv"dn,y — rz’ 
4rr 


1 A” 100 
=hv"w,P'r,'dx ———- ——.._ (2) 
hv’ 4rr? 177.3 


A” is the area of the diaphragm S, and P’=I/'A’ 
is the power in the primary Mo beam, and r the 
distance from the center of the U block to the 
center of the diaphragm S,. A’’/47r’ is the solid 
angle which the ionization chamber subtends at 
the uranium radiator. 

If the uranium is set at 45° to the primary 
molybdenum beam and thus at 45° to the 
ionization chamber in the position B, then the 
paths of the primary and the U fluorescence beam 
in the U radiator will be the same. Let yp’ and yu” 
be the absorption coefficients of the Mo Ka; 
and U La, respectively, in the uranium. Then 
the power in the fluorescence beam entering the 
ionization chamber will be 

A" Tr! 


100 100 





“sy (3) 





spectively. Similar expressions will give the con- 
tribution from the U Las, LB, LBs, LBs, LB:, 1; 
these added together give the total contribution 
due to the Mo Ka. To this must be added the 
total contribution due to the Mo Kae, Ki, KBe, 
and then after making the appropriate correc- 
tions for the absorption of each of these lines in 
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the gas of the ionization chamber, the sum total 
will give values proportional to the observed 
ionization currents. 

The absorption coefficients of these various 
wave-lengths in uranium need to be determined, 
as well as the value of r,. It is also necessary to 
test whether the relative intensities of lines of 
the Mo characteristic spectrum agree with those 
from the fluorescence spectrum. This must also 
be checked for the lines from the Zy1; level of 
uranium. 


ABSORPTION COEFFICIENTS IN URANIUM 


Because of the high absorption or uranium it is 
impractical to use foil, even if it were available; 
consequently a solution of uranyl nitrate was em- 
ployed. A parallel sided cell was made 0.320 cm 
thick, the front and back faces covered with mica 
and held in position with rubber gaskets. A single 
crystal ionization spectrometer was used, the cell 
being placed between the two defining slits. A 
motor generator set supplied current at 540 
cycles to the primary of a transformer; the 
secondary current being rectified by kenotrons 
and then passed into condensers to smooth out 
the ripple. The high voltage was read on an elec- 
trostatic voltmeter and kept constant during the 
experiment. The current through the x-ray tube 
was read on a milliammeter and also kept con- 
stant, so that the power through the tube re- 
mained the same during an experiment. This 
equipment was used in all the experiments de- 
scribed in this paper. 

Solutions of 2, 5 and 10 g of uranyl nitrate 
crystals dissolved in 100 g of water were used in 
the cell and the three served as checks. Since the 
absorption coefficient of water is also found this 
served as an additional check, the results agree- 
ing very well with the published data." 

The results shown in Table I were obtained for 
the atomic absorption coefficient for uranium. 
The value of the Ly1; discontinuity, or the ratio of 
ue on the short and long wave-length sides of the 
Li absorption limit is found from the graph to 
equal 2.27. Kustner™ gives the value of 2.2 for 
the Ly discontinuity. From these results it 
would appear that the wave-length of the Zi 

13 Dershem, Phys. Rev. 21, 30 (1923); Hewlett, Phys. 


Rev. 17, 284 (1921). 
4 Kustner, Phys. Zeits. 46, 845 (1932). 
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TABLE I. Alomic absorption coefficients in U. 
























































Line din A ma 
Sr Kay 0.8734 2.90 x 10-20 
U Lp, 7531 1.94 
Cb Ka; .7446 1.94 
U Lp, 7185 1.84 
U_ LBs 7088 3.81 
U Ly .6136 2.53 

069 ——— 

+L g Discontinuity; 

050 iy 

: 

0.40 

O30r— —— = _ ——— 

O20 . 

76 182 186 190 toa 
Log 


Fic. 2. Graph showing logarithm of atomic absorption 
coefficient for uranium plotted against logarithm of wave- 
length (log ue as a function of log \). 


absorption limit is less than 0.7208A, since the 
line U L8,=0.7185A lies on the long wave-length 
side of the limit. 

The quantity 1.27/2.27=0.56 represents the 
fraction of the radiation on the short wave- 
length side of the Ly limit which is absorbed 
photoelectrically in the Ly level. Since the scat- 
tering term may be assumed negligible, p’= 7’ 
and hence 7,/=0.56y’. 


RELATIVE INTENSITIES OF THE MoO Ka TO THE 
Mo K& IN THE FLUORESCENCE SPECTRUM 


A metal x-ray tube with a rhodium target was 
operated at 40 kv and 10 milliamperes, the molyb- 
denum scatterer being placed immediately in 
front of the aluminum window of the tube. The 
distance from the center of the rhodium target to 
the molybdenum scatterer was approximately 
2.5 cm. A single crystal ionization spectrometer 
was used, the molybdenum fluorescence beam 
being collimated by two slits 0.8 mm wide, 2.4 
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cm high, and 25 cm apart. The ionization 
chamber was made of copper 20.3 cm long with 
mica windows, and filled with methyl iodide at 
20.8 cm pressure. The ratio of the peak intensities 
of the ay to B12. read from the graph is 4.85 to 1. 
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Fic. 3. Graph of intensity against crystal angle for K 
fluorescence spectrum from molybdenum. 


The true base line occurs at 0.048 mm per sec., 
found by moving the ionization chamber from its 
correct position for reflection from the crystal. 
This intensity ratio has to be corrected for: (1) 
the relative absorption in the air between the 
scatterer and the window of the ionization 
chamber, a factor of 1.014; (2) the relative ab- 
sorption in the mica window of the ionization 
chamber, a factor of 1.014; (3) the relative ab- 
sorption in the methyl iodide, a factor of 0.815; 
(4) the relative absorption in the Mo scatterer, a 
factor of 1.375; and (5) the reflection coefficients 
of the two wave-lengths from the calcite crystal, 
a factor of 0.965. This gives a final result of 5.42 
for the relative intensities of the Mo Kay, 
Mo KB .2in the fluorescence spectrum. In order to 
check this value with that given by Meyer for 
these lines in the characteristic spectrum, curves 
were drawn for the Ka;, Kae with peak intensities 
of 100 and 50.6 and displaced relative to one 
another by 2.5 minutes of arc, and similarly for 
the KB,, KB. The value obtained in this manner 
gave 5.9, which would appear to be within the 
limits of experimental error. Measurements are at 
the present time being carried on in this labora- 
tory by Dr. J. H. Williams on relative intensities 
of lines in the characteristic spectrum, and pre- 
liminary results indicate that Meyer’s value of 8; 
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is somewhat low. The calcite crystal which was 
used here was one previously calibrated by 
Allison (crystal VA) by means of a double 
crystal spectrometer. 


THE RELATIVE INTENSITIES OF THE Lay TO THE 
LB, LINES IN THE FLUORESCENCE SPECTRUM 
FROM THE Ly;; LEVEL OF URANIUM 


The apparatus used in this experiment was the 
same as in the previous one, except the rhodium 
target was replaced by a molybdenum target, and 
the molybdenum scatterer by a uranium scat- 
terer. The x-ray tube was operated at 40 kv and 
20 milliamperes, and the ionization chamber 
filled with methyl bromide at atmospheric 
pressure. The true base line occurs at 0.05 mm 
per sec., found by moving the ionization chamber 
from its correct position for reflection from the 
crystal. Lines from the LZ; and Ly levels appear 
in the graph but their relative intensities, with 
respect to the lines from the Ly1 level, cannot be 
compared with that from primary x-rays from 
a U target. The ratio of the peak intensities of 
the ay to 62 read from the graph is 2.22. Both the 
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Fic. 4. Graph of intensity against crystal angle for the 
fluorescence spectrum from the Lin level of uranium. 


% Allison, Phys. Rev. 41, 1 (1932). 
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a and §-lines are of shorter wave-length than the 
K limit of bromine, so that the bromine K radia- 
tion will be excited in the ionization chamber and 
part of this will escape to the walls of the 
chamber. The method of correction for this fol- 
lows that given by Allison and Andrew,’ and 
introduces a factor of 1.09. 

The other corrections follow as in the previous 
experiment: (1) correction for absorption in the 
air and mica window, a factor of 1.053; (2) cor- 
rection for relative absorption in the U, a factor 
of 1.62; (3) correction for relative reflection from 
the calcite crystal, a factor of 0.903. This gives a 
final result of 3.73 for the relative intensities of 
the U LayLfe, in the fluorescence spectrum. 

The relative intensities of these lines in the 
characteristic spectrum were determined by 
Allison’® and his results give the value of 3.78. 
Thus it would appear that the same relative in- 
tensities are obtained for the lines from the Zyy; 
level in the characteristic spectrum as in the 
fluorescence spectrum. 


CALCULATION OF THE FLUORESCENCE YIELD 
FROM THE URANIUM Ly; LEVEL 


The diagram of the apparatus used is shown in 
Fig. 1. The fluorescence molybdenum radiation 
passed through two circular holes 1 cm in diam- 
eter and 20 cm apart, and was visible on a 
fluorescent screen. The ionization chamber was 
10.7 cm in diameter, and 13.2 cm long, and the 
windows at the front and back were made so that 
none of the uranium radiation struck the walls or 
the back of the chamber. 

A difficulty was encountered which was due to 
the radioactivity of the uranium, or rather to the 
B-particles given out by the radioactive products 
UX,, UX». The ionization caused by the §-par- 
ticles was eliminated by interposing a piece of 
paraffin 1.485 cm thick between the uranium and 
the ionization chamber. This of course reduced 
the intensity of the fluorescence beam from the 
uranium; but not too seriously. 

The intensity of the direct beam from the 
molybdenum was reduced to measurable values 
by using absorbing screens of silver of known 
thickness. Silver was chosen because it can be ob- 
tained pure, and also no silver K radiation will be 


16 Allison, Phys. Rev. 32, 907 (1926). 
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excited. The absorption coefficient determined in 
this manner will be less than the correct absorp- 
tion coefficient for silver since some of the scat- 
tered radiation can enter the ionization chamber. 

In a similar manner the absorption coefficient 
for paraffin cannot be taken directly from tables, 
but was found approximately by finding the ab- 
sorption coefficient of paraffin for the direct 
molybdenum beam. The value found in this ex- 
periment for \=0.698A, which is weighted mean 
wave-length from Mo, is 0.384. From data given 
by Allen!’ the absorption coefficient for paraffin 
for \=0.698A is 0.523. It was therefore assumed 
that the absorption coefficient for \=0.876A, 
which is the weighted mean wave-length from the 
uranium, should be, for this experimental ar- 
rangement 0.384X0.901/0.523=0.661 where 
0.901 is the value given by Allen for \=0.876. 
There may be here a large source of error, for the 
scattering of the molybdenum and the uranium 
beams will not be the same, but this is some ex- 
tent taken up by the difference in the geometry 
of the two beams. The molybdenum beam is more 
or less a parallel cylindrical one, while the 
uranium beam is a divergent one. 

The intensities observed for the uranium 
fluorescence beam must then be multiplied by 
e°-982 or 2.67. 

Five separate determinations were made with 
different voltages and currents through the x-ray 
tube, and different pressures of methyl iodide in 
the ionization chamber. Sufficient intensity was 
found when the tube was operated at 30 kv and 
10 milliamperes. 

From Eq. (3) the contribution from the Mo 
Ka, and U La; is given by substituting 0.56y’ for 
rz’. Similarly, the contribution from the Mo Ka, 
and U Lag will be given by replacing 100/149.8 
by 11/149.8 in Eq. (3) and letting u”’ be the ab- 
sorption coefficient of U Lain U. Similarly for all 
the lines. The distance r from the center of the 
uranium to the center of the diaphragm S: was 
9.54 cm, and A” the area of S, was 20.27 cm’. 
The intensity I’ of the fluorescence radiation 
from the uranium will be given by 


I” =wz(SitS2+53+S;) A" /4rr* 
where S;, So, S3, S; represent the contributions 


17 Allen, Phys. Rev. 28, 907 (1926) 
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made by the Mo Ka;, Kae, KB;, KBs lines, re- 
spectively. Typical values of these contributions 
are S,=0.210, S.=0.1065, S;=0.047 and S, 
= (0.008. These values include the corrections for 
the absorption of the various wave-lengths in the 
methyl iodide. Hence 


I” =0.371w,I’ /4rr’. 


I" and I’ are proportional to the ionization cur- 
rents observed for the uranium fluorescence beam 
and the direct molybdenum beam, respectively, 
after each has been corrected for absorption. 
After correction for absorption in the paraffin, 
_ 


= 25.6 mm per sec. After correction for ab- 


sorption in the silver, J’= 5822 mm per sec. Thus 
25.6 1 


4rr? 
I, =— 


5822 0.371 A” 








= 0.67. 


The value given above was also the average of all 
the readings. 

Another method of calculation, which is con- 
siderably simpler than the one given above is 
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that of taking average values for the wave- 
lengths and absorption coefficients. The weighted 
mean average of the wave-length of the molyb- 
denum fluorescence beam is \’=0.698A, and for 
the uranium beam is \’’=0.876A. These values 
are taken from the fluorescence spectrum re- 
ported in this paper. The values of yw’ and yp”, 
the atomic absorption coefficients of the molyb- 
denum and uranium beams, respectively, in 
uranium are 3.65 X10-*, 2.92 x10-*. The calcu- 
lation then gives w,= 0.66. 

The fact that these two calculations give values 
so close to one another, may be taken to show 
that the intensities of the lines in the character- 
istic and fluorescence spectra, used here, are very 
nearly the same. It should be noted even if 
Meyer’s value of 8: is as much as 5 percent too 
low, the change in the result for the final value 
for wz would not be more than 0.5 percent. 

I take great pleasure in thanking Professor 
S. K. Allison for the suggestion of the problem 
and his useful advice throughout the course of 
the work. 
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Attempts to Observe an Electron Affinity Spectrum 
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The failure of previous attempts to observe an electron 
affinity spectrum (emitted by the combination of neutral, 
electronegative atoms and electrons) is explained by 
inadequate experimental methods. On the basis of recent 
results regarding combination spectra of positive ions and 
electrons, new attempts have been made to observe the 
electron affinity spectrum of atomic iodine by using three 


methods (hollow cathode, positive column with addition 
of a rare gas, glowing filament). No new spectrum was 
observed. This failure is discussed theoretically. It seems 
that the capture of electrons by halogen atoms is an 
improbable process as compared with the combination of 
positive ions and electrons. 





I. PROBLEM AND PREVIOUS WORK 


HE electron affinity of the halogen atoms, 
i.e., their ability to form negative ions, is 
known from various observations: the existence 
of polar molecules and crystals, the low mobility 
of the negative particles in their vapors, the nega- 
tive values of e/m in the “positive ray.’’ The value 
of the electron affinity (i.e., the energy difference 
between the atom and electron separated and the 
negative ion) has first been derived by Born! from 
the grating energy of the alkali halogen crystals. 
Another determination of the electron affinity of 
the iodine atom, based on the thermal equilib- 
rium at high temperature of CsI vapor, has been 
made by J. E. Mayer.? The iodine atom has an 
electron affinity of 74.2+1.5 k cal. 

These results establish the existence of the 
negative ions and its energy of formation. They do 
not give evidence, however, by what process this 
ion is formed. Franck*® suggested that a neutral 
gas atom as well as a positive ion, capturing a 
free electron, might radiate the energy of forma- 
tion plus the kinetic energy the electron had be- 
fore this process. The corresponding spectrum, 
therefore, would consist of a continuous band 
starting from the frequency given by the energy 
of electron affinity, spreading to short waves. The 
intensity distribution within this continuous 

1M. Born, Problems of Atomic Dynamics, Cambridge, 
1926, p. 170. 

2J. E. Mayer, Zeits. f. Physik 61, 798 (1930); J. E. 
Mayer and Lindsay Helmholz, Zeits. f. Physik 75, 29 
(1932). 

3 J. Franck, Zeits. f. Physik 5, 428 (1921); cf. K. Fajans, 
Verh. d. D. Phys. Ges. 21, 714 (1919). 


band would depend upon the kinetic energy dis- 
tribution of the electrons (modified by some 
probability factor for the capture also depending 
on the kinetic energy). The observation of this 
spectrum would be important in that it would 
give evidence of the elementary process, and in 
addition yield by its long wave-length limit a 
straightforward determination of the energy of 
electron affinity. 

All previous attempts to observe electron 
affinity spectra of gaseous atoms or molecules 
have failed. Two continuous bands observed 
in the positive column of the electric discharge 
through iodine vapor have been ascribed to 
the process under discussion. Both of them 
turned out later to belong to the band spec- 
trum of the iodine molecule.‘ These failures in 
observing the electron affinity spectrum in the 
positive column can be understood on the basis of 
more recent information. In this source of light 
even the combination of positive ions and elec- 
trons—in spite of the Coulomb force acting— 
does not play an important part, the reason being 
that only very slow electrons have an appreciable 
chance to be caught by ions.® By analogy we con- 
clude that in the positive column most electrons 
are too fast to be captured by neutral iodine atoms. 

Von Angerer and Mueller® investigated absorp- 


4 W. Steubing, Ann. d. Physik 64, 673 (1921); J. Franck, 


Zeits. f. Physik 5, 428 (1921); W. Gerlach and F. Gromann, 
Zeits. f. Physik 18, 239 (1923); O. Oldenberg, Zeits. f. 
Physik 25, 136 (1924). 

5 R. Seeliger, Phys. Zeits. 30, 346 (1929). 

6E. von Angerer and A. Miiller, Phys. Zeits. 26, 643 
(1925) and L, A. Miiller, Ann. d. Physik 82, 64 (1927). 
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tion spectra of alkali halogen compounds. In their 
vapors at high temperature a considerable con- 
centration of negative halogen ions is to be ex- 
pected. The neutralization of these negative ions 
(I->I+) requires the energy of electron affin- 
ity. The question is: does the neutralization take 
place by the action of light ? Continuous bands 
were observed, which the authors interpreted 
later, however, as belonging to the neutral alkali 
halogen molecules. Their failure to observe elec- 
tron affinity spectra might be due to the masking 
of a faint absorption by another stronger one. 
Franck and Scheibe’ have interpreted continu- 
ous absorption bands observed in the solutions of 
the alkali halogen compounds in water as electron 
affinity spectra of the negative halogen ion losing 
the excess electron under the influence of light. 
This hypothesis has recently been modified, 
however, by Franck and Haber.* They assume 
that the electron is not liberated but that it 
jumps from the halogen ion to a water molecule 
which in the same process splits up into H and 
OH-. The spreading of the continuous spectrum, 
therefore, is not related to the velocity distribu- 
tion of the electrons liberated but is the same 
phenomenon as the broadening of all spectral 
lines in liquids. This result does not solve the 
original, simple problem: the electron affinity 
spectrum of the undisturbed atom. This can be 
solved only by observing the gaseous state.°® 


7 J. Franck and G. Scheibe, Zeits. f. physik. Chemie (A) 
Haberband 22 (1928); G. Scheibe, Ber. d. Dtsch. Chem. 
Ges. 59, 1321 (1926). 

8 J. Franck and F. Haber, Ber. d. preuss. Ak. d. Wiss; 
Phys. math. Klasse XIII, 8 (1931). 

* In this argument it has been assumed that the electron 
affinity spectrum consists only of the continuous band 
without showing sharp lines, in: other words that the 
negative ion has no definite excited states. As the electron 
affinity is understood on the basis of the great stability of 
a completed group of electrons, there seems to be no 
reason to assume a stable negative ion with the excess 
electron in an excited state. It might be as stable as the 
negative ion of any other atom having no pronounced 
electron affinity. No corresponding atomic spectrum (line 
spectrum of a negative ion) has been observed. Bartlett 
(Nature 125, 459 (1930)) estimates that the Cl~ ion has 
no stable excited states at all below its neutralized state. 
Since he applies an extrapolation method, it might be 
safer to conclude that excited states, if there are any, 
should be very near the ionized state and therefore almost 
unstable. This result does not agree with Mecke’s tentative 
interpretation of certain intense band spectra of BeH and 
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So we come to the conclusion that as yet the 
problem has been attacked under conditions that 
from our present standpoint seem to be unfavor- 
able. Therefore, the failure to observe the spec- 
trum under discussion is not necessarily to be 
interpreted as due to the nonexistence of this 
spectrum. It will be worth while to search for it 
under more favorable conditions, 


II. EXPERIMENTS 


(a) Hollow cathode method 


Recent successes in producing combination 
spectra of positive ions and electrons” suggest 
various methods for producing the electron 
affinity spectrum of neutral halogen atoms. In the 
inside of a hollow cathode of a discharge tube 
containing helium Paschen" discovered marked 
continuous spectra connected with well de- 
veloped line series of the helium atom. As only 
slow electrons are able to recombine with ions it 
must be concluded that the concentration of slow 
electrons in the hollow cathode is considerable. 
This is largely due to the fact that no rapid ac- 
celeration exists, because, from symmetry, the 
electric field will be small along the axis of the 
hollow cathode. The corresponding experiment 
carried through with iodine vapor completely 
dissociated by heat should give a chance to ob- 
serve the electron affinity spectrum—slow elec- 
trons being present due to the properties of the 
hollow cathode, and iodine atoms due to the 
thermal dissociation of the vapor. 

The iodine atom is particularly suited for this 
experiment, because iodine is the only halogen 
vapor that can be dissociated to nearly 100 per- 
cent at easily accessible temperatures." Any new 
spectrum, therefore, that might come out under 
special conditions will be likely to belong to the 
iodine atom. The only remainder of the molecular 
spectrum that is known to appear in spite of 
complete dissociation is the band with the long 


MgH as due to negative molecular ions (Zeits. f. Physik 72, 
155 (1931)). 

10 A survey has been given by R. Seeliger, Phys. Zeits. 
30, 329 (1929). 

1 F, Paschen, Sitzungsber. d. Preuss. Akad. d. Wiss., 
p. 207 (1927); cf. R. A. Sawyer, Phys. Rev. 36, 44 (1930). 

2G, E. Gibson and W. Heitler, Zeits. f. Physik 49, 465 
(1928); cf. W. G. Brown, Phys. Rev. 38, 708 (1931). 











536 


wave-length limit 3460A. This is well outside the 
region between 3740 and 3910A in which, 
according to Mayers determination, we must 
search for the spectrum under discussion. It 
happens that this spectral range is appreciably 
free from intense bands even for a large concen- 
tration of molecules. 

The discharge tube was made of quartz with 
tungsten electrodes." Care has been taken that 
the inside of the hollow cathode was observed not 
against the glowing wall of the oven as a back- 
ground but, instead, against a quartz window. 
As it is known to be difficult to keep the iodine 
pure in the electric discharge, a continuous flow 
of iodine vapor was supplied. Iodine vapor reser- 
voirs with different pressures were attached one 
on each side of the discharge tube. The resulting 
flow of vapor washed all impurities away from 
the region of the electrodes and discharge. From 
time to time the whole system was evacuated 
after freezing the iodine. 

The results have been negative. Under no con- 
ditions (pressure varied up to 8 mm, current up 
to 90 m.a.) has a new continuous spectrum been 
observed. The plates showed mainly the line 
spectrum of the iodine atoms with faint traces of 
tungsten lines. In addition the I, band 3460 
came out with high intensity. It will be discussed 
in another note. Several plates showed traces of 
the known diffuse bands of the I, molecule. 

An argument against the hollow cathode 
method is that it might fail under the present 
conditions, because most of the free electrons 
might be captured by the atomic iodine imme- 
diately after leaving the cathode so that no free 
electrons would be present near the axis. In order 
to test this possibility, all parts of the hollow 
cathode discharge, particularly those parts in the 
neighborhood of the metal, have been carefully 
investigated by the spectrograph without any 
positive result. Only intensity differences within 
the line spectrum appeared." 

The presence of freeelectrons was checked by an- 
other method. Their disappearance by immediate 


18 For technical detail of the electric discharge in iodine 
cf. H. F. Fruth (Phys. Rev. 31, 614 (1928). Special stop- 
cock grease has been used as advised by M. Bodenstein. 

4 E. Miiller (Ann. d. Physik 14, 851 (1932)) proved that 
glowing platinum in iodine vapor emits electrons, not 
negative ions, 
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capture would produce a fundamental change in 
the discharge, as compared, for example, with 
the discharge through hydrogen or nitrogen. 
The conductivity is largely due to ionization by 
impacts of free electrons. This process cannot be 
simply replaced by impacts of ions only, particu- 
larly at high pressure. The reason for this is that 
only the light electrons keep their kinetic energy 
in the numerous elastic impacts not leading to 
excitation and ionization and thus can gradually 
gain energy from the electric field. The heavy 
ions, on the other hand, give up some kinetic 
energy in every elastic impact. In order to store 
up in a heavy ion the kinetic energy required for 
ionization in spite of collisions, enormous fields 
are necessary particularly at high pressure. More- 
over for ionizing an atom by impact, the heavy 
ion must have a much larger kinetic energy than 
the electron which requires just the limiting value 
of the ionization energy.'® Hence the potential 
difference between the electrodes in a discharge 
tube from which all free electrons are taken away 
will be very large compared with the potential 
difference, e.g., in the hydrogen or nitrogen dis- 
charge. The comparison of potential differences, 
therefore, for atomic iodine and hydrogen or 
nitrogen will give evidence whether or not an 
appreciable number of free electrons still exist in 
the iodine discharge. In the same discharge tube 
for the same pressure (5 mm) and temperature 
(1000°C), the current-potential-curve for iodine 
yielded even smaller values of potential difference 
than in hydrogen. We must conclude that there 
was some concentration of free electrons left. It is 
true that this observation was concerned not only 
with the potential drop in the positive column 
itself but that in the measurement the cathode 
drdép was included. Actually in the positive 
column the potential drop in iodine (although 
mainly molecular) has been found'® to be some- 
what higher than in nitrogen but of the same 
order (75 instead of 43 volts/cm). Hence it seems 
certain that in the discharge through atomic 
iodine the concentration of free electrons might 
be reduced by the formation of negative ions, 
but it cannot be very small as compared with the 
better known discharges through hydrogen or 

6 J. Franck, Zeits. f. Physik 25, 312 (1924); G. Joos and 


H. Kulenkampf, Phys. Zeits. 25, 258 (1924). 
16 Matthies, Ann. d. Physik 18, 493 (1905). 
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nitrogen.” The failure to observe the electron 
affinity spectrum cannot be due to the complete 
absence of free electrons. 


(b) Positive column in a mixture with helium 


It still might seem questionable whether the 
negative results of the experiments just discussed 
is due to the nonexistence of the electron affinity 
spectrum or to unknown properties of the hollow 
cathode discharge in the halogens. It is desirable, 
therefore, to search for the same spectrum under 
other widely different conditions that are also 
known to be favorable for recombination spectra 
of ions and electrons. Intense spectra of this kind 
have been discovered recently by Krefft!* in the 
positive column of the electric discharge through 
metallic vapors with the addition of rare gases. 
For a group of metals Krefft was able to compare 
the probability with which they emit the re- 
combination spectra. Under most conditions the 
radiation emitted from the positive column of the 
ordinary discharge is known to be due to previous 
excitation by electron impact but not to recombi- 
nation of ions and electrons,” the reason being 
that the concentration of slow electrons, which 
alone have a good chance to recombine, is too 
small. Under the conditions of Krefft’s experi- 
ment, however, it is not easy to understand the 
large concentration of slow electrons indicated 
by the intense combination spectra of ions and 
electrons (continuous bands connected with series 
limits). It seems possible that with He added the 
conditions are particularly favorable for the 
slowing down of electrons, because, as Sommer- 
meyer”? found out, the energy input in the electric 
discharge through He is largely used up for the 
numerous purely elastic collisions between the 
electrons and He atoms.*! In the present experi- 


7 This conclusion is supported by the small probability 
of the formation of negative ions in (molecular) iodine 
derived from electrical measurements by W. Hey and A. 
Leipunski (Zeits. f. Physik 66, 669 (1930)); cf. F. L. 
Mohler, Phys. Rev. 26, 614 (1925) and L. B. Loeb, J. 
Frank. Inst. 197, 45 (1924). 

*H. Krefft, Phys. Zeits. 32, 948 (1931); Zeits. f. Physik 
77, 752 (1932). 

9 Cf. J. Franck and P. Jordan, Anregung von Quanten- 
Spruengen durch Stoesse, Berlin, 98 (1926). 

*°Sommermeyer, Ann. d. Physik 13, 315 (1932); ef. H. 
Kessel, Zeits. f. Physik 70, 614 (1931). 

*t Results of Mohler (Bull. Am. Phys. Soc. 7, No. 7, 
p. 8 (1932)) just published suggest another reason for the 


ments, therefore, He of 7 mm pressure has been 
added +o atomic iodine (iodine vapor heated to 
1050°C). The positive column was projected 
“end on’”’ on the slit of the spectrograph. The 
pressure of the iodine was varied between 0.2 
and 1.1 mm. Not a trace of a new continuous 
spectrum was observed although the exposure 
was long enough to bring out some traces of the 
bands of molecular iodine. This negative result 
points again to a low intrinsic probability of the 
combination of the neutral iodine atoms with free 
electrons taking place with radiation. 


(c) Glowing filament method 


In order further to test the negative result of 
the experiments discussed it seemed desirable to 
try still another experiment. A straightforward 
method for producing electrons with very low 
velocity at a considerable concentration has been 
very successfully applied by Mohler.” He ob- 
served recombination spectra in the close neigh- 
borhood of a glowing filament, where the velocity 
distribution of the electrons is mainly determined 
by the temperature of the metal from which they 
have just been evaporated without being ac- 
celerated as yet by the electric field. In a quartz 
tube with iodine vapor heated to 1000°C a glow- 
ing tungsten filament was used; tungsten inter- 
feres least with the atomic iodine and, therefore, 
gives the purest conditions. 

Again no new continuous spectrum came out. 
This observation, however, does not give as good 
evidence as the hollow cathode method because 
the exposure could not be extended to a consider- 
able length of time because of the superimposed 
spectrum of the glowing filament. It is not pos- 
sible to screen it off so completely that in a very 


slowing down of electrons in the positive column, Mohler 
observed combination spectra of positive ions and electrons 
in the positive column of a discharge through pure caesium 
at a considerable pressure. It might be that the exception- 
ally low value of the excitation potential of caesium (1.38 
volts) together with the large probability of the lowest 
quantum transition in the alkali atoms is responsible for a 
large concentration of slow electrons. Also from this point 
of view it seems worth while to investigate the positive 
column in atomic iodine at various pressures because the 
iodine atom has an excitation potential as low as 0.94 volts 
(L. A. Turner, Phys. Rev. 27, 406 (1926)) although with 
a much smaller probability of transition. 
*2 F. L. Mohler, Phys. Rev. Sup. 1, 221 (1929). 
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long exposure it gives no continuous background 
even in the near ultraviolet. 

From the three groups of experiments dis- 
cussed it will be concluded that the neutral 
iodine atom has a smaller chance than a positive 
ion to capture a free electron and radiate the 
energy liberated although the energy value of 
this process is about the same as the ionization 
energy of caesium whose ions show an intense 
recombination spectrum. It must be admitted, 
however, that this negative result cannot be 
proved by experiments with the same certainty 
as a positive observation. The mechanism of the 
discharge is not known so completely that we can 
with full certainty predict the most favorable 
condition for a certain elementary process. In 
particular it might be that still slower electrons 
than those appearing in the discharges just dis- 
cussed show a larger tendency of combining with 
neutral halogen atoms.” 


III. THEORETICAL DISCUSSION 


In order to explain the failure in all attempts to 
observe the electron affinity spectrum two rea- 
sons will be discussed, one dealing with the veloc- 
ity of the electrons in the discharge, the other 
with the process of combination itself. 

A great difference between the two cases of 
combination—positive ions and neutral halogen 
atoms capturing electrons—will be due to the 
long range over which the Coulomb force of the 
ion affects free electrons, whereas the halogen, 
being a neutral atom, is surrounded by a field 
which, although unknown in detail, certainly 
becomes negligible within a very small distance. 
Therefore classically the electron to be captured 
must approach the halogen atom much more 
closely than the positive ion with the result that 
the group of electrons which can be captured by 
the halogen atoms will be comparatively small. 

This comparison, however, between combina- 
tion processes, the energy of which is radiated, 
must take into account the selection rules. Ac- 
cording to these rules the angular momentum of 
the electron can change in the radiation process 


*3 Recombination spectra of positive iodine ions and free 
electrons were not observed either. This is not surprising 
as Krefft, too, failed to observe this type of spectrum for 
certain groups of elements. 
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only by the quantized amount /:/27. Let us con- 
sider electrons with the average value of the 
kinetic energy. The capture of an electron by a 
positive ion will frequently happen over a com- 
paratively large distance because of the long 
range of the Coulomb field. In this case the elec- 
tron will have a large angular momentum with 
regard to the ion. After emitting the quantum 
h/2zx it will still retain an appreciable value and, 
therefore, go to a quantum orbit of the neutral- 
ized atom with large azimuthal quantum number, 
A frequency will be emitted belonging to a high 
subordinate series. The continuum connected with 
a principal series, however, will be emitted by the 
smaller group of electrons passing the ion with 
smaller angular momentum, i.e., at a smaller 
distance. This smaller number of processes will 
give rise to a fainter spectrum. This argument 
holds only if we deal with electronic tempera- 
tures such that the electrons have angular mo- 
menta of the order of magnitude #/2z7 when 
passing the ions at distances of the order of mag- 
nitude of atomic dimensions. Actually for an 
electron temperature of 5000°C, an electron 
which has the most probable velocity and an 
angular momentum with respect to the ion of 
h/2x, will pass the ion at a distance of 210-8 
cm; this is well within the range of atomic dimen- 
sions. On this basis we understand the fact, re- 
ported by Krefft, that in the combination spectra 
of positive ions and electrons the continuum con- 
nected with the principal series is much less in- 
tense than the one connected with the first sub- 
ordinate series. 

In the electron affinity spectrum of the halogen 
atoms, as discussed in Section I, we expect to ob- 
serve only the continuum representing the 
principal series. It is not surprising that because 
of the more constricted field of force of the 
halogen atoms, the probability of this process 
comes out still smaller than for the positive ions, 
although the energy with which the halogen keeps 
the electron is of the same order. 

This interpretation, based on the velocity dis- 
tribution of electrons in the discharge, does not 
preclude the inverse process, the separation of the 
excess electrons from the negative ions by ab- 
sorption of light, provided the ions are produced 
in sufficient concentration. For the initial state of 
this process, the negative ion, is independent of 
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the velocity distribution of electrons. Therefore, 
it does not contradict the theory that in absorp- 
tion in solutions the electron affinity spectrum has 
actually been observed by Scheibe, although 
modified by the water to such an extent that no 
conclusion can be drawn regarding the same ab- 
sorption process (its wave-length and its proba- 
bility) in the vapor. The concentration brings in 
another difference between vapors and solutions 
in water. For example, in a 1/100 normal solution 
of KI the I~ ions have a partial pressure of 1/4 
atmosphere; that is a very large pressure as com- 
pared with the value in electric discharges. 

The theoretical result discussed, regarding the 
halogen atom in the gaseous state, might be modi- 
fied in one respect by a recent result derived from 
wave mechanics by Jen.** As yet it has been as- 
sumed from analogy with the behavior of positive 
ions that the electron affinity spectrum of certain 
neutral atoms it if exists at all will have a large 
intensity near its long wave-length limit. A differ- 
ent result, however, has been derived by Jen. In 
any case (positive ion or neutral atom) the effec- 
tive cross section for the capture of electrons ap- 
proaches infinity as a limit with decreasing kin- 
etic energy of the electrons; the concentration of 
slow electrons, however, approaches zero as a 
limit. The intensity observed (product of probabil- 
itv factor and concentration factor) comes out 


**C. K. Jen, Phys. Rev. 42, 588 (1932). Cf. the following 
paper. His results are partly supported by the probability 
for the formation of negative ions as derived from electrical 
measurements by W. Hey and A. Leipunski (Zeits. f. 
Physik 66, 669 (1930)). Their results, however, do not 
apply directly because their experiments were performed 
with molecular iodine vapor. 


with a finite limiting value. For the electron 
affinity spectrum of hydrogen—contrary to the 
spectrum of combination of positive ions—the 
theoretical intensity curve shows a pronounced 
maximum at a considerable distance from the 
limit. Hence it might be difficult to observe the 
theoretical long wave-length limit (defined by the 
energy of electron affinity) and to identify the 
spectrum, even if the maximum were intense 
enough to be photographed. 

Jen’s theory certainly does not lead to the re- 
sult that the capture of an electron by a neutral 
atom resulting in radiation is a “forbidden” 
process. It seems hardly possible, however, to 
estimate from wave mechanics the actual inten- 
sity of the electron affinity spectrum of hydrogen; 
and we certainly must not correlate Jen’s theoret- 
ical results quantitatively with observations in 
the halogen vapors. 

The main result is that the chance for observ- 
ing electron affinity spectra in certain parts of the 
electric discharge seems to be a great deal smaller 
than the chance for observing combination 
spectra of positive ions and electrons.*® 


*% It is interesting to compare the negative result of the 
present spectroscopic investigation with the results of 
Hogness and Harkness (T. R. Hogness and R. W. Hark- 
ness, Phys. Rev. 32, 784 (1928)), who investigated the 
ionization processes of iodine by the mass-spectrograph. 
Working with molecular iodine they had no chance to 
investigate the atomic process under discussion (I +«—I~). 
They found out, however, that negative atomic ions I~ are 
formed by neutral molecules capturing electrons and 
breaking up at the same time into a neutral atom and a 
negative ion (I2+e«~I+I-). Not even I;~ ions are formed 
simply by neutral molecules I, capturing free electrons but 
by a more complicated process 
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The Continuous Electron Affinity Spectrum of Hydrogen 
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A wave-mechanical calculation has been made of the 
intensity distribution in the continuous emission spectrum 
due to the capture of electrons by normal hydrogen atoms 
and the coefficient of absorption in the corresponding 
absorption spectrum due to the neutralization of negative 
hydrogen ions. Approximate wave functions for the con- 
tinuous range of upper energy levels have been combined 
with the discrete state wave function for H~ given by 
Hylleraas to give the matrix components for the desired 
transitions. The computed results are graphically illus- 


trated and their significance discussed. It is shown that 
an experimental attempt to find the electron affinity 
emission spectrum of normal hydrogen atoms may en- 
counter at least two difficulties: (1) insufficient intensity 
compared to that in experiments on the recombination 
between H?* ions and electrons, and (2) spectral distribution 
such that the spectral limit is too diffuse for positive 
identification. The relation of these results to the interpre- 
tation of experiments on the electron affinity spectrum of 
halogen atoms is discussed. 





INTRODUCTION 


OLLOWING the first proposal by Franck,! 
many attempts have been made to determine 
the electron affinity of electronegative elements 
by spectroscopic means. Certain unidentified 
continuous spectra have often been taken as 
those related with electron affinity, but such 
identifications have proved to be either false or 
uncertain.? The recent experiments of Oldenberg 
on the emission spectrum of iodine atom also 
point to a negative result.* A theoretical explana- 
tion for this experimental failure is obviously 
demanded. However, a complete answer to the 
problem is difficult because of the complexity of 
the atoms involved. As a start we propose to 
treat the electron affinity spectrum of the hydro- 
gen atom, which is the simplest of such problems, 
and hope to furnish probable explanations for the 
more complicated halogen atoms. 

If free electrons are present in the midst of 
normal hydrogen atoms, there may be a proba- 
bility of combination between them to form 
negative hydrogen ions. The change of energy in 
such a transition is released as radiation, which, 


* O. Oldenberg, Phys. Rev. 43, 534 (1933). 
1 J. Franck, Zeits. f. Physik 5, 428 (1921). 
2W Steubing, Ann. d. Physik 64, 673 (1921). 
O. Oldenberg, Zeits. f. Physik 25, 136 (1924). 
E. V. Angerer und Miiller, Phys. Zeits. 26, 643 (1925). 


if spectroscopically observed, is expected to form 
a continuous emission spectrum extending toward 
the short wave-length side from an edge whose 
frequency is determined by the electron affinity 
of the atom. Conversely, a continuous absorption 
spectrum will be formed when negative hydrogen 
ions are broken up into neutral atoms and elec- 
trons under the influence of light. 

The electron affinity of a hydrogen atom has 
been calculated by Hylleraas* and Bethe* froma 
general treatment of the two-electron problem. 
In spite of the discordant theoretical results' 
which had appeared previously and the absence 
of sufficient experimental data, their calculations 
seem to be quite definite and convincing, 
especially in view of the remarkable success of 
Hylleraas’ method on other examples of the two- 
electron problem.® With the help of their char- 
acteristic wave function of a negative hydrogen 
ion and our present approximate solution of the 
wave equation in the continuous state, we shall 
calculate the matrix component of transition for 
the continuous spectrum, from which the emis- 
sion intensity and absorption coefficient can be 
readily deduced. 


8 E. A. Hylleraas, Zeits. f. Physik 60, 624 (1930). 
4H. Bethe, Zeits. f. Physik 57, 815 (1929). 
5 L. Pauling, Phys. Rev. 29, 285 (1927). 
H. Briick, Zeits. f. Physik 51, 707 (1928). 
6 E. A. Hylleraas, Zeits. f. Physik 65, 209 (1930). 


540 


it 


“< 


orn Ss 


of 


nn en ee! 





THEORY OF ELECTRON AFFINITY SPECTRUM 541 


I. APPROXIMATE SOLUTION OF THE WAVE EQUATION IN THE CONTINUOUS STATE 


The Schrédinger differential equation for two electrons acting in the field of a single fixed proton 
charge is 
Viv tveet [Ert+2(1/n+1/re—1/riz) pW =9, (1) 


where Ez, the total energy, is measured in terms of Rh and the 7’s are in terms of a= h?/4r?me?. 
Since the above equation cannot be solved exactly by any ordinary process, we shall resort to a 
method of approximation. It may be assumed that the electron in the normal state of hydrogen is 
not appreciably disturbed by the external electron so that the interaction term is considered to affect 
the potential function of the latter only. Separation of variables is thus made possible. 
Let 
¥ = Fo(a)f), Er=Eot+«, (2) 


where F,(i) and Ey denote respectively the characteristic wave function and energy value of the 
electron belonging to the hydrogen atom in the normal state; f(j) and ¢« denote respectively the con- 
tinuous wave function and energy of the external electron. The letters 7 and j are used to designate 
the “‘inner’”’ and “‘outer”’ electrons; either of them may be numbered as 1 and the other as 2. 

Eq. (1) is then separated into 


V2 Fo(t) +(Eo+2/r;) Fo(i) =0, (3) 
Vif(A) + [e4+2U(r;) Lf(Z) =0. (4) 

The normalized characteristic solution of Eq. (3) for the lowest energy level is 
Fy(i) = abe", (5) 


The solution for Eq. (4) may be sought first by the usual method of separation of variables in 
spherical coordinates. The angular component of f(j) has the usual expression P,” (cos @;)e'"*i while 
the radial component, R,,; (r;), must satisfy the differential equation 


d?(r;R)/dr7?-+ | e+2U(r;) —l(l+1)/r7 }(r7;R) =90, (6) 
where / is an integer. 
Now we must find an expression for U(r;), which consists of the potential energy due to nuclear 
attraction and the interaction energy. To a close approximation, the interaction term is the average 
of 1/r;; over all the configuration space of the inner electron. The result is 


U(r) =-— | ——. (7) 


’; fos) Vij 


1 F(i)dr; 


The integral is readily carried out by expanding 1/r;; in spherical harmonics with the result 
U(r;) =(14+1/rje". (8) 


Even with the above expression, a direct solution for Eq. (6) is still quite difficult. However, Eq. (8) 
can be very roughly represented by 
U(r;) ~Z*/1r; (9) 


where Z* is a constant (less than one in this case) and is to be interpreted as the effective nuclear 
number with respect to the external electron. Z* can be adjusted to make the expression in Eq. (9) 
fit a limited portion of the curve represented by Eq. (8). As it will be shown later, / in Eq. (6) must 
be 1 by the selection rule for the transition under consideration, so that the term /(/+1)/r;? pre- 
dominates for r;<1. This means that we will only have to secure reasonable fitting for large values 
of r;. The numerical determination for Z* will be discussed later. 
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By our above assumption Eq. (6) becomes 
d?(r;R)/dr?+ {e+2Z*/r;—l(l+1)/r7?} (7;R) =0. (10) 


It is well known that the above differential equation is satisfied by the Laplace contour integra]? 
which can in turn be transformed into 


ene 

xe" pt —} - 

—2iR., i(7r;) =n efi(E+i_h) IZ € “(€—te}) i € tat. (11) 
—i(e)* 


The solution for f(j) is therefore simply the continuous wave function for the hydrogen atom with 
an effective atomic number Z* instead of 1. The normalization of such a function can be carried by 
the method of Fues® as follows: 





. ey + Ae, 
farte@ [ fe(jde=1 or 0 (12) 
according as ¢,<@«<e,+Ae or not. The result is identical to that of Fues: 
— PQl+1) @-m) 1} e8x/20(2" 74) . 
fe, 1, m(J) -| | ar : —— R,, i(r;)Pi™(cos 0;)e*%i, — (13) 
2 (l+m)!5 w[l2(6)) )44|TU+14+72(2*/ )}) | 


Fo()f., 1, m(j) is a particular solution for the wave function in the continuous state. Since 7, 7 can be 
either 1, 2 or 2, 1, there are two linearly independent solutions, the proper normalized linear combina- 
tions are, of course, the symmetric and antisymmetric ones,’ 


We, l, m(H°+e) == 2-1 { Fo(1)fe. l, m(2) + Fo(2)fe, i, m(1) a 14) 


The bracket (H°+e) is now added to denote the wave function of a normal hydrogen atom and an 
external electron before capture, in order to distinguish it from the discrete wave function. In the 
present special case, however, the characteristic wave function for H~ (see next section) is symmetric; 
the antisymmetric function represented by the lower sign in Eq. (14) can be omitted because it contrib- 
utes nothing to the matrix component. 


II. CHARACTERISTIC WAVE FUNCTION OF H~ OBTAINED BY HYLLERAAS 


By a variational method Hylleraas® obtained for the lowest energy level of H~ an approximate 
characteristic function in the following form, 


¥(H~) = Ce~eat) {1 +Briety(ri—re)*}, (15) 


where a, 8, y are constants and C is the normalization constant which can readily be shown to be 





; , 358 3y 77 By 9777} 
C=a'/r — 2 ss . 


Sa a 8a’ at 


The presence of the factor riz in Eq. (15) is somewhat troublesome when we come to the integration. 
But we can eliminate the difficulty by the following expansion 





co ag 1 
rne=? > ( - )"P.(c0s®), (16) 
v=0 \2v+3 2v—1 


where 








7L. Schlesinger,” Einfiihrung in die Theorie der gewéhn- | 8 E. Fues, Ann. d. Physik 81, 281 (1926). 
lichen Differential-gleichungen auf funktionentheoretischer | 9 E.g., A. Sommerfeld, Wave Mechanics, p. 229. 
Grundlage, § 67. 
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» (v—|pl)! 
P,(cosO)= 3} ——— P,4(cos 01) P,*(cos 02) ¢*#(%1-%) , (17) 
w=» (y+ |u|)! 


cos 0 =cos 6; cos 62+sin 4; sin 62 cos (¢:—¢2), 


¢= larger of the two radii, 7; and 72; x=ratio of the smaller to the larger radius. 


III. THe SQUARE OF THE MATRIX COMPONENT OF TRANSITION 


The matrix component of transition formed by y, 1, m(H®+e) and ¥(H~) in the directions x, y, z 
of the Cartesian coordinates may be defined as 


po l, m tee 
Fa tar™ [ve t, m(H®+e)5 yitye bW(H~)dxidyidz,dxed yodze. (18) 
Ze. 1, m = LS + Se 





Let A/?, 1, m be the square of the amplitude of the total matrix component per unit energy range. 
Then the sum over all allowable values of /, m is 
, ® M?, l, m= { [X., l, m|?+ | Y,, l, m|?+ \Z«. l, mi*}- (19) 
l,m l,m 
The three components can be evaulated when we transform Eq. (18) in terms of spherical coordinates 
and make use of the expressions derived for y¥., 1, m(H®+e) and ¥y(H~-). We immediately find that 
M?, 1, m, vanishes except for ]=1 and m=0 or +1. Therefore the external electron must have an 
angular momentum quantum number 1 in order to have a chance to be captured by the atom. 
Summing up the contributions due to the three different values of m, we arrive at the result 





Rr C e3rZ*/et 
a ee —|W)?, 1l=1;m=—1,0, +1 (20) 
ie 9 & |P[2+i(Z*/e)]|? 


where 


w= [fev a)rytare} [370 {1-+-y(ry—re)?} +67 {ri(x3/5 —x) +re(x2+3) }] 


—i(e)? 


rite! a e 2 
| ebre(F+-a(€) #4 le? (¢ —G(¢)3)!-i2 “ag tartare (21) 


The integration of Eq. (21) is effected first by inverting the order of integration and then applying 
the theorem of residues for the evaluation of complex integrals. The process is essentially similar 
to that used by Fues* and Sugiura," except that here we have to consider also the residue at infinity. 
Omitting all detailed steps, we write the final result as follows: 


. 24407 Z*(Z*+.6) | B(a+Z*) (- 4y ) (Z* —2a) 
Pn  netaneneiis —______—_—_—— ( -+——_ ‘on iclitesiaiihainaians 
(l1ta)® 1-74 /« (1+a)*Z*(Z* + €) 3 (1+a)*7 (a’?+e)3 





1307 —58(1+ a) } 2 : 307 2 2 : 
+ -—_-__—__—__-_- — [ett 22%a: a-z) }+ —— | —- a®—Z*a*+- (Z* — ela 
3(1+a)(a*+ €)? 5 (a?+ €)>L3 5 


zZ* ‘ | 2Z* € | a+Z* 2(2a+3-—Z*) 
+— (7e—22Z") | exp| ~~ — tan <|-9 ee ——__——— 
45 








l(1+a)!Z*(Z*4+ 6) (1+a)(Z*+.€) 
2 2(4a+2—2Z*) 22° é! || 
» —_— i. ( 


e a. 


+—_—_—— — ex 
(1ta)?{(2a+1)?+e}2 3(1+a) {(2a+1)?+}3) 








tan-! ——— (22) 
é 2a+1 


” Y, Sugiura, Scient. Pap. of Inst. of Phys. and Chem. Res. (Tokyo), No. 193 (1929). 
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IV. INTENSITY OF EMISSION, EFFECTIVE CROsS 
SECTION OF COMBINATION, AND THE ATOMIC 
COEFFICIENT OF ABSORPTION FOR H7- 


(a) To calculate the total intensity of emission 
one must assume a definite distribution of elec- 
trons while they are mixed with the atoms. Since 
we are here dealing with neutral atoms only, we 
shall assume that their effect on the distribution 
of electrons can be neglected for the present pur- 
pose. This cloud of free electrons may be con- 
sidered to behave exactly like an ideal gas in a 
fixed enclosure, having for the aggregate an elec- 
tron temperature, 7.. Hence, we can apply the 
results of the quantum statistics for free elec- 
trons. 

From the preceding section it can be shown" 
that (a?/Rh)>°:. nM? 1, m@E represents’ the 
amplitude squared of the matrix component for 
the transition of an electron having its energy 
lying between E and E+dE. Here a*/Rh is the 
conversion factor from the units used in the pre- 
ceding sections into c.g.s. units; E=Rhe. 
Therefore (e?a?/Rh)>°:, mM. 1.mdE is to be 
interpreted as the square of the amplitude of 
electric moment. If gz is the number of free elec- 
tron states in the energy interval between E and 
E+dE, then (e2a?/geRh)d 1, mM%, 1, ndE is the 
weighted mean of the amplitude squared of the 
electric moment. We need only to multiply this 
by the total number of electrons in dE to obtain 
the total amplitude squared of the electric mo- 
ment. By the well-known quantum statistics," 
the total number of electrons in the energy inter- 
val between E and E+dE is equal to gz-§ where 
® is the average number of electrons per state. 
Therefore the total amplitude squared of the 
electric moment is 


(ea Rh) D1, mM2, 2, mE. (23) 


Now the electric moment vibrates with a fre- 
quency v determined by the equation 


hv=E+e, (24) 


where $= electron affinity of the normal hydro- 
gen atom. By referring to the classical electro- 
dynamics for the radiation from a dipole oscilla- 
tor, we can get the following wave-mechanical 


1 E.g., G. Birtwistle, The New Quantum Mechanics, 
Chap. XXIX. E. Fermi, Zeits. f. Physik 36, 902 (1926). 
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equivalent for the intensity of the continuous 
spectrum: 


647+ vte*a? 


I(v)dv= N°. RS >, = Af? l, nd E, (25) 
3c°Rh 


where idv=dE, by Eq. (24); I(v)dv=time rate 
of emission energy in ergs per sec. per cc in the 
frequency interval between v and v+dvp; N° 
= number of normal hydrogen atoms per cc. The 
y-function may be assumed to be that of an ideal 
electron gas obeying the new statistics of Fermi 
and Dirac, namely 


F=1/ce*Te—1, (26) 


where 7.=electron temperature, and C= func- 
tion of 7. and electron density. Fermi!* showed 
that for very large T.,* C is so large that Ce®/*? 
>1. Hence we have the Maxwell’s distribution 


and the ¥-function is accordingly 
iy = LNeh?/(2emkT,)? Je-*#'*, (26a) 


where N* is the number of electrons per cc. 
Combining Eqs. (25) and (26a), we have 


I(v)  642tvte*a? 


NoNe  3¢3R 








h® . ye 
~e-FElkTe F M2 pm. (27) 


(2rmkT,)' adie 


(b) The effective cross section of combination 
for an atom (gq) may be defined as the cross sec- 
tion of a small sphere such that the number of 
electrons with velocity v striking the surface of 
the sphere will be equal to the number of elec- 
trons with the same velocity actually captured 
by the atom. If n(£)dE is the number of electrons 
per cc in the energy interval between E and 
E+dE (according to Maxwell’s distribution), 
then 


n(E)dE=[2rN*E}/(akT.)' Je~*'* dE, (28a) 
I(v)dv = N°*hv- qun(E)dE. (28b) 


Solving for g, we get 
qg=[16rhv?/3(mv)?c? Jea?d 1, m7, i,m. (29) 


 E. Fermi, Zeits. f. Physik 36, 902 (1926). 
* We are here interested only in large T.. (See Section V.) 
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(c) The atomic coefficient of absorption (k) 
for H~ may be similarly defined as the effective 
cross section of H~ for the capture of quanta of 
radiation. The number of quanta captured, or 
the number of H~ atoms dissociated, per second 
per cc in the frequency interval between vy and 
y+dy is clearly 


kcpdv/hv 


where pdv= radiation energy density in the same 
frequency interval. 

Assuming thermal equilibrium between emis- 
sion and absorption processes (extension of Ein- 
stein’s method to transitions between continuous 
and discrete states), Kramers’ and Milne’ 
derived the following relation between g and k 


g/k =2gh?v*/mvCc*, (30) 


where g is the statistical weight of the lower 
energy state in this case, g(H~)=1. The explicit 
expression for k in terms of the square of the 
matrix component of transition is 


k= (845v/3cRh)ea2D 1, mM2,1,.m- (31) 


V. NUMERICAL COMPUTATIONS 


The three constants in the H~ wave function 
given by Hylleraas are in the units here used: 
a=0.770, 8=0.308, and y=0.119. To fit the ex- 
pression in Eq. (9) with that in Eq. (8) for large 
values of the radial variable, Z* may take any 
value ranging between 0 and 0.5, depending upon 
which portion of the curve is to be fitted best. 
Various values of Z* in this range have been 
tried in the present calculation. Except for 
slight differences in numerical magnitudes, they 
give approximately the same distribution of 
Dd), m-lZ?, 1, m as a function of the energy param- 
eter «. For illustration, we have chosen for Z* a 
value of 0.42, although much smaller values 
will equally serve the purpose. The values for 
Y), » 2. 7, m are not critical to the variation of Z*, 
so that the general interpretation of results 
about to be obtained will not necessarily depend 
upon any particular value of Z*. 

The electron affinity of the normal hydrogen 
atom is given by Hylleraas as 0.715 v, with 


13H, A. Kramers, Phil. Mag. 46, 836 (1923). 
ME. A, Milne, Phil. Mag. 47, 209 (1924). 
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which we can calculate the frequency of radiation 
for any kinetic energy of the electron (by 
Eq. (24)). 

Now, we will consider the order of magnitude 
of the electron temperature (7.). Langmuir and 











a 02 03 0. 0506 €(Rh) 
17254 3607 2014 1397 A(A) 


Fic. 1. Intensity distribution at different 
electron temperatures. 


Mott-Smith” find a temperature of about 
30,000°C at 210-+ mm pressure in discharges 
through mercury vapor. Kopfermann and Laden- 
burg’ find convincing evidence for statistical 
equilibrium in discharges through neon at 1 mm 
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Fic. 2. Effective cross section of combination (qg) and 


coefficient of absorption for H™~ (&). 


pressure with an electron temperature of about 
22,000°C. From these evidences and others, it 
seems certain that the electron temperature 
under electrical discharge conditions must be 


4 J. Langmuir and H. Mott-Smith, Gen. Elec. Rev. 27, 
449 (1924). 

‘®H. Kopfermann and R. Ladenburg, Naturwiss. 5, 
513 (1931). 
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quite high. In the present calculations of the 
emission spectrum we shall choose several dif- 
ferent values of electron temperature to illustrate 
its effect upon the general distribution of emis- 
sion intensities. 

The accompanying graphs (Figs. 1 and 2) 
show the results of computation for the emission 
intensities at different electron temperatures, the 
effective cross section of combination, and the 
coefficient of absorption for H-. 


CONCLUSION 


(a) Discussion of results 

The computed results show that the intensity 
per unit frequency interval in the emission spec- 
trum has a low but finite value at the long wave 
limit. There is a point toward the short wave side 
of the limit where the intensity has its maximum 
value, which is displaced more and more away 
from the limit as the electron temperature is 
raised. Thus, at a very high electron temperature 
(e.g., 20,000°K) the intensity at the limit may be 
more than ten times smaller than the maximum 
value whose spectral position is considerably 
away from the limit (see Fig. 1). If the abso- 
lute magnitude of the intensity is low, then the 
limit would be so faint that it is beyond recogni- 
tion. The identification of the continuous spec- 
trum and the subsequent determination of elec- 
tron affinity by the position of the limit would 
become very difficult. 

The curve for the effective cross section of com- 
bination (see Fig. 2) is similar in appearance to 
that of the recombination between H* and elec- 
trons obtained by Stueckelberg and Morse,” ex- 
cept that the order of magnitude of the former is 
much smaller and it does not diminish as rapidly 
for high-velocity electrons. It must be mentioned 
here that the present calculations are based on 
random distribution of electrons while Stueckel- 
berg and Morse assumed a unidirectional stream 
of electrons. For a more exact comparison, we 
have calculated the effective cross section in 
exact accordance with Stueckelberg and Morse’s 
method. The resulting value is not shown in the 
figure, but is about 100-500 times smaller than 
for their case. If the concentration of the corre- 


17E, C. G. Stueckelberg and P. M. Morse, Phys. Rev. 
36, 16 (1930). 
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sponding reaction quantities (i.e., the number of 
H~- atoms and electrons in one case and the num- 
ber of Ht atoms and electrons in the other) are 
the same, the emission intensities would bear the 
same ratio. Hence the intensity of the electron 
affinity spectrum for normal hydrogen atoms js 
expected to be generally very low. 

The curve for the atomic coefficient of absorp- 
tion (see Fig. 2) has a maximum near the limit of 
the continuous absorption spectrum. The magni- 
tude of the coefficient is generally higher than the 
values for caesium given by Mohler'® (possibly by 
as much as 100-fold). If there is sufficient concen- 
tration of H~ ions present, it should be possible 
(apart from the experimental difficulties) to ob- 
tain a continuous absorption spectrum due to the 
photo-ionization of H™~ ions. 


(b) Application of the results to the electron 
affinity spectrum of halogen atoms 


The halogen atoms are known to have great 
affinity for electrons and their negative ions have 
a stable external configuration similar to the inert 
gases. If the negative ions are formed by an ele- 
mentary process of combination between neutral 
atoms and electrons, the emission spectrum of 
halogen atoms should be very roughly similar to 
that of hydrogen. On this qualitative basis we 
expect that the conclusions obtained for hydrogen 
may partially apply to the continuous spectrum 
of halogen atoms. 

In his recent experiments, Oldenberg found no 
trace of an electron affinity spectrum for the 
iodine atom under widely varying experimental 
conditions. The possibility of having a distribu- 
tion of intensities such that the maximum inten- 
sity is far separated from the spectral limit be- 
comes less for the iodine atom because of the 
slower change of the v‘-factor at a larger electron 
affinity (about 3.2 volts). It seems that the only 
other alternative is to assume that the intensity 
of emission is too faint to be experimentally de- 
tectable; or in other words, the combination be- 
tween neutral iodine atoms and electrons with 
simultaneous emission of radiation is compara- 
tively an improbable process—a _ conclusion 
reached by Oldenberg. 

Strictly speaking, however, we are not justified 
to make quantitative comparisons between the 


18 F, L, Mohler, Phys. Rev. Sup. 1, 2 (1929). 
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two cases for at least two reasons. First, the 
present calculations for hydrogen are based 
upon very rough wave functions in the continu- 
ous state; hence the solution is far from being 
rigorous and exact. The results may be expected, 
at best, to be only approximately correct for 
hydrogen and must not be applied directly to the 
halogen atoms, which presumably have quite 
different potential functions. Secondly, the dif- 
ference in magnitudes of electron affinity between 


the halogen atoms and hydrogen is quite large— 
the halogen atoms have electron affinities in the 
range of 3-4 volts against about 0.7 v for hydro- 
gen. The influence of the magnitude of electron 
affinity comes in when one deals with the v*- 
factor in the expression for intensity. Under these 
circumstances we cannot expect to have a simple 
quantitative relation between the two problems. 

The writer is very much indebted to Professor 
E. C. Kemble for invaluable advice. 
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Note on Kowalewski’s Top in Quantum Mechanics 


Otto Laporte, University of Michigan 
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It is well known that in classical mechanics algebraic 
integrals of the top equations only exist in the cases of 
Euler (asymmetric top, no electric moment), of Lagrange 
(symmetric top, electric moment parallel to the axis of 
figure) and of Kowalewski. The latter case is that of a 
symmetric top whose two equal moments of inertia are 
twice as large as the third (A = B=2C) with an electric 
moment perpendicular to the axis of figure. The quantum 
mechanical analogue to Euler’s and Lagrange’s cases being 
well known, Kowalewski’s case was tried. If 3, y, ¢ are 


the Euler angles, Qr=Q:+7Q2, Q1r=Q:—iQ» linear com- 


. binations of the momenta around the principal axes and 


U=Qr+4Cz sin de-t¥, 
then Kowalewski’s integral becomes: 
UU*+ U*U+4h?(Q1011 + 01101) = Diag. Matrix, 


which differs from the classical result by the symmetriza- 
tion and by the last term proportional to h?. 





I. INTRODUCTION 


HERE are two cases where the classical 

equations for the rotation of a rigid heavy 
body around a fixed point may be integrated 
readily by separation of variable of the Hamil- 
tonian in question: Firstly the case of Lagrange, 
the case of a symmetric top (A = B¥C) with the 
center of gravity on the axis of figure. Secondly 
the case of Euler, the case of a totally asym- 
metric top (A# BHC) with the center of gravity 
coinciding with the fixed point of rotation. How- 
ever in 1888 Sonya Kowalewski' succeeded in 
solving another case of top motion, namely that 
of a symmetric top whose two equal moments of 
inertia are twice as large as the third (A = B= 2C) 
and whose center of gravity lies in the plane of 
the two equal moments, or perpendicular to the 
axis of figure. Kowalewski only was able to 
integrate this case because she discovered (in 
addition to the energy and momentum integrals) 
a third algebraic integral which brought the 
number of first integrals up to three, as in the 
cases of Euler and Lagrange. 

It is of interest to note that in 1907 it was 
proved by Husson? that a third algebraic integral 
only exists in the aforementioned cases of Euler, 
Lagrange and Kowalewski or in special cases of 
these. 


1S, Kowalewski, Acta Mathematica 12, 177 (1888). 
2 E. Husson, Ann. Fac. Science Toulouse 8, 73 (1906). 


The quantum mechanical analogue of La- 
grange’s case, namely that of a symmetric top 
molecule with the electric moment parallel to the 
axis of figure was treated by Reiche and Rade- 
macher,® by Kronig and Rabi,+ and by Manne- 
back.® Euler’s case, the asymmetric top, was 
treated quantum-mechanically by Wang,® by 
Kramers and Ittmann,’ and by Klein.’ In the 
present paper we show that the quantum 
mechanical analogue of the Kowalewski top also 
possesses a third algebraic integral. Thus it is to 
be expected that a complete solution of the 
problem is possible in spite of the inseparability 
of Hamiltonian of the problem. 


Il. THE FUNDAMENTAL EQUATIONS OF 
KKOWALEWSKI’s Top 


For the sequel the reader is referred to the 
Leiden dissertation of H. B. G. Casimir® where 
more detailed derivations of the equations of 
rotation in quantum mechanics are given. 

The independent variables are the Euler angles. 


3F. Reiche and H. Rademacher, Zeits. f. Physik 39, 
444 (1926); 41, 453 (1927). 

*R. deL. Kronig and J. J. Rabi, Phys. Rev. 29, 262 
(1927). 

6 C, Manneback, Phys. Zeits. 28, 72 (1927). 

6S. C. Wang, Phys. Rev. 34, 243 (1929). 

7H. A. Kramers and G. P. Ittmann, Zeits. f. Physik 53, 
553 (1929); 58, 217 (1929); 60, 663 (1930). 

8Q. Klein, Zeits. f. Physik 58, 730 (1929). 

9H. B. G. Casimir, Leiden Dissertation, 1931. 
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The notation of Whittaker’s Analytical Dynamics 
is used so that if OX YZ is a right-hand system 
fixed in space and Oxys a right-hand system 
fixed in the body and if OK is the nodal line then 


{s0Z=3; +«YOK=9; +«yOK=y. 


It is not necessary to write down all nine 
direction cosines (cf. Whittaker, page 10), all we 
need are the cosines of the angles formed by 
the vertical Z 


with x: y,;=—sin # cos y, 
with y: ye= sin # sin y, (1) 
with s: y3= cos v. 


As axis of figure we choose the z axis. Let the 
electric moment be in the direction of the x axis, 
so that the potential energy is 


V=yvyi=—xsIin 0 cos y. (2) 


As momenta we have Q;, Qe, Q3, the momenta 
around the axes of the principal moments of 
inertia. Now it is essential for the following, 
that both for the Q and for the y we always use 
the following complex combinations: 


Or =Qi+7Q2; 
On = Qi—1Q2; 


As operators upon functions of the Euler angles 
the Q look as follows"®: 


0 ra) d | 
(cos v———)+i—|, 
dy de av! 


v1 =Yitiy2=—sin de~*¥, 
(3) 


yu = 71—ty2= —sin det*¥, 





aT | 1 
Q; = —the~*¥ 
sin 3 


Qn = theti¥ 





a) a) _¢ 
(cos i ——— }—71—};, (4) 


sin 3 ay ay av) 


0 
Q; =—ih—. 
oy 


From these formulae the exchange relations of 
the Q with each other and with the y may be 
derived ; or they may be obtained by generalizing, 
in the well-known way their Poisson brackets 
of classical mechanics. Using the abbreviation of 
Born and Jordan" 


1° Casimir, p. 57, Eq. (18). 
1M. Born and P. Jordan, Elementare Quantenmechanik, 
Berlin, 1930, p. 23, Eq. (7). 


(i/h)(xy—yx) =[x, y] (5) 


we put down, for future reference, the following 
table of exchange relations:'? 


LOr, Q3]=7Q1, (a) LQ: v.J=0, (h) 
LOn,QsJ=—iQn, (b) ([Q1, vs]=iv1, (i) 
LQr, QxJ=—27Q;, (c) [Qu,vs]=—tyn, (k) 
LOr?,QsJ=2707, (dd) [03 nJl=—-im, (1) 
LOn’, QsJ=—27Qn’, (e) [Qs,vnJ=iyvn, (m) 
[Qr0n, Qs]=0, (f) [Q1,mjJ=—2iys, (n) 


[vis 7: ]=0, (g) (On, v1] =2i7s, (o) 
CO?, vn] =—2t(Qrvs+73Q1), (p) 
[On?, v1] = +27(Quys+ysQn). (q) 


If we put A = B= 2C we obtain as Hamiltonian": 
I= (1/8C)(QrQun+ QnQn 
+(1/2C)Q2+3u(vi+71). (6) 


H put equal to the diagonal matric E is the first of 
the three algebraic integrals. The second is due to 
the fact that 7 does not contain the Euler angle 
¢ explicitly: thus 


Pe = S(yunQr +7101) +7303 = | diag. Matrix. (7) 


III. KowALEwsk1's INTEGRAL 


We now have to work with the Euler equations 
of the top which are obtained from the con- 
sideration, that the time derivative of any 


matric function f is given by f= (LH, f]. We then 
obtain, using repeatedly the exchange relations: 


4001 = —i(QsQ1 + 010s) +41 Cus, 

4COn = +7(0;:0n+QnQs)—4iCuys, (8) 
4CQs= 2iCu(yi— 101), 
4Cyr = — 2i(Qsyr +0) +i(Orvs +7300), 

4Cyn = +2i(Qsyn+nQs)—i(Quys+72n), (9) 
4Cys=—}i{ (Qrm+7Q) — (Quu+n0n) I. 


2 Casimir, p. 44, Eqs. (1) to (6). 
8 Casimir, p. 45, Eqs. (10). 
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Except for the symmetrizations Eqs. (8) and (9) 
are exactly the same as in classical mechanics. 
We multiply the first Eq. (8) with Q; from the 
left and then again with Q; from the right and 
thus obtain the time derivative of the square of 


Q1: 
4C(Q;")* = —71(07Q3+2010301+ 0301") 
+41Cu(Qryst+ysQ1). 
The last term of this equation is eliminated by 
means of (9.1). Applying exchange rule (a), (d) 
and (Il) we have: 
C(Qr?—4uCy1)° 
= —103(Qr?—4uCy1) —th(Qr— 2nCy1) 
= —1(Q—4uCy1)03+2h(OP—2uCy1). 
We shall use both forms later on. We now treat 


the second Eq. (8) and the second Eq. (9) in a 
completely analogous manner and obtain: 


(10.1) 
(10.2) 





OTTO LAPORTE 


C(Qrn?—4uCymn)° 
= i03;(On?— 4uCyr) —_ th(Qn?— 2uCyir) 
= i(Qn’?— 4uCyn)Q3+ th(Qy?— 2uCy 11). 
We introduce the abbreviations: 
OP—4uCyr = U, 
QOn?—4uCyn = U*, (12) 
UU*+ U*U=T. 


(11.1) 
(11.2) 


We compute 

T= UU*+UU*+U*U+ U*U. 
It now simplifies the calculation of JT con- 
siderably if in the first term of this expression one 


uses (10.2) for U and in the last term (10.1) for 


U. Similarly in the second term of (13) we will 
use (11.1) for U* and in the third term (11.2) for 
U*. We then have: 


—iCT = {—(QP?—4uCy1) Qs +h(Qr—2nCy1) } (Ou?—4uCrn) 
+(Qr?—4uCy1) {Q3(Qn?—4uCyn) — h(Qn?— 2uCyi) } 
+ {(Qu?—4uCym) Qs +h(Qn?— 2uCym) } (QP? — 4uCr1) 
+(Qn?—4uCym) {| — Q3(Qr?—4uCr1) —h(Qr—2uCy1) }. 


Evidently all terms which are not proportional to / cancel. We are left with: 


CT =h?{[QP—2uCy1, On?—4uCyn ]+(On?—2uCyn, Or—4uCri J}. 


The terms containing Q in the fourth order 
cancel each other in the two exchange brackets. 
So do the terms containing products of y. We 
then have: 


= —2h?u{(Qr, vn ]+COn’, v1 J}. 


Applying exchange rule () and (q) we finally 
get: 


T =4ih?u(Qrvst+7301— Ouv3— 730). (14) 


Leaving this equation for the moment we now 
return to the Euler Eqs. (8). We propose to 
multiply (8.1) with Qy first from the left and 
then from the right and similarly (8.2) with Q; 
both from the left and right sides. These four 
equations are then added: 





4C(Q10n + QOnQ1)” 
=1(— Q1Q301— On Q103s— O3Q10n | 
— Q103001 + 010300 + Q10n?s + (15) 
+Q30n01+ OnQs01) | 
+4tuC(Qny3st+y3sQu— O1y3s— ¥3Q1)- | 
In the first parenthesis on the right side the first 
and last as well as the fourth and fifth term 


cancel. The remaining four terms of the first 
parenthesis may be written: 


[LOr, Oud, Qs] 


which vanishes because of (c). Combining (14) 
and (15) we find: 


T= — 4h?(QrQn + QnQr)* 


or integrating and using (13): 
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UU*+ U* U+4h?(Q:0n+Q0nQr) = Diag. Matrix. 


(16.1) 


Taking into account (12) and (3) one may write this: 


(QP? —4uCy1) (On? —4uCym) + (Qn?—4uCym) (QP? —4uCy1) +447(Q10n + OnQr) = Diag. Matrix. 


or 


(QP°— Qe? +4uC sin 3 cos ¥)?+(Q102+020:—4uC sin 3 sin y)?+4h?(02+0,*) = Diag. Matrix. 


Eqs. (16) represent the quantum mechanical 
analogue of the Kowalewski Integral, which in 
classical mechanics simply is: UU*=const. It 
is interesting that the simplest quantum-like 
generalization, namely symmetrization of U and 
U* is not sufficient to obtain an integral; the 
term proportional to /? has to be added. Nor is 
it possible to embody this correction term in the 
definition of U as one might attempt in order to 


(16.2) 


(16.3) 


preserve the classical form as closely as possible. 

Eqs. (6), (7) and (16) form the three first 
integrals of the problem. By means of these (or 
rather their classical counterparts) Kowalewski 
showed that the problem could be completely 
solved and reduced to hyperelliptic integrals. 
It remains to be seen whether or not a 
similar reduction is possible also in quantum 
mechanics. 
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The theory of gases in nonstationary states, given by 
Lorentz and Enskog, is generalized for the quantum 
statistics to give the hydrodynamical equations and the 
distribution function in first and second approximation, 
and formal expressions for the viscosity and heat con- 
ductivity coefficients. These results are valid for all sta- 
tistics and for all degrees of degeneration. Two essential 
points contribute to this generality: (a) Exact expressions 
independent of statistics and degeneration can be given for 
the coordinate and time derivatives of the coefficient A of 
the equilibrium distribution function in terms of the 
pressure and temperature gradients and time derivatives, 
though a closed expression for this coefficient as a function 


of vand T is known only in limiting cases; (b) The function 
W of the general equation of state for ideal gases in all 
statistics 


pu =(RT/M)W(v'T) 


is adiabatically invariant. Numerical values of the viscosity 
and heat conductivity coefficients, which should come out 
of the formal theory on the introduction of special assump- 
tions about the molecular forces, have not yet been 
obtained. It is our hope that these results, when found, 
may furnish an experimental test of the existence of 
Einstein-Bose statistics in real gases, as is required by 
theory. 





INTRODUCTION 


ar 


1 


One may distinguish broadly two methods of 
treating the problem of viscosity and heat con- 
ductivity of an ideal gas. The first of these is an 
inexact method according to which an equilib- 
rium distribution function is considered to be 
descriptive of the state within the gas, subject 
only to the condition that the five quantities, the 
density, the temperature, and the three com- 
ponents of mass velocity are functions of posi- 
tion. The conception of the mean free path is 
fundamental in this theory. 

According to the second and more nearly exact 
method, which is due mainly to Maxwell and 
Chapman on the one hand, and to Lorentz and 
Enskog! on the other, one recognizes the fact 
that the above assumption regarding the equi- 
librium distribution function is not valid, since 


of 





Ox Orr 


the latter does not fulfill the fundamental Boltz- 
mann continuity equation. We propose to gen- 
eralize this method so as to include the Einstein- 
Bose and the Fermi-Dirac as well as the \ax- 
well-Boltzmann statistics. 


§2 

In this generalization we will follow Lorentz 
and Enskog, rather than Maxwell and Chapman. 
An understanding of the method used may be ob- 
tained from a consideration of the electron theory 
of metals of Lorentz.? The essential points in this 
theory follow. 

The fundamental Boltzmann equation in this 
case, where one assumes the existence of a steady 
state, the dependence on one coordinate only, 
and elastic collisions between the electrons and 
fixed ions with no interaction among the elec- 
trons themselves, takes the form: 


of 
ts—t+X —=nRr f {f(re'ry'r 2!) —f(r2tyr2) } cos ddQ, (1) 


where # is the angle between the direction of the line of centers in collision and the electron veloc- 
ity r; dQ is the element of solid angle within which the direction of the line of centers must lie; 7 is 


1 Maxwell, Collected Works II, p. 1; Lorentz, Theory of | (1916); Enskog, Kinetische Theorie der Vorgiinge in missig 


Electrons, Note 29, p. 266, Ges. Abh. I, p. 72, Vortrige 
tiber die kinetische Theorie der Materie und Electrizitét, p. 


185; Chapman, Phil. Trans. 216 A, 279 (1915); 217 A, 115 


verdiinnten Gasen, Dissertation, Upsala, 1917. Also see 
Jeans, Dynamical Theory of Gases, Chapters VIII and IX. 
2 Lorentz, Reference 1. 
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the number of ions per unit volume; R is the sum 
of the radii of the electron and ion; and X is the 
acceleration due to the external electric field. 
Primes denote quantities which are to be taken 
after collision. 

In the determination of the distribution func- 
tion which fulfills this equation one can distin- 
guish the following steps: 

(a) The principle of solution. One assumes: 


f=fA1+¢), (2) 


where f is the equilibrium distribution func- 
tion, and ¢ is a small perturbation term. Intro- 
ducing (2) into (1), and retaining only first order 
quantities, one obtains, since f® is a function 
only of r: 


of of™ 
fs ——_ +t X —=nRf™ [te —¢) cos ddQ., (3) 
Ox Ore . 


(b) Conditions imposed on ¢. Eq. (3) is in an 
inhomogeneous integral equation of the second 
kind with a symmetric kernel. Because the homo- 
geneous equation has a solution (namely ¢ 
=const.), the necessary condition that the in- 
homogeneous equation has a solution is that the 
left member of the equation is orthogonal to ¢1. 
This condition is identically fulfilled. To make 
the solution of (3) perfectly definite, one may im- 
pose on ¢ the auxiliary condition: 


[ eio=o. 


One finds that this implies: 


(4) 


[0 edradrgdr.= [Hs edradrydr = 0, 


e 


which means physically, that the density and 
temperature are determined by the mean values 
formed with the equilibrium distribution function 





Of 


dl 


+D(N)= f dor f ge(og)ae IL AA+F OP A+) fill +of)(1+efr’) }, 
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just as in the equilibrium state. 


(c) Reduction of the left member. Introducing 
f@ =A (x)e~* (z)r2 


where h(x)=m/2kT(x), one obtains for the left 
member: 


r,(d In A/dx—r'dh/dx—2hX)Ae~*”. (5) 
(d) “‘Ansatz” for ¢ and its verification. The form 
of the left member suggests the solution: 


g=rzx(r), 


which fulfills condition (4). Introducing this ex- 
pression into the right member of (3), one finds 
upon integration that the latter becomes propor- 
tional to 7,, giving the final result: 


1 ( dh dinA 
r—— 
an R? dx dx 


from which all further conclusions regarding the 
heat and electric conductivity follow.* 

This description is intended to show only the 
analogy between the method of the Lorentz 
theory of metals and the generalized Enskog 
theory of gases. Points (a) to (d) of the Lorentz 
theory given in §2 correspond to points §5 to §8 
in the theory of gases generalized to include the 
quantum statistics which is to follow. A strict 
analogy does not exist. One would have to con- 
sider in the Lorentz theory also nonstationary 
states, and extend the theory of gases to include 
also mixtures in order to see the former as a 
special case of the latter. 








1 
x(r)= +2hx),, 


r 


THE GENERAL HYDRODYNAMICAL EQUATIONS 
$3 

For Einstein-Bose and Fermi-Dirac gases in a 
nonsteady state, one may write the general 
Boltzmann equation in the form: 


(6) 


3A well-known defect of the Lorentz theory is that it | following difficulty. If one attempts to find the next 


gives no account of the Joulian heat, since only elastic 
collisions between the electrons and immovable ions are 


considered. This is probably the reason, also, for the | 


approximation to the distribution function, one obtains an 
integral equation of the same type as (3). But now the 
left member is not orthogonal to the solution of the homo- 
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where 


D(f)=raldf/axe)+Xaldf/dre). (7) 


In these expressions 7;(7,7,7.) are the three com- 
ponents of the absolute velocity, X ;(X YZ) are 
the three components of the external force per 
unit mass, and the usual summation convention 
has been used.‘ The function w(dg)dQ represents 
the effective cross section for a collision which 
changes the direction of the relative velocity g of 
two molecules by the angle #, such that this 
velocity after collision g’ lies within an element of 
solid angle dQ= sin ddddy, where ¢ is the azi- 
muthal angle about g. The subscript 1 denotes 
functions and variables pertaining to the second 
molecule in the collision, over which the integra- 
tion takes place, and the primes denote functions 
and variables which are to be taken after the col- 
lision. The differential d¢ has the usual interpre- 
tation: 


do= VG(m/h)*dr,dr,dr., (8) 


where G is an eventual weight-factor, and V is 
included for dimensional reasons, but is to be 
considered as having the value unity to agree with 
the meaning of f as giving the number of particles 
per phase cell in coordinate-momenta space of 
which the space part has unit volume. 

Eq. (6) includes (namely for 9=0) the Boltz- 
mann equation of the classical statistics® as a 
special case. It differs from the latter in the Ein- 
stein-Bose statistics (@= +1) and in the Fermi- 
Dirac statistics (@= —1) only in the “collision” 
term of the right member. The essential points of 
difference in this term are: (1) the use of the 
appropriate ‘‘Stoszzahlansatz,” and (2) the 
necessity of determining the transition probabil- 
ity function, w(g), quaantum-mechanically, and 
of taking into account the identity of the mole- 
cules in this determination, as was first done by 
N. F. Mott in his theory of collisions between two 
electrons.® That the “streaming”’ terms of the 
left member need not be changed in the quantum 
theory, has been proved in general by Nordheim 


geneous equation as a consequence of which the in- 
homogeneous equation has no solution. In the theory of 
gases one does not encounter an analogous difficulty. 

* Greek letters are used for dummy indices. 

§’ Compare e.g., Jeans, Dynamical Theory of Gases, p. 210. 

* Proc. Roy. Soc. A126, 259-267 (1930). 


and Kikuchi.’ As usual, the principle of micro- 
scopic reversibility has been used in the deriva- 
tion of the collision term in Eq. (6). 

Just as in the classical case® one can derive from 
Eq. (6) a general transport equation for a func- 
tion F(xyzr,r,r.t). Specializing F to be succes- 
sively m, mr,, mr,, mr., and $mr’, one obtains the 
general hydrodynamical equations: 








dp/dt+ p(du./IdXa)=90, al 
p(du;/dt) = pXi—OPia/IXa, b 
(9) 
d(Q/p) a 
p + _ — PapS af; Cc 





dt OXa 


4 


where u;(uvw) and q;(q.9,g:) are the components 

of the mass-velocity and heat current, respectively: 
(a) w=ri, (b) gi=1pUiU? (10) 

and where, furthermore 

W=UP/+U/7+U? (11) 


Us="i- Ui; 


denotes the relative velocity. The functions p;; 
and S;; are the components of the stress and the 
rate of pure strain tensors, respectively: 





— Ou; Ou; 
(a) py=pU,U;, (b) Sum 3 +=), (12) 


Xj Ox; 


Finally, Q denotes the energy of heat motion per 
unit volume: 


Q=3pV?=3 (pret Pyyt P22) =3p/2 (13) 
and d/dt denotes the total time derivative: 
d/dt=0/dt+u.(0/dxq). (14) 
$4 
One obtains the hydrodynamical equations in 
first approximation by calculating the ten mean 
value quantities, (10b), (12a) and (13) with the 


aid of the equilibrium distribution function f™. 
Since the latter depends only on U®, one obtains’: 


Pii= bis, 

7 Nordheim and Kikuchi, Zeits. f. Physik 60, 652 (1930). 

8 Compare Enskog, reference 1, pp. i4-19; Jeans, ref- 
erence 1, Chapter IX. 

*The equation: Q=3p/2 or pu=2e/3, where « is the 
energy per unit mass, and v is the specific volume, follows, 
also, from an application of the general virial theorem to 
an ideal gas. It is, itself, often called the virial theorem. 


gi= 0, 
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whereupon (9) reduces to: 
dp/dt+ p(du_/dxa)=0, al 
p(du;/dt)= pX;—(dp/dx,), b 


d *) 2 p Ota 
—{-)=--- ‘ c 
di \p 3 pOXa 








J 


The energy equation (15c) can be put into an- 
other form by introducing the equation of state: 


b/p=(RT/M)W(Tp-5), (16) 


which holds in all statistics. M denotes the 
molecular weight, p is the density, and W is a 
function of the argument Tp~! which is a different 
function for the different statistics, but is in each 
case an adiabatic invariant.” The latter property 
makes it possible to remove W from the differ- 
ential operator d/dt when (16) is introduced into 
(15c). Thus the energy equation becomes: 


dT /dt= —2T(du,/dxa). (15c’) 





where 
x= o/(1 +68) 


1(x)= f don f gorda fof +f) (1+ Of") (xxx — x). 


One can show that Eq. (18) has the following 
properties: (a) It is reducible to the form of an 
inhomogeneous integral equation of the second 
kind"; (b) the kernel of this equation is symmet- 
ric; (c) the kernel is invariant with respect to 
rotations in U ;-space. Furthermore, the operator 
I has certain properties, which will be required 
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In Note 1 it will be shown that an expression for 
the velocity of sound in first approximation, that 
is, neglecting viscosity and heat conductivity, 
follows immediately from (15). 


THE NONSTEADY STATE DISTRIBUTION FUNCTION 


§5. The principle of solution of the Boltzmann 
equation 

In order to proceed with the theory of gases, 
one needs to know the deviation of the actual 
distribution function from the equilibrium state. 
For this determination one must obtain a solution 
of Eq. (6). As in the Lorentz theory outlined in 
$2, one assumes a solution of the form: 


f=fO(1+e9), 


which, upon introduction into Eq. (6) and sim- 
plification with the aid of the principle of de- 
tailed balancing," yields the inhomogeneous in- 
tegral equation: 


af /at+D(f) = —1(x) 


(17) 


(18) 


(19) 


(20) 





later. Using the notation: 
[F, G]= f FI(G)d¢, 


one finds by the same reasoning as is used in 
deriving the 17-theorem: 


[F, G]= tf dof dds f ew(ogaa- sof (1 +f) (1+ 0/1) -(F+F,— F’ — Fy')(G+G,—-G'—-G,’), 


 Uhlenbeck and Uehling, Phys. Rev. 39, 1014 (1932). 

"This principle follows, as is well known, from the 
previously assumed principle of microscopic reversibility 
and a generalized H-theorem. Comp. Nordheim, Proc. 
Roy. Soc. A119, 689 (1928); Pauli, Sommerfeld Fest- 
schrift, p. 30. 

" This reduction is accomplished by a series of trans- 
formations of variables, in a manner analogous to that 


used by Enskog in the classical case (Reference 1, pp. 
140-148). 


(21) 
from which follows: 
[F, G]=(G, F], a 
[F, G+H]=(F,G)]+[F, 7], b (22) 
[ F, bG]=0b(F, G], c 
[F, F]=0, d 


where b is a function independent of all variables 
of integration. 
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§6. Conditions imposed on ¢ 
Consider the homogeneous equation following 
from (18): 
I(y)=0. 


Multiplying with Yd@ and integrating, one ob- 
tains: [y, y]=0. From the properties of this 
operator this equality exists only if: 


¥t+yi-Y—y'=0. (23) 


Eq. (23), according to the number, momentum 
and energy conservation theorems, has as its 
solution a linear combination of the functions": 


Vi=1; Yo= rz; ¥3,=7y3 Way = P23 Y5=r. (24) 


In order, therefore, that the inhomogeneous Eq. 
(18) may have a solution, the left member must 
be orthogonal to the five functions (24). As conse- 
quences of this orthogonality one finds that the 
five parameters p, 7, u, v and w on which the 
equilibrium distribution function f® depends, 
must, for the purpose of solution of Eq. (18), be 
functions of x, y, z and ¢ which fulfill the hydro- 
dynamical equations in first approximation (15)." 
To make the solution of Eq. (18) definite one 
may impose on ¢ or x the five conditions: 


(18a) 


f Vif edo= f VfL + OF) xdb=0, 
mt, 2,--- 5) (5) 


which interpreted physically express the fact 
that the parameters p, 7, u, v and w of the 
equilibrium distribution function may be inter- 
preted physically as the density, temperature, 
and the three components of mass velocity, re- 











fo Pe 


Ae-*L dt T dt oc OX 


Now one may prove, as is done in Note 3, that in 
all statistics and for all degrees of degeneration 
the following relations exist: 


dinA dW idp 31 dT 
+0/ a eo —) a 
2T dt 


> 


dt dt 


dln A p (- Op 51 ~) ; 
Seer f<e—<“oc-— =e 2. ) 
Ox; mnkT \pox; 2T Ox; ; 


Since W is adiabatically invariant, and since it 





L (29) 








rdT 1 dInA OUe 
+— —+— i ~2e( X.— ) 
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spectively, just as in the equilibrium state. One 
can show, in a manner analogous to that used by 
Enskog in the classical case (reference 1, p. 27- 
31), that if one has determined the arbitrary con- 
stants in the general solution of Eq. (18) such 
that the conditions (25) are fulfilled, at the time 
t=0, they will remain fulfilled for all values of the 
time. 


$7. Introduction of new variables and reduction 
of the left member 


One may introduce the equilibrium distribu- 
tion function: 
1 
{= 26) 
(1/A) exp (m/2kT) U?—6 





into the left member of Eq. (18). The five param- 
eters, p, 7, u, v and w on which this function de- 
pends are defined by the equations: 





p=m [ f%d6, 
1 

u;= - f mr fae, 
, 

(27) 

3p 3RT 

e=- —=- — W(Tp-3) 
2p 2M 
1 1 

=— f- mU?f©dg. J 

p 2 


Transforming, also, to dimensionless variables 
E(émo) where §:=(r;—uc}; P=P+n°+e and 
c=m/2kT, Eq. (18) becomes: 


r oT | Ole 
}>———~—} + aaks 
T ax.) 








|- —I(x). (28) 


ot Oxg 





has been proved that a solution of Eq. (18) exists 
only if p and 7 fulfill the hydrodynamical equa- 
tions (15a) and (15c’) in first approximation, one 
finds that dln A/dt vanishes. Furthermore, one 
obtains from (15b) and (29b) the result: 


‘8 That (24) are the only solutions of (23), one can 
prove strictly for spherically symmetric molecules (com- 
pare Jeans, Reference 1, p. 25), to which this theory is 
restricted. 

4 See Note 2 for the proof of this result. 
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is obtained.'® 

Dimensionless variables and functions will be 
introduced, also, into the operator I(x). One may 
decompose the transition probability function in 
the following arbitrary manner 


om 3 


w(dy) = fs?-a(y)-u(dy), 
where 
g=(2/c)}, 


sis a constant having the dimensions of a length 
and in the case of spherical molecules being the 
actual collision radius, o(y), is a dimensionless 
function, and u(y) is a dimensionless probability 
function representing the scattering into an ele- 





where 


1 
I(x) = ol [ d2-y0"(a)u(oy) (xtx1— x! — x1) ff (1 + Of) (1+ Of). 


f™ (11 £4od@T 5 p 7 Oa OUg 
{= = & — #--—-)+(et—— ian) ( +—*)}=-100 
Ae" lei'T Axe 2 nkT 3 Oxg OX 
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Introducing also the value of d7°/dt from (15c’) 
the equation 





(30) 





ment of solid angle per unit angle divided by the 
average scattering per unit solid angle. In the 
case of a classical collision between two spherical 
molecules the functions o and u are simply unity. 
In the case of a classical collision between mole- 
cules repelling inversely as the vth power of the 
distance, where y>1, u is independent of y, and 
oa depends on y and on the temperature 7, 
Finally introducing the variables £;, the differ- 
ential becomes 


do,= VG(m/h)8c~ dw; dw = dé,dnid¢\. 
In these variables and functions 


I(x) = (2!s?A2 VG/c?)(m/h)3I'(x), (31) 


(32) 


One finds that the bracket expression [ F, G] defined with the operator J’ has the same properties 


(22) as the expression defined with the operator J. 


§8. The ‘‘Ansatz” for x and its verification 


Since it is linear, and the integration in J’(x) is over only the velocity components £;, Eq. (30) 


may be decomposed by choosing for the solution: 


c (- 1 0T 


2's®AVG \m/J T Ox, 


where 7; and 7;;=7;; are nine new unknown 
functions, which fulfill respectively the nine inte- 
gral equations: 


fo 5p 
§; (#-= )=rw0, a 
A*e-” 2 nkT 


(0)2 7 
— (s8-— iu )= Tas). b 
A*e-” 3 


~ 





> (34) 








In a manner analogous to the Lorentz theory 


%In the classical case, p=nkT and f\=Ae-™; the 
left member of this equation then agrees completely with 
the result of Enskog. (Reference 1, Eq. 53.) 


— —_ 7,— 





C h\* /0uq Oug 
Seal) 
2: AVG \m OxXg OXe 





(33) 





(see d in §2), one can show that 7; and 7;; have 
the forms: 


i= §:x0(7), a| 
(35) 
mij= (tits — 726;;/3) x0(17). I 


The problem is reduced in this way to the deter- 
mination of the two unknown functions xa(r) 
and x,(r), The procedure differs, however, from 
that in the Lorentz theory in that it is not pos- 
sible in this case to obtain closed expressions for 
these functions. They may be obtained only by a 
method of successive approximations after defi- 
nite assumptions regarding the nature of the 
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molecular interactions have been introduced; i.e., 
after the functions o(y) and u(dy) have been 
specified. 

One can prove that the functions 7; and z;; 
have the forms given by (35) as follows. One veri- 
fies that the left members of both (34a and b) are 
orthogonal to the five functions y; of (24), and, 
therefore, orthogonal as well to the functions 


Vi=1; eo=§; vVs= 0; We= S53 Vs=7*, (24a) 


which as before represent the solutions of the ho- 
mogeneous equation. The forms of 7; and 7;; may 
be made definite by imposing on both the five 
conditions: 


f Vf (140%) xdw=0 (36) 


as consequences of which x will fulfill conditions 
(25). As further consequences of conditions (36) 
Eqs. (34a and b) each have one and only one 
solution. Finally, one uses the fact that Eqs. 
(34a and b) have properties analogous to those of 
Eq. (18) discussed in $5; namely: (a) they are 
linear integral equations of the second kind, (b) 
the kernel is symmetric, (c) the kernel is in- 
variant with respect to rotations in &,-space. In fact, 
one can show that: 


I'(n)=R(1) (E+ f dw'n(&!)K(r1'e), (37) 


where ¢ is the angle between the two vectors 


cher? hry? 1 dT 
Peng (-) won moe Sake) — 
2's?A2VG \m/ T Ox 





f=f 1 - 


This expression reduces to that of Enskog" for 
the Maxwell-Boltzmann statistics in which: (1) 
the equilibrium distribution function is 
f=Ae-®; (2) xa(r) and xo(r) are determined by 
Eqs. (34a and b) with p= nkT, 6=0, and classical 
values of o(y) and u(gy); (3) account is taken of 


16 A second method is by means of the resolvent of (37). 
One ...n show that the resolvent as well as the kernel is a 
function only of 7, 7’ and ¢. Writing down the formal 
solutions of (34a and b) and using (38) one finds again 
that the solutions must have the form (35). 

17 Enskog, Reference 1, Eq. (65), p. 41. 
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£; and £,’, which have the magnitude 7 and 7’, 
This is the essential point in the proof. One may 
proceed now in different ways. The most simple 
method perhaps is verification. This is accom- 
plished when one has proved that integrals of 
the type: 


f awees(or'o) = §F(r) (38a) 


and of the type: 


aw (seer—— iu) fore 
r 
= (s/-7 iu) Gl). (38b) 


The proof of these results can be given geomet- 
rically (see Note 4).'® 

The auxiliary conditions (36) are fulfilled by 
(35b) identically, and by (35a) when x,(r) ful- 
fills: 


f PFOA Of xalr)de=0. (39) 


It is always possible to satisfy this condition by 
adding to x.(7) a suitable constant. 


§9. The distribution function 

In terms of the unknown functions x.(r) and 
x»(7) the nonsteady state distribution function 
may be written with the aid of Eqs. (17), (19), 
(33), (35a and b) in the following manner: 


cf (0) er2 


<< _(-) (+ “*) ( r ) «) |. 40) 
Eat sesnathins 28 ba Xx | 40) 
252A? VG Oty Xe i i a 








the weight factor by which the A of Enskog’s 
equations differs from the A used here. 


FORMAL EXPRESSIONS FOR THE VISCOSITY AND 
HEAT CONDUCTIVITY COEFFICIENTS 


$10. The heat conductivity coefficient 


In order to obtain this coefficient it will be 
necessary to derive one of the components of the 
heat current. By using the definition (10b), the 
expression (40) for the distribution function, and 
transforming to dimensionless variables, one 
obtains 





Oe ———— 
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gz= —X(OT/dx), (41) 


where: 
m 


h=———_ — | fer t272x,(r)dw. (42) 
(2c)'s?A? T 

For the purpose of the numerical calculation of 

this coefficient, one may transform this expres- 

sion; first, by adding a term which is identically 

zero by virtue of the auxiliary conditions on 

xa(T), giving 


m ? er 5 p 
{eo — piotere( 2 —— ) Exel) 
(2c)'s?A? T 2 nkT 





and then, by introducing a simplification ob- 
tained when one substitutes for a part of the inte- 
grand its value given by the integral Eq. (34a) 
in xa(r). Thus, one arrives at the final result: 


m 


1 
A=—— pre LExa(7), Exa(7) ]. (43) 
(2c) 4s? T 


That is always positive follows from the 
property (22d) of the bracket expressions. 


§11. The viscosity coefficient 


In order to obtain this coefficient, one may de- 
termine either p,, or the deviation of p,, from p. 
Using the definitions (12a) and (13), introducing 
the distribution function (40), and transforming 
to dimensionless variables, one obtains in the 
latter case: 


3 
where: 
3m 7 - 7 
eo fs mer(e ar ) (© -<) xo(7)dw. 
2icts?A* 3 3 
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By using the integral Eq. (34b) for x», this can 
be written: 


2m 


(eZ). (e-2) 00} 





— 


from which one finds, as in the case of the heat 
conductivity coefficient, that the viscosity coeffi- 
cient is always positive. 

In the other case one obtains in the same man- 





ner: 
Ov OU 7 
Pry= -1(—+—), (47) 
Ox oy 
where: 
2m ae 
= | fer ea? ya r)dw 
(2c)*s?A?e 
2m 
= —— [Enxo(r), Enxe(7) J. (48) 
(2c) s? 


A partial integration of these expressions shows 
that (46) and (48) are identical. 


§12. The hydrodynamical equations in second 
approximation 


The preceding derivation shows that the stress 
tensor components have the same form in the 
quantum as in the classical statistics. Eqs. (44) 
and (47) for these components may be combined 
in the form: 

Pis=(P+(2u/3) (Ou e/IxXa)|bi;—2wSi;. (49) 
Introducing these expressions and Eq. (41) for 
the heat current components into the general 
hydrodynamical Eq. (9). one obtains the hydro- 
dynamical equations with viscosity and heat 








(45) conduction: 
dp/dt+ p(du,/dxa)=9, a) 
du; Op 1 Pua Oru; 2 OUq Op _ Op ; 
p ——= pX,—-—+- ——+u——_—- — — 42S ta ——, b| (50) 
dt Ox; 3 Ox\OXe OXeOXa 3 OXg OX; OXe > 
3 d(p/p) 29 oT 
Sf reenter gees (a —) = — £5.04 Del SeaSea— Sa. c 
2 dt OXa OX } 
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Ordinarily 4» may be considered as a constant, 
since the variations in T and p, on which yu de- 
pends,'® are generally small. When this is the 
case, the last two terms of Eq. (50b) vanish, 
and this equation then takes the usual form.’ 
Introducing Q=3/2 into the energy Eq. (50c), 
one finds that it gives the rate of increase of heat 
energy per unit volume due to: (a) conduction of 
heat from adjacent elements, (b) the work done 
per unit of time by the internal (scalar) pressure, 
and (c) the work done in the distortion of the 
volume element. 


§13. Conclusion 


The formal part of the theory is now com- 
pleted. To proceed further, and thus obtain 
explicit expressions for the viscosity and heat 
conductivity coefficients as functions of the 
temperature and density, one must introduce 
special assumptions regarding the molecular 
forces. With the introduction of these assump- 
tions, the transition probability function w(dy) 
may be determined from the quantum-mechan- 
ical collision theory, whereupon the functions 
o(y) and u(dy) of the integral equations are 
completely specified. The integral Eqs. (34) may 
be solved then numerically. Enskog has de- 
veloped a method of solution in the classical 
case which may be extended, also, to the quan- 
tum statistics. These questions will be considered 
in a following paper. 


Note 1 


In order to derive the equation for the velocity of sound 
from (15), one assumes that the conditions are such that 
the velocity components and the deviation of the density 
p from its equilibrium value are small quantities of the 
first order, and that no outside forces exist. Then the 
density may be written in terms of a condensation factor s: 
p=po(1+s), where po is its equilibrium value. Retaining 
only quantities of the first order (15a and b) become: 





8s /dt = —OuUg|dXa, (51) 
po(du,;/dt) = —dp/dx;=— K(ds/dx;), (52) 
af dua af 


fol 


If these conditions are valid for the five ¥; of (24), they 
must be valid, also, for ¥; equal to m, mUz, mU,, mU,, 
and mU?*. Substituting these values successively into this 


aUg dt “OXa 


18 Only in the classical case is » independent of p. 
19 See Lamb, Hydrodynamics, Fifth Edition, p. 546. 
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where K is the coefficient of volume elasticity p(dp/dp). 
Eliminating u; or s, one obtains the wave equation for 
$ or 4; with the velocity: 


C? =K [po = (dp/dp) p=p,- (53) 


Now assuming that the frequency is sufficiently high so 
that one may neglect heat conduction, this differential 
quotient may be taken under strictly adiabatic conditions, 
The first order energy equation gives such a determination 
of this expression. Under the conditions assumed it 


becomes: 
Oz--nO 
Xe ~ Qa at 
Now using Eq. (51), integrating, and solving for p, one 
obtains: 


(54) 


b= 4(b/p)op(1+2p/po), (55) 


from which follows: 


C? = (dp/dp) p=, = (5/3) po/po (56) 


in accordance with the thermodynamic calculation of 
(dp/dp) aaiav * 


Note 2 


From the theory of integral equations the conditions: 


Jvlafojat+OYyide=0 G=1,2,- 


must exist in order that the inhomogeneous Eq. (18) may 
have a solution. This follows, also, upon multiplying both 
sides of (18) by yid¢ and integrating. The right member 
then becomes [y¥i, x] which is zero according to (21) 
and (23). 

In Eq. (57) f represents the equilibrium distribution 
function in the six-dimensional phase space (x;, 7;). Since 
it is a function only of the relative velocities, U;=r;—ui, 
it is convenient to change to axes which move with the 
volume-element of the gas. The differential quotients then 
become: 


* 5) (57) 











of af OU a’ f 
“at att BU’ 
of a’f OUe af 
Oe; ” bee ~ 8s Ue" 
af af 
ars OU," 


and Eq. (57) may be written, omitting the primes: 


vf Ita 
Up, if axg 


af 


“50, (¢=1, 2, +++ 5). 


= |de= 0 (58) 





equation, and using (27), one obtains the first order hydro- 
dynamical Eqs. (15). 


27Compare Uhlenbeck and Uehling, Phys. Rev. 39, 


1014 (1932). 
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Note 3 
From the definitions (27) one obtains 
n=p/m=VG(2xmkT/h*)' Fy, a | 


(59) 
p=40=4ep=VGET(2emkT/h)'Fy, bf 


where F, is the Sommerfeld integral: 


1 nee) 


s=r?=(m/2kT) V*. 


mo 8s gtdz 


1 
F(A) = sre4 (1/A)e*—0 








The result: 


b/nkT = F;/F, (60) 


follows as an immediate consequence of these equations. 
Also, from the definition of F,, one obtains the relation: 


A(dF,/dA) = F,-1. (61) 


Considering » and T as functions of the time and the 
coordinates, one obtains from (59b) the equation 


EW ES 51 *). 


dt pdt 2T dt 
On the other hand, Eq. (61) gives: 
dF; ,dinA dF; .dinA 
d°)Cdt( (CA Tt 


Combining these results and using the equation of state 
(16), one obtains Eq. (29a). In a similar manner one 
obtains Eq. (29b). 


Note 4 


Interpret f(rr’e) as a density distribution in the £,’- 
space. (See Fig. 1.) This function will have rotational 
symmetry around the axis OP, the direction of 7, and it 
will depend only on this direction and not on the special 
choice of the coordinate axes £,;’. To prove (38a) we 
remark that the left member represents the coordinates 
of the center of gravity, which, because of the rotational 
symmetry, must lie on the line OP. Therefore these 
coordinates, and accordingly the left member of (38a) 
must be proportional to &;. Since f(rr’e) is invariant with 
respect to rotations of the &;’ axis, the proportionality 
factor must be a function only of r. 

In order to prove (38b) one may consider the integral 


J dio" § 6/9 (o"d) = Ais 
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as representing the moments and products of inertia of the 
density function f(rr’e) in the £;’-space. Because of 
rotational symmetry, the ellipsoid of inertia is a rotational 
ellipsoid about the axis OP. Taking O as the center of this 
ellipsoid, and p and gq as the lengths of the main axes, one 
may write for its equation 


OR*/P+RE/¢ =1. 


s ; 
/ 


/ 
‘P(E ns) 


—@e'ns') 








n’ 
Fic. 1, 


Designating the direction cosines of OP by \4(AiA2A3) one 
obtains OR = \,ta’, whereupon the equation of the ellipsoid 
becomes, since RO? = OQ*—OR? 


(rar+ af bas) fa’! -£o 


But the equation of the inertial ellipsoid is also 


Aasta’ tg’ =const. 
Therefore 
, , , , p* 
fdw £5 E;'f(rr Im Nt bij 





~ tit 7°54 


J ao'( ee! < bi) f(rr'e) ~ (est = iss) 


and the proportionality factor as before can depend only 
on rf. 
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Effect of Etching on the Relative Intensities of the Components of Double Laue Spots 
Obtained from a Quartz Crystal 


M. Y. CoLBy AND S1pon Harris, Department of Physics, University of Texas 
(Received January 20, 1933) 


Laue patterns obtained from a quartz crystal of dimen- 
sions 100 6X2 millimeters taken through the 2 millimeter 
direction are given. One face of the crystal was highly 
polished and the other face was etched. The pattern ob- 
tained when the etched face of the crystal was next to the 
photographic film exhibited a weaker inner component for 
spots close to the origin; the pattern obtained when the 
etched face of the crystal was next to the source of x-rays 


exhibited a stronger inner component for spots close to 
the origin. The pattern obtained when the etched face 
of the crystal was next to the photographic film and the 
crystal was oscillating piezoelectrically along the direction 
of the long dimension showed no appreciable difference 
from the pattern obtained with the crystal in this same 
position when at rest. An explanation of the recent results 
obtained by J. M. Cork is offered. 





NUMBER of investigations on quartz 

plates oscillating piezoelectrically along the 
direction of thickness of the plates caused by the 
longitudinal effect! have been recently re- 
ported.” 3.4.5» 6 The results obtained showed 
that the intensities of the Laue spots were greatly 
increased by the piezoelectric oscillations. J. M. 
Cork found that the Laue pattern obtained when 
both sides of the quartz plate were etched with 
hydrofluoric acid exhibited double spots with the 
inner components showing more intensity for the 
spots close to the origin. He found that the spots 
farther out exhibited components of more nearly 
equal intensity. One purpose of the present work 
was to further investigate this phenomenon. As 
Fox and Cork had failed to observe any effect on 
the lines obtained by the Bragg method caused 
by piezoelectric oscillation produced by the 
longitudinal effect, and as the authors’ have suc- 
ceeded in obtaining an effect on the lines obtained 
by this method caused by piezoelectric oscilla- 
tions in a long quartz crystal produced by the 


1The terminology here used is taken from Piezo- 
Electric Terminology, by W. G. Cady, The Inst. of Radio 
Eng. 18, No. 12 (1930). 

2 Y. Sakisaka, Jap. Jour. Phys. 4, 171 (1927). 

3G. W. Fox and P. H. Carr, Phys. Rev. 37, 1622 (1931). 

4G. W. Fox and J. M. Cork, Phys. Rev. 38, 1420 (1931). 

5C. S. Barrett and C. E. Howe, Phys. Rev. 38, 2290 
(1931); 39, 889 (1932). 

6 J. M. Cork, Phys. Rev. 42, 749 (1932). 

™M. Y. Colby and Sidon Harris, Phys. Rev. 42, 733 
(1932). 


transverse effect, it was thought that a Laue 
pattern obtained from this same crystal oscillat- 
ing in the manner described above would be of 
interest. 

The crystal used was an X-cut crystal of 
dimensions 100X6X2 millimeters. The long 
dimension was in the Y direction. A detailed ac- 
count of how the crystal was cut is given in a 
previous paper.’ The Laue patterns were made by 
passing x-radiation obtained from a Mo target 
excited by 30 kv through the crystal in the 2 
millimeter direction. The photographic films 
were placed at a distance of approximately two 
centimeters from the crystal. One face of the 
crystal was highy polished and the other face 
was etched with hydrofluoric acid. 

The pattern obtained with the etched face of 
the non-oscillating crystal next to the film is 
shown in Fig. 1A. It is readily seen that the spots 
close to the origin are double in nature and that 
the inner component of most of the spots thus 
situated is weaker in intensity than the outer 
component. This result would indicate that the 
reflecting power of the surface closest to the film 
was weaker than the reflecting power of the sur- 
face closest to the source of x-rays; this result was 
to be expected. The spots farther out from the 
origin show considerably less difference in the 
intensities of their components. This result may 
be explained by the fact that the radiation re- 
flected from the surface of greater reflecting 
power has to pass through more of the crystal in 
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Fic. 1. A, pattern with etched face of non-oscillating 
crystal next to film. B, same with etched face of crystal 
next to x-ray source. C, pattern with crystal in position of 
A, but oscillating. 


reaching the photographic film than does the 
radiation which is reflected from the face of 
weaker reflecting power. It will be observed that 
the inner ends of the remote spots are stronger 
than the outer ends. Here, the radiation reflected 
from the surface of greater reflecting power has to 
traverse so much more of the crystal in reaching 
the photographic film than does the radiation re- 
flected from the surface of weaker reflecting 
power that the radiation reaching the film re- 
flected from the former face is actually weaker in 
intensity than the radiation reaching the film re- 
flected from the latter face. 


The pattern obtained with the etched face of 
the non-oscillating crystal next to the source of 
x-rays is shown in Fig. 1B. It is readily seen that 
the relative intensities of the components of the 
spots close to the origin are in the reverse order in 
which they occur in Fig. 1A. This fact shows that 
the differences in intensities of the components of 
spots close to the origin is due to differences in the 
reflecting powers of the two surfaces of the 
crystal. Thus the result of J. M. Cork might be 
explained by supposing that, although both sides 
of his crystal had been etched, the reflecting 
power of the face nearest the film was stronger 
than the reflecting power of the face nearest the 
source of x-rays, and that the reflecting power of 
both faces was still greater than the reflecting 
power of the inner layers of the crystal. His result 
might be explained in this manner without 
elaboration of the usual Laue theory. 

The pattern shown in Fig. 1C was obtained 
when the crystal was in the same position that it 
was when the pattern shown in Fig. 1A was 
taken; however, in the case of Fig. 1C, the crystal 
was oscillating piezoelectrically in a direction 
parallel to the long dimension under the in- 
fluence of the transverse effect. Very little differ- 
ence, if any, can be observed in the two patterns. 
When the crystal is oscillating in this fashion, a 
loop of pressure exists through the crystal in the 
center of the long dimension. The x-rays were 
passed through the crystal at this point. It might 
be supposed, then, that the reflecting power of 
the etched face should be increased more than 
the reflecting power of the polished face, and that 
the difference between the two components of the 
spots close to the origin should be decreased by 
oscillations of this fashion. It might be further 
expected that the reflecting power of the inner 
layers of the crystal would be uniformly in- 
creased and consequently the blackening between 
the two components of the double spots would be 
relatively increased by oscillations of this fashion. 
This point requires further investigation. 
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On the Interpretation of Heat in Relativistic Thermodynamics 


R. C. Totman, California Institute of Technology ANv H. P. ROBERTSON, Princeton University 
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This article investigates the interpretation of the right- 
hand side of the relativistic second law of thermodynamics 
dx* d 
(o—) (—p)ide'detavtat a 
ds ” To 
and shows that the quantity dQ» can be interpreted as the 
heat—measured by a local observer at rest in the fluid 


at the point of interest—which flows relative to the fluid 
into an element of the fluid having the instantaneous 
proper volume d Vo during the proper time dfo, where these 
quantities are chosen so as to satisfy the numerical equality 
d Vodty = (—g)\dx'dx2dx'dx' and the quantity T> is taken as 
the temperature ascribed to this heat by the local observer, 





§1. INTRODUCTION 


HE analogue in relativistic thermodynamics 

of the usual second law of thermodynamics 

may be conveniently written in either of the two 
equivalent forms! 


dx* dQ 
(+ —) (—g)'dx'dx?dx'dx*t =— (1) 
ds Jy, To 
or 


fs] dx* dQo 
—— (+ — (-—g) \)asatastata—, (2) 
Ox" ds 


0 


where the sign “‘is greater than’’ (>) applies to 
irreversible processes and the sign ‘‘is equal to” 
(=) applies to reversible processes. 

The quantity ¢ occurring in these expressions 
is the proper density of entropy as measured by a 
local observer at rest at the point of interest 
with respect to the thermodynamic fluid or work- 
ing substance which is under consideration. The 
quantities dx*/ds are the components of the 
macroscopic “‘velocity”’ of this fluid at the point 
of interest. And the quantity dQ)/T>) may be 
described as the heat which flows into the element 
of the fluid and during the time denoted by 
dx'dx*dx*dx‘, divided by its temperature—both 
of these quantities being measured in proper 
coordinates. 

In the applications of relativistic thermody- 
namics which have so far been made, there has 


1 Tolman, Proc. Nat. Acad. 14, 268 (1928); ibid. 14, 
701 (1928); Phys. Rev. 35, 896 (1930). 


been no flow of heat relative to the fluid under 
consideration, so that the term dQ)/T7> has 
actually been zero. These applications have been 
in the first place to cases of static thermodynamic 
equilibrium? where the heat flow was zero on ac- 
count of the condition of static equilibrium. In 
the second place, applications have been made to 
the reversible expansion and contraction of cer- 
tain models of the universe* in which there 
has been no flow of heat from one portion of 
the material filling the model to another owing 
to uniformity of conditions throughout the 
models considered. In the third place, applica- 
tions have also been made to the irreversible 
expansion and contraction of models of the uni- 
verse‘ in which, however, there was still no flow 
of heat from one portion of the material to 
another,—again on account of the assumed 
homogeneity of the models. 

In general, nevertheless, the expression of the 
relativistic second law given by (1) and (2) 
should be valid when the heat flow does not hap- 
pen to be equal to zero. And the purpose of the 
present note is to make as clear as may be a cor- 
rect method in the general case for the specifica- 
tion of the quantity dQo/T». 


2 Tolman, Proc. Nat. Acad. 14, 353 (1928); ibid. 17, 
153 (1931); Phys. Rev. 35, 904 (1930); Tolman and 
Ehrenfest, Phys. Rev. 36, 1791 (1930). 

Tolman, Phys. Rev. 37, 1639 (1931); ibid. 38, 797 
(1931); ibid. 38, 1758 (1931). 

4 Tolman, Phys. Rev. 39, 320 (1932); Tolman and Ward, 
Phys. Rev. 39, 835 (1932). 
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§2. INVARIANCE OF dQ)/T») WITH RESPECT TO 
CooRDINATE TRANSFORMATIONS 


Examining the expression for the relativistic 
second law in the form given by (1), it will be 
noted that the two factors (¢ dx*/ds), and 
(—g)'dx'dx*dx'dx*, which are multiplied _to- 
gether to give the left-hand side of the expression, 
are both of them invariant scalars. The first of 
these factors (¢ dx*/ds), is a scalar since it is the 
contracted covariant derivative of a vector and it 
will have a numerical value which is the same in 
all coordinate systems. The second of the factors 
(—g)'dx'dx*dx*dx* is the known expression in 
natural measure for the four-dimensional volume 
specified by the range dx'dx*dx*dx‘, and hence is 
also a scalar, having a value proportional to the 
designated range, but otherwise independent of 
the coordinate system. 

As a result of this invariant scalar character of 
the left-hand side of (1), it is evident from the 
principle of covariance that the right-hand side 
of the expression must also be a similar invariant 
scalar, since otherwise the postulated law would 
not lead to the same results when applied in 
different coordinate systems. Hence the quantity 
dQo/T> must itself have a value, which is of 
course proportional to the four-dimensional 
volume (—g)!dx'dx*dx°dx‘, but is otherwise inde- 
pendent of the coordinate system. This result is 
very useful, since it permits us to obtain a general 
determination of the quantity dQo/7> by first 
employing a specially chosen coordinate system 
in which the interpretation will be quite simple 
and clear. 


$3. DETERMINATION OF dQ)/T) BY UsING 
‘““Co-MoviING”’ COORDINATES 


To obtain a coordinate system which will 
make the interpretation easy, we may take for 
the space-like coordinates x!, x*, x* a three- 
dimensional network of lines permanently con- 
necting adjacent macroscopically identifiable 
points of the fluid under consideration and mov- 
ing with the fluid, and take for the time-like co- 
ordinate x‘ the readings of a set of natural clocks 
which have been distributed to a sufficient 
number of different points throughout the fluid 
and are then allowed to move therewith. The 
possibility of obtaining such a coordinate system 
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would appear guaranteed by the specifications 
for setting it up, and such systems which might 
be called ‘co-moving coordinates” are often very 
useful. 

With such coordinates, it is evident that the 
fluid will always be permanently everywhere at 
rest with respect to the space-like coordinates, so 
that we can write 


dx'/ds =dx?/ds =dx*/ds =0 (3) 


for the spatial components of its macroscopic 
“velocity.”’” Furthermore, for the temporal com- 
ponent we can write 


dx*/ds=1 (4) 


since increments in coordinate time and proper 
time have been made the same by our specifica- 
tion of x* with the help of natural clocks moving 
with the fluid. 

Substituting (3) and (4) in the relativistic 
second law as given in the form (2), we then find 
that this reduces to the simple expression 


fs] dQo 
— (¢0(—g)*)dx'dx*dx*dx‘=— (5) 
ox* To 


and because of the mutual independence of the 
coordinates this can be rewritten in the form 


7] dQo 
_ (d0(— g) tdx'dx*dx*)dxt=—. (6) 


x 0 


To see more clearly, moreover, the significance of 
this result we may now introduce the evidently 
valid substitutions 


dty=dx*; dVo=(—g)*dx'dx*dx', (7) 


where dt) is an element of proper time as meas- 
ured by a local observer at rest in the fluid at the 
point of interest, and dVp is the proper spatial 
volume, associated with the four-dimensional 
region dx'dx*dx*dx*, as measured by this same 
local observer. Introducing (7) in (6), we can 
then write 


a] dQo 
— ($0d Vo) dt5=—. (8) 

Oto To 
The interpretation of dQ)/T> with the help of 
this result is then quite simple. Since ¢o is the 
density of entropy as measured by a local ob- 
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server at rest in the fluid, and dV) is the spatial 
volume as measured by this observer of a small 
region through the boundaries of which no fluid 
is passing, the left-hand side of the expression 
gives the change which the local observer finds in 
the time dé in the entropy of a definite small ele- 
ment of the fluid. In accordance with the principle 
of equivalence, however, this change which a 
local observer finds in the entropy of a system 
small enough so that gravitational curvature can 
be neglected must be related to the heat that 
flows into the system in the way given by the 
ordinary second law of thermodynamics. Hence 
we can now conclude that dQ» is the heat as 
measured by the local observer that flows in the 
time dé) into the element of fluid instantaneously 
contained in the proper volume dVo, and Ty is 
the temperature associated with this heat. Fur- 
thermore, as a result of our previous discussion 
in §2 of the invariance of dQ)/7> to coordinate 
transformations, we are thus furnished with a 
correct value to use for this quantity in any co- 
ordinate system, provided of course that we allow 
for the proportionality between the magnitude of 
dQo/T> and the magnitude (—g)'dx'dx*dx*dx* of 
the four-dimensional volume associated with the 
specified coordinate range. 

Hence, returning to our general expression for 
the relativistic second law, which is valid in any 
coordinate system 


dx* dQ 
(6 —) (—g)*dx'dx*dx*dx4=— (9) 
S7 ps 


Ty 


we may now state that the quantity dQ) occurring 
in this expression can be taken as the heat— 
measured by a local observer at rest in the fluid 
at the point of interest—which flows into an 
element of the fluid having the instantaneous 
proper volume dV during the proper time dh, 
where these quantities are chosen so as to satisfy 
the numerical equality 


dV pdt = (— g) \dx'dx*dx*dx4 (10) 


and the quantity 7) is taken as the temperature 
ascribed to this heat by the local observer. We 
are thus provided with a perfectly definite inter- 
pretation of all the quantities occurring in the 
relativistic second law which will satisfy the 
principle of covariance by leading to the same 
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results no matter what coordinate system is used, 
and which satisfies the principle of equivalence by 
reducing, as we have seen with the help of a 
specially convenient coordinate system, to the 


ordinary second law for an infinitesimal portion 


of the fluid. 


$4. DETERMINATION OF dQo/To BY UsING 
NATURAL COORDINATES 


In accordance with the above specifications 
for determining a correct value to use in the rela- 
tivistic second law for the flow of heat dQ, into 
the element and during the time denoted by 
dx'dx*dx*dx', it will be noted that this quantity is 
to be obtained from measurements by a local 
observer of the rate at which heat flows into a 
specified element of the fluid. Hence this quantity 
is to be regarded as corresponding to a flow of 
heat measured relative to boundaries at rest in 
the fluid, rather than in some way relative to 
boundaries at rest with respect to the spatial 
coordinates that are actually being employed, as 
might at first sight have seemed plausible. Since 
this conclusion was obtained, however, with the 
help of a coordinate system specially chosen so 
that boundaries at rest with respect to the fluid 
were also at rest with respect to the spatial co- 
ordinates, it will be illuminating to show that 
we should also be led to the same interpretation 
for dQo/T» using a system of coordinates with re- 
spect to which the fluid is not at rest. 

To obtain such a demonstration, we shall now 
use a set of coordinates x, y, 2, t with respect to 
which the fluid is not at rest, but which are so 
chosen as to be natural coordinates for the point 
of interest. In accordance with the principle of 
equivalence such natural coordinates for any de- 
sired spacetime point can always be found, and 
with respect to them the principles of special 
relativity as expressed in their usual form may be 
taken as valid in the immediate neighborhood of 
the point in question. Hence, by choosing these 
coordinates, we shall be able to employ the 
principles of thermodynamics as developed for 
the special theory of relativity by Planck and 
Einstein in arriving at the desired interpretation. 

Making use of these natural coordinates x, y, 
z, t, we shall evidently have 


(—g)'=1 
and ' 
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0 a) fe] 0 
—(—g)'=—(—g)'=—(—g)!=—(—g)!=0 (11) 
Ox oy Oz ot 


at the point of interest. Hence by substituting 


) dx 7) dy 0 
|— (60) +— (#2) +— 
Ox ds} dy ds} @z 





dz 0 
do —) +— 


( 
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into the relativistic second law of thermodynam- 
ics in the form given by (2), we can now write the 
left-hand side of that expression in our present 
coordinates in the form 


dt 
(+ —) fesayae 
ds 


(12) 


ds ol 


where dx/ds, dy/ds, and dz/ds are the components of the macroscopic velocity of the fluid at the 
point of interest with respect to proper time, and by a simple substitution we can rewrite this in the 


form 
dt ~*) fs) 


dt “*) 0 
” ds di} av 


te) 

[— (+ (o moe 

Ox ds dt 
where dx/dt, dy/dt and dz/dt give the components 
of the velocity of the fluid as ordinarily expressed 
in terms of the coordinate time ¢. 

In accordance with the special theory of rela- 
tivity, however, entropy is an invariant for the 
Lorentz transformation, and hence entropy den- 
sity will be affected by the Lorentz factor of con- 
traction ds/dt in such a way that we can substi- 


tute 
$= do(dt/ds), 


where ¢ is the density of entropy at the point of 


dg 
oonnel 

at 
Finally, noting the significance of the various 
terms in (16), we may now rewrite this expression 


for the left-hand side of the relativistic second 
law in the simple form 


(14) 





0p d¢ 0p 
u—+v—+w—+¢ 
Ox oy 02 


(do/dt)d Vdt+-$(d/dt)(dV)dt, (17) 
where we have written dV for the volume of the 
element of fluid, instantaneously contained in the 
region dxdydz, and have written d¢/dt for the 
total rate of change in entropy density as we 
follow the moving fluid. Hence in natural co- 
ordinates the left-hand side of the relativistic 
second law is seen to be equal to the change that 
takes place in time dt in the entropy of the ele- 
ment of fluid instantaneously contained in the 
region dxdydz, and we can write for these co- 
ordinates 


02 


dt dz ra) dt 
(+ = -) Pie (« —) feed, (13) 
ds dt al ds 





interest taken with respect to our present system 
of coordinates x, y, 2, ¢. 

Substituting (14) into (13), and writing for 
simplicity u, v, and w for the three components of 
velocity dx/dt, dy/dt and dz/dt, we then obtain 


] re] fe) 0d 

[— (pu) +— (ov) +— (ow) +— |dxdydedi (15) 
Ox oy 02 ot 

and by performing the indicated differentiations 

and rearranging the order, this can be rewritten 

in the form 


Ou ov 


ow 
+= +—) Jisayasar (16) 
Ox Oy dz 





dx* d 
(66 —) (—g)*dx'dx*dx'dxt=— (odV)dt. (18) 
ds Jy dt 


In accordance, however, with special relativ- 
istic thermodynamic theory we can relate this 
change in the entropy, of the little thermody- 
namic system consisting of this element of fluid, 
to heat flow and temperature by the expression 


d/dt (¢dV)di=dQ/T, (19) 


where dQ is the heat absorbed by the element at 
temperature 7 and in the time dé referred to our 
present coordinates. Furthermore, since the ratio 
of heat to temperature is an invariant for the 
Lorentz transformation, we can also take 


dQ/T =dQo/To, 


where dQ» is the absorbed heat as measured in 


(20) 
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proper coordinates and 7) is the proper tempera- 
ture of the element. Moreover, in accordance 
with the Lorentz contraction for volume elements 
and Lorentz dilation for time intervals we can 
write 


dVdt=dV dbo, (21) 


where d Vp is the volume of this element as meas- 
ured in proper coordinates and df) is the proper 
time during which the heat absorption takes 
place. 

Hence combining (18), (19) and (20), we have 
now obtained, also using our present coordinates, 
the result 


(22) 


0 


dx* dQy 
(+ —) (—g)8dx'da2dx8dxt=—, 
ds, T 


where, in accordance with (21), dQ) is to be taken 
as the heat—measured by a local observer at 
rest in the fluid at the point of interest—which 
flows into an element of the fluid having the 
instantaneous proper volume dV» during the 
proper time df, these quantities being chosen so 
as to satisfy the numerical equality 


dV dtp =d Vdt = (—g)'dx'dx’dx'dx*. (23) 


Thus, using natural coordinates, we have satis- 
factorily obtained the same interpretation of the 
quantities occurring in the relativistic second law 
as was obtained in the preceding section using 
co-moving coordinates. And the treatment gives 
the desired illustration of the fact that our inter- 
pretation of dQ) as heat flowing through a 
boundary stationary in the fluid was not a result 
of our original special choice of coordinates in 
which the fluid was at rest. 


§5. CONCLUSION 


Before concluding, two further remarks may 
be made concerning the interpretation of the 
relativistic second law, which may be illumin- 
ating. 

In accordance with the treatment given in the 
last section it is to be noted that the quantities 
¢, ou, dv and ¢w must be regarded as fofal densi- 
ties of entropy and entropy flux at the point of 
interest, since otherwise, for example, the left- 
hand side of expression (19) could not have been 
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interpreted as the change in the total entropy of 
the specified element of fluid and correlated as 
was done with the influx of heat. This implies 
that we are to treat our thermodynamic fluid 
from a macroscopic point of view and take ¢ as 
the total entropy associated and moving along 
with unit volume of that material, rather than to 
try to employ a microscopic point of view and 
take ¢ as entropy belonging in some way solely 
to the molecules of the fluid with an additional 
quantity belonging to the radiation in the space 
between them. Similar remarks apply to ¢» and 
godx*/ds. It will be noted that the correct pro- 
cedure is in entire agreement with ordinary 
thermodynamic practice and with the spirit of 
thermodynamics as a macroscopic phenomeno- 
logical science. 

In laying down a correct method for deter- 
mining the magnitude of the heat flux divided by 
temperature dQo)/7T>) which corresponds to the 
infinitesimal element of fluid and time denoted 
by dx'dx*dx*dx*, it will have been noted that we 
have merely specified that the local observer 
should measure the proper heat dQ) entering an 
element of the fluid of instantaneous proper 
volume dV, during the proper time dt, where 
dV» and dt& are chosen so as to satisfy the 
numerical equality 


dV odty = (—g) 'dx'dx2dx3dx' (21) 


and have not introduced any specifications as to 
the shape of the four-dimensional region d Vodh. 
This, however, will be seen to be legitimate since 
the heat flux is directly proportional to the prod- 
uct of volume and time, independent of shape 
except for higher order differentials which can be 
neglected on shrinking the four-dimensional re- 
gion down to the point of interest. 

It is hoped that the treatment which we have 
given in this article will help to make the interpre- 
tation of the relativistic second law clear. The 
considerations presented are in agreement with 
the results obtained in previous applications of 
relativistic thermodynamics where there was no 
flow of heat relative to the fluid under consider- 
ation, and should be specially helpful when 
further applications are made to systems in 
which a flow of heat relative to the fluid does 
take place. 
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The adiabatic tension coefficient of resistance of single 
antimony crystals has been determined for various orienta- 
tions. Since antimony crystals have the same type of 
symmetry as bismuth, the curves connecting the coefficient 
with the primary and secondary orientations have the 
same general shape as in the case of bismuth and again 


‘uphold the theory put forth by Bridgman. The observed 


values of the coefficients for antimony when the tension 
and current are parallel to the trigonal axis and per- 


pendicular to it are different both in magnitude and in 
sign, whereas in the case of bismuth the two were both 
negative and very nearly equal in magnitude. The six 
coefficients necessary completely to determine the behavior 
of the resistance when deforming forces are applied to the 
antimony crystals have been found (1) without correcting 
for strain and (2) correcting for the changes arising from 
the strain produced by the tension. 





N a previous paper' the author has studied 

experimentally the effect of tension on the 
electrical resistance of single bismuth crystals. 
The experimental dependence of the tension 
coefficient on the principal and secondary 
orientations of the crystals with respect to the 
cylindrical axis was later found by Professor 
Bridgman? to be in agreement with a geometrical 
theory based on the symmetry of the crystal. The 
object of the present paper is to examine the 
effect of tension on the resistance of single 
crystals of antimony which has the same type of 
symmetry as bismuth and to see whether the 
results obtained are consistent with the theory. 
The values of the six tension coefficients are 
found numerically, both as observed directly and 
as corrected for the change in resistance due to 
the strains in the crystal produced by the 
tension. 


PROCEDURE 


Experimentally the procedure is somewhat 
more difficult than in the case of bismuth. The 
specific resistance of antimony is roughly one 
third that of bismuth so that the resistance of a 
similar piece of antimony is approximately one 
third as great. Furthermore, since the tension 
coefficient of resistance of polycrystalline anti- 
mony is of the order of one tenth’ that of 


! Mildred Allen, Phys. Rev. 42, 848 (1932). 
?P. W. Bridgman, Phys Rev. 42, 858 (1932). 
*P. W. Bridgman, Proc. Amer. Acad. 57, 41 (1922). 


polycrystalline bismuth, it is to be assumed that 
the tension coefficients of antimony will be of the 
order of one tenth of those of bismuth. Conse- 
quently, the absolute change in resistance to be 
measured will be about one thirtieth that studied 
in the case of bismuth. This is counteracted in 
some degree by the fact that twice and even four 
times as great a tension (force per unit area) as 
had been necessary to use with bismuth was here 
safely applied. On the other hand, antimony is 
more brittle and so is more difficult to handle. 
Therefore the apparatus had to be studied 
afresh to see in what way the sensitivity could be 
increased—and preferably to about ten times its 
previous value. Considerable increase in sensi- 
tivity was obtained by using twice as large, and 
sometimes four times as large, a current as in the 
work with bismuth. The most obvious way 
further to increase the precision of the method 
was to use a more sensitive galvanometer, but 
even the new Leeds and Northrup galvanometer 
was rated to have only twice the voltage 
sensitivity of the older galvanometer already in 
use. It was pointed out in the previous paper! 
that when the deflection method was used the 
voltage sensitivity of the galvanometer depended 
on the resistance r of the compensating circuit Q 
across which the galvanometer was connected 
and was in fact proportional to the sum of this 
resistance r and the resistance of the galvanome- 
ter itself. It is thus of advantage to use a small r. 
Moreover, in comparing small changes of re- 
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sistance it seemed desirable to attempt to use the 
galvanometer so that its voltage sensitivity 
remained constant. This was effected by choosing 
r as 3 ohms plus half the length of the slide 
wire in the auxiliary circuit, and by roughly 
balancing the potential fall across the crystal 
against the auxiliary potential difference by 
varying the large resistance governing the cur- 
rent in the auxiliary circuit. The final precise 
balancing could then be carried out by moving 
the contact on the slide wire a small distance. 
For this reason the large resistance in the auxili- 
ary circuit Q was chosen to be a 10,000 ohm 
variable radio resistance. The voltage sensitivity 
of the galvanometer, with the scale at 5 meters as 
before, was 2.15107 volts per cm. Readings 
then were satisfactorily large as they varied from 
0.15 cm to 4.00 cm. With this deflection method it 
proved to be possible to measure the tension 
coefficients to within about 110-* in absolute 
value where the coefficient itself varied from 
+9.0X10-* to —35.010-®. The null method 
used before was discarded, since the changes of 
resistance were too small to give more than a 
millimeter’s displacement on the balancing slide 
wire of the main circuit; its only use in this set of 
experiments was to check the sign of the change 
of resistance. 

A further precaution was to decrease the size of 
the air chamber containing the crystal and the 
surrounding copper cylinder, so that temperature 
equilibrium was obtained more rapidly. The 
substitution of a mercury thermostat for the 
metal one and of a Burgess vacuum switch for the 
tilting mercury switch made the temperature less 
variable and entirely independent of the room 
temperature. This greater stability of tempera- 
ture was of importance in measuring the 
markedly smaller effect and was reflected in the 
character of the deflection runs made with 
alternate readings with and without tension, 
Fig. 1, which in this work with antimony usually 
showed a definite periodic variation. This peri- 
odicity must arise primarily from the periodic 
changes in temperature produced by the regular 
increases and decreases in heating current gov- 
erned by the thermostat. In addition there is 
usually to be noted a small gradual displacement 
of the whole curve up or down, and this drift in a 
given direction can reasonably be ascribed to the 


Deflection 





Time 


Fic. 1. a without tension; } with tension. 


zero drift of the galvanometer and to the unequal 
changes of voltage in the two balancing circuits 
arising very probably from the polarization of 
the dry cell of the auxiliary circuit Q. As is seen 
in Fig. 1 this drift is very small when temperature 
equilibrium has been thoroughly established 
before beginning to take readings. 

The crystals used were of Kahlbaum antimony. 
Some of them had been made and used by 
Professor Bridgman and were generously put by 
him at the writer’s disposal; they had diameters 
of approximately 1/8’. Others were crystallized 
by the writer, following the procedure indicated 
by P. W. Bridgman‘; most of these latter had a 
diameter of approximately 1/16’. The advantage 
of the smaller diameter was the greater resistance 
giving a correspondingly greater change of 
resistance; the disadvantage, the greater difficulty 
in handling this brittle metal. 


RESULTS 


I. Without correcting for strain 

Twenty-one separate individual single crystals 
of antimony have been studied and their tension 
coefficients measured. Because of the smallness of 
the effect it proved not to be sufficient to 
measure the coefficients once with the current 
flowing in but one direction. With most of the 
crystals the coefficient was determined with one 
tension and then under identical conditions with 
the current reversed; to check this result this 
entire process was then repeated with either the 
magnitude of the current or of the tension 
changed. The two determinations were likely not 
to agree if the crystal were slightly cracked. 
Furthermore the accuracy with which the coeff- 


4P. W. Bridgman, Proc. Amer. Acad. 63, 351 (1929). 
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cient could be determined under given conditions 
without reversal of the current was found by 
taking two determinations with only the galva- 
nometer leads reversed in the two cases. In only 
a few cases did the difference between these two 
readings differ by more than 0.3 X 10~* in absolute 
value. On the other hand, the readings with the 
current reversed were likely to differ by a 
considerably larger amount, usually of the order 
of 1.0X10-*. These differences were thus larger 
than the error to be expected as intimated with 
the galvanometer leads reversed. It seemed at 
one time that these discrepancies were pro- 
portional to the magnitude of the current used, 
and so were to be ascribed to a tension coefficient 
of the Peltier effect, but further work showed this 
explanation to be untenable, so that probably 
they arise from some irregular heating effect and 
are eliminated by taking the average of the 
readings with the direct and reversed currents. 

The values of the tension coefficient were then 
to be studied on the basis of Professor Bridgman’s 
theory” as functions of the primary orientation of 
the crystal 6 and of the secondary orientation ¢. 
The angle @ was defined as the angle between the 
trigonal axis of the crystal and the cylindrical 
axis. The angle ¢ was defined as the angle 
between the projection of the cylindrical axis on 
the principal cleavage plane (which is perpen- 
dicular to the trigonal axis) and the projection on 
the same plane of the normal to that secondary 
cleavage plane which makes an angle of about 71° 
with the principal cleavage plane. The relation 
derived by Professor Bridgman was the following: 


B=(1/pe) { p11 sin* 6+ p33 cos‘ 0 
+(2pis+ p44) sin? 6 cos* 6 


—2pi,sin* @cos @cos3¢}, (1) 


where 8 is the tension coefficient, which is the 
relative change in resistance AR/R produced by a 
tension 7 of 1 kg/cm’, the p,,.’s are five of the six 


5 Since pi; and p33 determine the tension coefficients for 
6=90° and @=0°, respectively, it appears from Eq. (2) 
that the tension coefficient for @=90° is known with less 
error than for @=0°. The coefficient p;4 gives the distance 
between the curves for ¢=0 and for ¢=60° for a given 
primary orientation; this ‘‘spread’’ of the curves is very 
accurately known. The remaining constant (213+ 4.) 
determines the exact shape of the curve for ¢=30° and 
the error here is seen to be comparatively large. 


tension coefficients necessary completely to define 
the behavior of the resistance of a trigonal 
crystal under the action of deforming forces, and 
pe is the specificresistance at theorientation 6. The 
expression is the same as that used for bismuth, 
since bismuth and antimony crystals have the 
same type of symmetry. The term in ¢ shows the 
trigonal symmetry. Since none of the twenty-one 
readings of the coefficient probably were entirely 
consistent with any of the others and the tension 
coefficients were quite small, it seemed worth 
while to use least squares in evaluating them. 
The values of the four constants and their 
probable errors’ as found in this way directly 
from the experimental values are 


pu=(+ 0.24+0.09) x 10-™, 
ps3 = (—12.20+0.33) X 10-™, 


2 pint pas =(4+5.72+0.53) KX 107", 
Pis=(+5.10+0.16) K10-". 


(2) 


In Fig. 2 the full lines give the curves computed 
from these values for ¢= 60°, ¢= 30°, and y= 0°. 


+12 
+8 
+4 
0 
-4 
-8 
eo 
So -l2 
a” 
“loc -16 
a 


-20 
-24 
-28 
-32 
-36 


-40 





0 410 20 30 40 50 60 70 80 90 
8 
Fic. 2. Dependence of tension coefficients on orientation. 
a, ¢=0°; b, e=30°; c, ¢e=60°. Full curves, experimental 
curves. Dotted curves, corrected for strain. 
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This single set of curves is sufficient entirely to 
describe the behavior of the electrical resistance 
of antimony under tension, since from Eq. (1) 
it is known that the coefficient always varies 
between these indicated outside limits pro- 
portionately to cos 3y. These curves are like the 
similar ones for bismuth in that the coefficient 
has a single constant value for 6=90° and for 
6=0° (as it must according to the theory) and 
varies in sign according to the orientation; but 
they differ in that for bismuth the values of these 
constants have the same order of magnitude and 
the same sign, whereas for antimony the absolute 
magnitude at 6=0° is more than fifty times that 
for 2=90° and the signs of the two are different. 
A numerical comparison of the observed experi- 
mental values and the computed values of the 
tension coefficient with the constants of (2) is 
given in Table I and it is seen that the difference 
between the two only once exceeds 1.110~, 











TABLE I. 

6 @ R X 10%) expt. theor. | ©XP*: _ 
and eeemmeieiedl corrected for 
size strain 

90° 16° BL 2.09 |/+ 1.03 + 0.54/— 1.92 — 2.40 
83° 46° BL 2.02 |+ 3.04 + 2.74/4+ 0.39 + 0.11 
83° 24° As 1.98 |+ 0.96 — 0.14/— 2.08 — 3.16 

79° 17° As 2.48 |— 2.03 — 1.70})— 5.25 — 4.92 

77° 55° BL 1.96 |+ 5.05 + 5.68/+ 2.83 + 3.49 

75° 40° BL 2.18 |+ 4.02 + 3.87/+ 1.58 + 1.44 

75° 36° As 5.74 |+ 1.79 + 2.83|— 0.78 + 0.26 
73° 03° BL 1.99 |— 4.52 — 4.61}/— 8.13 — 8.20 
72° 06° AL 5.51 |— 443 — 4.57;/— 8.00 — 8.16 
67° 51° BL 3.10 |+ 6.76 + 7.88|+ 4.97 + 6.09 
65° 38° BL 5.96 |+ 3.99 + 3.24/+ 1.74 + 0.83 
64° 34° As 4.32 |+ 2.644 + 3.06/+ 0.21 + 0.65 
61° 17° BL 1.50 |— 3.46 — 3.75|— 6.75 — 7.07 
62° 07° BL 1.16 |— 6.39 — 5.99|—10.02 — 9.64 
62° 04° BL 2.01 |— 5.644 — 6.34/— 9.32 —10.03 
68° 01° As 2.18 |— 6.33 — 5.59|— 9.99 — 9.27 
54° 53° As 1.93 |+ 7.25 + 7.06/+ 5.62 + 5.45 
50° 51° AL 4.37 |+ 5.04 + 4.98/+ 3.29 + 3.22 
53° 36° BL 1.41 |+ 3.47 + 1.92/+ 1.16 — 0.38 
44° 06° As 1.64 |—11.5 —1046)—15.2 —14.1 
22° 56° As 0.86 |—21.7 —21.77/—25.6 —25.6 














B indicates crystal was made by P. W. Bridgman. 

A indicates crystal was made by the author. 

L designates a large crystal of approximately 1/8’’ diam- 
eter. 

s designates a small crystal of approximately 1/16’’ diam- 
eter. 
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which considering the smallness of the effect 
indicates satisfactory agreement. The one reading 
which differs by more than this amount lies on a 
very steep part of the curve where a small error in 
measuring ¢g necessarily produced a large error in 
B. 
Only two readings were made for @ less than 
50°. This came about for two reasons: it is 
comparatively difficult to grow crystals with the 
principal cleavage plane nearly perpendicular to 
the cylindrical length, and, more important, 
crystals with such orientations are so exceedingly 
brittle that it is very nearly impossible to mount 
them in the apparatus without sufficient jarring 
to break them along a principal cleavage. The 
point for @= 22° fitted within experimental error 
on the curve computed from the points for 6 
between 44° and 90°, the value of 8 computed on 
this basis being —22.910~-* as compared with 
the experimental value of —21.710-*. This 
furnishes a check on the correctness of the general 
trend of the curve which seems entirely satis- 
factory, although in the final calculation this 
lowest point was included and the values of the 
constants modified somewhat thereby. 
Completely to determine the six coefficients 
which are necessary to describe the behavior of 
a trigonal crystal under the action of forces, the 
data derived from the effect of tension on anti- 
mony crystals must be supplemented by the 
values of the hydrostatic pressure coefficients 
when the current is flowing parallel to the 
trigonal axis and when it is perpendicular to it. 
Professor Bridgman’s best values for these 
pressure coefficients at 30°C are: 


[(1/p)(AR/R) Joo? = +3.7 X10, 
[(1/p)(AR/R) Joo = +19.4X 10-. 


6=90° 
6=0° 


These are taken from a recent paper‘ in which 
the pressure coefficients are carefully studied as a 
function of pressure; these values are extrapolated 
a short distance to zero pressure to be comparable 
with the very small forces used in the tension 
experiments. We have therefore the additional 
equations as given by Bridgman? 


+3.7 X10 = —(piitpi2tpis)/(penso?), (4) 
+19.4X 10% =—(2pi3+p33)/(pomc)- (5) 
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This gives finally for the six tension coefficients of 
resistance uncorrected for the changes of re- 
sistance arising from the change in dimensions 
caused by the application of the forces: 
pu= +0.2 x 10-", 
net 12.2 x 10-", 


pia= +5.1 x 10-", 


pi2= — 4.6 X10-", 
Pi3> +2.7X10-", 
pas= +0.2X10-”. 


(6) 


It is to be noted that these are of the order of 
magnitude of one tenth those found in the case of 
bismuth as was anticipated from the poly- 
crystalline tension coefficients. This polycrystal- 
line coefficient was found by Professor Bridgman 
to be approximately 5 <10~* and this lies within 
the limiting values found in this set of experi- 
ments. One would be inclined to explain his one 
reading of 20.0 10-* as arising from a crack in 
the antimony which was not otherwise detect- 
able, particularly as this metal is so very brittle. 


II. Correcting for strain 


In the case of bismuth the changes of resistance 
arising from the deformation of the crystals were 
so small compared with those actually produced 
by the tension that, being within experimental 
error, they could very properly be neglected.® In 
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the case of antimony, however, with the in- 
creased precision of the method, the deformation 
of the crystals gives changes greater than the 
experimental error, sometimes of the order of 
magnitude of the tension coefficients themselves, 
and should be considered, although in practice 
the observed change of resistance will always 
include that arising from the change of dimension. 
As a matter of fact, since these corrections do not 
follow the same law of variation with the 
orientations as do the tension coefficients, Pro- 
fessor Bridgman’s equation is rigorously only 
applicable when these corrections have been 
made. 

Correcting for strain is not an entirely simple 
matter. As Professor Bridgman has pointed out 
in the case of hydrostatic pressure’ the change of 
resistance due to the strain produced by the 
tension will consist of three parts: (1) that arising 
from change in length; (2) that from change in 
cross section; and (3) that from change in the 
orientation 6. These will be considered separately. 

(1) To find the relative change of resistance 
arising from the change in length produced by the 
application of a tension of 1 kg/cm? as a function of 
the primary orientation 0 and of the secondary 
orientation ¢. The elastic behavior of antimony’ 
is defined by the equations 


rz =SuXztSi2Vyt+sisZ2t+s5u¥,+0Z, +0X,, 
Cyy=S12Xet+5i1Vyt+si3Z2—514Y,+0 +0, 

C€22=S13X2t+513Vy+533Z,+0 +0 +0, (7) 
Cyz=SuX2—SuVytO +54,Y,+0 +0, 

éez=O0 +0 +0 +0 +4544Z2 +25uXy, 

€ry=O0 +0 +0 +0 = 4254Z2+2(Si1i—Si2)Xy, 


* This has been shown quantitatively by applying the 
method of least squares to the bismuth data. The numerical 
values of the coefficients are, giving Professor Bridgman’s 
values for comparison, 








Least Squares 





Corrected 
Bridgman Uncorrected for strain 
pu — 7.7X1079 — 7.2107 — 7.7X107° 
P33 — 6.6 — 6.2 — 6.7 
pia +15.7 +13.5 +13.2 
pir + 5.6 + 4.3 + 4.2 
pis + 1.8 + 1.6 + 1.9 
Pua —12.3 —10.4 —13.1 








The uncorrected values differ somewhat from those given 
by Professor Bridgman in his theoretical paper* since his 
values were computed from six graphically interpolated 
points and these are based on all forty-five directly ob- 
served points. The average difference between the experi- 
mental values of the tension coefficients of resistance and 
those computed from these uncorrected constants is 
0.6 X 10-5; in only four cases is this difference greater than 
1.2 10-5, Since the correction of 8 for strain for bismuth 
lies between 0.36 107 and 0.83107, depending on the 
crystal orientation, it is evident that the strain correction 
is of the same order of magnitude as the experimental 
error. 
7P. W. Bridgman, Proc. Amer. Acad. 60, 305 (1925). 
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in which Z represents the direction of the 
trigonal axis. The constants for this case have the 
values, T being expressed in dynes/cm?, 


Su=+17.7X10-",  s33=+33.8X10-%, 
S=— 3.8X10-%, sy=+41.0X10-%, (8) 


Sis=— 8.5X10-%, 8.0X 107%. 


i tlio 


This immediately gives us, for 6=0°, when 7 is 
given in the practical units of kg/cm’, 


(1/T)(Al/l) =533X 980 X 1000 = +3.3110-*, (9) 


which gives the order of magnitude of the 
correction to be expected from this source. 

To find the most general dependence of 
(1/T)(Al/1) on @ and ¢, it will be convenient to 
transform from the coordinates X YZ defined by 
the trigonal axis and the crystalline properties of 
the antimony to X’Y’Z’ axes defined by the 
characteristics of the cylinder in which the metal 
is cast. The relative elongation e.’,’ will then give 
the change in length desired. We must first 
define these two sets of coordinates more pre- 
cisely. The Z-axis, Fig. 3, is to be parallel to the 








Fic. 3. 


trigonal axis of the crystal, the Z’-axis to the 
cylindrical axis of the casting. The angle between 
them is then by definition @. The Y-axis lies in 
the elliptical cross section of the principal 
cleavage plane and is the projection on this plane 
of the normal to the secondary cleavage plane 
making an angle of about 71° with the principal 
cleavage plane; thus the Y-axis makes an angle ¢ 
with the major axis OQ, by definition. The 
Y’-axis cannot of course lie in this plane, but it is 
to be in the plane QOZ’, so that OQ is its pro- 
jection on the principal cleavage plane. The nine 
direction cosines necessary completely to define 
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the transformation are then known and are 
recorded in Table II. In the new set of coordi- 











TABLE II. 
ve of Zz 
P 4 cos ¢ sin ¢ cos 0 sin ¢g sin 6 
Y —sin ¢ cos ¢ cos 6 cos ¢ sin 6 
Z 0 —sin 0 cos @ 








nates the equations of elasticity have the form, 
since the tension is always along the Z’-axis, 


m J} — yi 
Cx's’ =S'13Z ss €y eg ™S a3Z 2’, 


vy , = c¢’..7' 
Cy'’y’ +S e222’,  C2's' =S' 532 2’, 


(10) 


—— a se: 
C2’2'=S'33Z 2’, Czy’ =S 63Z 2, 


where the six s’’s define the behavior of the 
crystal in the primed set of axes and are known 
functions of the s’s, the transformation of the 
s’s being given by Voigt.’ Carrying out the 
calculation for s’3;, we get for the relative 
change in length along the axis of the cylinder 
the equation 


(1/77) (Al, l) =S1) sin 6+ S33 cos‘ 6 
+ (2513+ 544) sin? 6 cos" 0 


— 25,4 sin® 6 cos @ cos3y. 


(11) 


The term in cos3¢ is in agreement with the 
known trigonal symmetry. 

(2) To find the relative change of resistance 
arising from the change in cross section produced 
by the application of a tension of 1 kg/cm’. 
It is to be noted that OY’ is the major axis of the 
horizontal elliptical cross section into which the 
tension deforms the circular cross section. This 
cross section originally had the area rr* and 
when strained has the area rr*(1+-5’13)(1+5'23) so 
that 

(1 T)(AA/A) =8'i3+5'23, (12) 


where it is to be expected that one or both of 
these s’’s will be negative. Evaluation of 4A/A 
in terms of the s’s gives for the correction due to 
change in cross section 
(1/7 )(AA /A) =5)3(cos* 6+sin* 6+ cos? @) 
+5 2 sin? 6+ (53;+533— 544) sin® 6 cos? 8 
+2514 sin*® 6 cos @ cos 3¢. 


(13) 


®W. Voigt, Lehrbuch der Kristallphysik (pp. 589-592), 
1910. 
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The values computed from this formula are 
negative as was to be expected. 

(3) To find the relative change of resistance 
arising from the change in the angle 0 produced by 
the application of a tension of 1 kg/cm?. Since the 
specific resistance of a crystal is given by the 
relation 


Pe= Po cos* 6+ P90 sin? 6 (14) 


we have 
dp 2(ps0— po) 


— =——————- sin @ cos 6d0. 
p pe 


(15) 


It is immediately evident, as inspection shows to 
be true also in the case of the previous two 
corrections studied, that the change in resistance 
produced by the application of tension is inde- 
pendent of the secondary orientation at @=0° 
and at 6= 90°. In the case of this third correction 
dp vanishes at both these limiting points. It 
remains to find the value of d@ produced by the 
tension; d@ signifies a rotation about the X’-axis. 
In determining d@, it must be borne in mind 
that the displacement of a point during the 
deformation of a body is compounded of a pure 
strain and of a rotation, the former being defined 
by the six elastic constants and the latter by 
three additional constants. Let us consider a 
small section in the middle of the cylindrical 
casting of the crystal, and let 0 be the very middle 
point which, since the tension consists of two 
equal and opposite forces, must remain un- 
changed in position. Let Q2 be a point on the 
principal cleavage plane through 0 and in the 
plane Y’OZ’, Fig. 3. Q: is a point on the cylin- 
drical axis having the same Z’ coordinate as Qe 
relative to 0; this may be called ¢. The angle 
?,0Q2 is by definition (90—@). Similarly we may 
call » the relative Y’ coordinate of Qo. 
According to Love, the most general displace- 





1 Ap 
— —=2 sin? 6 cos? 0 
T p 
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ment of a point Q after deformation when 0 is 
kept constant is given by the three expressions: 


Cort t+ dezynt perf —Nwstfw2 =dx’, 
2CryEt Cyynt dye —Swittws=dy’, 
b€r2€ + 3ey.nt+ e::5— Ewetnw =dz', 


(16) 


where £, n, ¢ are the coordinates of Q relative to 0 
and the w’s are the three additional constants 
defining the rotation. Applying these equations 
to the point Q,; when a tension is acting along 
OQ), i.e., along OZ’, we have the conditions, since 
f:=m=0, 


dx;'=0 or 4¢z','+we=0, 
dy;/=0 or 3e,'.'—w,=0, (17) 
dz,’ =e,’,'¢. 


Applying Eq. (16) to the point Q2 gives the 
relations, £2 being zero, 


diva! = hes'y't cot O+}e2"s'S—F cot Ows+iws, 


dye’ = ey'y’f cot 06+-4e,'.°%—fw, (18) 


dze’ = hey. cot O+ e226 +f cot Ow. 


The assumption that there is no motion in the X’ 
direction leads to a perfectly consistent result, 
i.€., 3€:’y’=3. From the last two conditions in 
Eq. (18) we can get the new angle between ‘00: 
and OQsz, since 


5s (+dz’ 
tan (6+d6é) =———— (19) 
nt+dye’ 
so that we have the relation 
d0=sin 0 cos 0{e2'2'—Cy’y’ tey’z’ cot A}. (20) 


This means that the final change in resistance 
due to change in angle is 


Ps0— Po 


{ S33’ — So3’ +543" cot 6} ’ (21) 


and the problem again reduces to the determination of the s’’s involved. On this basis, the final 
expression for the correction of the specific resistance of the crystals arising from the change in the 


angle @ is 


1 Ap 
——=2 sin? @ cos’ 6 
T p pe 





{ (Sir Sis) sin? @— (S33—Si3— S44) Cost 6+ (53:—5S33+544) sin® 6 cos? 0 


—2s,,4sin cos @cos3¢}. (22) 








576 


Thus the total effect on the resistance produced 
by the strain caused by the application of tension 
to trigonal crystals is given by the Eqs. (11), 
(13) and (22). Since R= pl/A, the relative change 
of resistance is 


AR/R=Al/I—AA/A+Ap/p, (23) 


the correction for the change in cross section 
being subtracted. The numerical values of these 
corrections in the case of antimony are plotted 
in Fig. 4 in the same way as the tension coeffi- 


1 aR 6 
T R™ 10 


a( 





0 10 20 30 40 50 60 70 80 90 
8 


Fic. 4. a, ¢=0°; b, e=30°; c, ¢=60°. Strain corrections of 
the tension coefficients. 


cients were plotted in Fig. 2. These data are to be 
used to correct the four constants determined 
from the experiments with tension, the strain 
corrections being subtracted from the experi- 
mental values and least squares being applied 
again to these corrected observations. The 
resulting constants and their probable errors are 


pu=(— 1.07+0.09) x 10-™, 
pss = (—13.90+0.54) X10-™, 
2pistou=(+ 4.34+0.34) x 10-, 
pu=(+ 5.81+0.16) x10-. 


(24) 


The curves in which these corrections have been 
made are indicated in the dotted lines of Fig. 2 
and the corrected values of the tension coeffi- 
cients in the last two columns of Table I. 

The two constants derived from pressure 
measurements must also be corrected for strain. 
The analysis for this correction for pressure has 
already been carried out by Professor Bridgman.’ 
The correction for angle is again zero for the 
limiting cases 6= 0° and @=90° and that arising 
from the combined changes in length and cross 
section is 
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((1/p)(AR/R)) =2(siit+siz) cos? 0 
+2513 sin? 6+533(1—2 cos? 6). (25) 


Evaluating this correction in the two limiting 
cases which are needed, gives the values 


@=90° ((1/p)(AR/R))o0= +5.4X 10-, 
@=0° —((1/p)(AR/R))o= +20.0X10-8, | 


With all these corrections for strain the final 
values of the tension coefficients of resistance are 


pu=— 1.1X10™™, pw=—3.3X10-%, 
p33 = — 13.9X10-”, pis= +2.5X10-", (27) 
Pu=+ 5.8X10-, pua=—0.6X10-", 


These may be considered the tension coefficients 
of specific resistance. 


26) 


CONCLUSION 


In this work with antimony, a crystal of 
trigonal symmetry, Professor Bridgman’s theory 
of the change of resistance has again been found 
to be in agreement with experimental results, as 
in the case of bismuth. The values of the six 
coefficients defining the change of resistance 
under the action of stresses have been found 
directly and have also been corrected for the 
change arising from strain. The latter procedure 
is the more defensible from a logical point of 
view, since the corrections for strain do not 
depend on the crystal orientations in exactly the 
same way as do the tension coefficients, so that 
the uncorrected readings can be expected to 
satisfy the Bridgman theory only approximately. 
However, in practice, since the corrections 
arising from change in cross section and from 
change in angle differ only slightly in their 
functional dependence on @ and ¢ from that of the 
tension coefficients and the correction for length 
not at all, both the uncorrected and the corrected 
observations can be represented within experi- 
mental error by Professor Bridgman’s equation, 
the constants involved being of course different 
in the two cases. This view is upheld by the fact 
that the probable errors are very nearly equal in 
the two cases. 

It is again a pleasure to thank the Director of 
the Laboratory, the authorities of Harvard Uni- 
versity and particularly Professor Bridgman for 
the privilege of carrying on this work in the 
Research Laboratory of Physics. 
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New Phenomena in the Change of Resistance of Bi Single Crystals in Magnetic Fields 


O. StierstapDt, Norman Bridge Laboratory of Physics, California Institute of Technology 
(Received August, 1932) 


T is well known that of all metals Bi shows the 
greatest change of resistance of electrical con- 
ductivity in a magnetic field. These changes in 
conductivity are strongly dependent upon the 
orientation of the crystal and upon the direction 
of the crystal axes with respect to the directions 
of the field and of the current. In the present in- 
vestigation it has been possible for the first time 
by means of a specially constructed crystal- 
goniometer to investigate the conductivity of a 
metal-crystal in every desired configuration of 
current, field and crystal axes, which has not 
previously been done. 

The Bi single crystal rods of 5 cm length and 
2-3 mm diameter were rotated about an axis, 
which was always perpendicular to the long-axis 
of the rod and to the magnetic field. The crystal- 
lographic orientation of the single crystal rods 
gave the following three cases: (1) The principal 
axis (111) perpendicular to the axis of the rod 
and parallel to the rotation axis (P;-crystals) ; 
(2) The principal axis (111) perpendicular to the 
axis of the rod and perpendicular to the rotation 
axis (Ps-crystals); (3) The principal axis (111) 
parallel to the axis of the rod and perpendicular 
to the rotation axis (P3-crystals). 
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Fic. 1. 7=3000 A W/cm. 


* Besides in each of these three cases one had to 
discriminate, whether the traces of the imperfect 
cleavage planes (111) in the (111)-plane were 
perpendicular or parallel to the directions of the 
current or of the field and whether such extreme 
positions could be obtained by the rotation of the 
crystal in the field. This led to a subdivision of 
the three groups in the following three sub- 
groups: P,;* and P,'!, P.* and P,!', P;* and P;''. 

Fig. 1 shows the behavior of the conductivity 
of a P;* crystal at a rotation in the constant 
magnetic field of H = 3000 A W/cm; Fig. 2 shows 
the same for the P,''-type. (g=angle between 
current and magnetic field.) 

It is evident that there is a fundamental dif- 
ference between the resistance changes of a single 
crystal and of a polycrystal in a magnetic field. 
In polycrystalline materials it is well established 
except in the case of ferromagnetic metals that 
the change of resistance is greatest, when the 
direction of the current is at right angles to the 
magnetic field (¢=90°) and is least, when it is 
parallel or antiparallel to the field (¢=0°, 180°). 
These curves, however, show that this is not al- 
ways the case for single crystals. 

Measurements on polycrystals, therefore, are 


, 
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not in general capable of elucidating the real 
laws of magnetic resistance changes of metals. 

The shapes of the curves Figs. 1 and 2 which 
repeat themselves with a period of 180°! sug- 
gested the calculation of the Fourier coefficients 
of these periodic functions. The analysis gave the 
following results: 


P,* (Fig. 1) Ao=7.7, A:=—1.35, 
Az=1.5, A3=0.12. 

P,'! (Fig. 2) Ao=9.2,  Ai=—1.7, 
Az=—1.16, A3=—0.15. 


The coefficient Ao is of little interest as it is 
independent of the rotation of the crystal in the 
field. 

The frequency A, is determined by the varia- 
tion of the angle between the current and the 
field and is in no way connected with the sym- 
metry of the crystal. It characterizes only the 
polycrystalline behavior. 

The crystallographic significance of the third 
Fourier coefficient, As, is very simple. It is due 
to the sixfold symmetry of the Bi lattice (six im- 
perfect cleavage planes of the type 111).? 

Both frequencies A; and A; with reference to 
their respective symmetry elements obey a 
cos? law. Further investigation showed that all 
the observed phenomena of the magnetic re- 
sistance variation of Bi can be expressed as a 
sum of a series of cos® terms. 

The last period, Ag, is very difficult to attribute 
to the primary structure of the Bi lattice, since it 
goes through not one or three but two cycles upon 
rotating the crystal through 180°. We shall see 
from the measurements upon crystals into which 
small amounts of impurities have been introduced 
that the Fourier analysis of the curves shows that 
the secondary structure of the Bi lattice is re- 
sponsible for the frequency Ae. 

The P2 series shows a much simpler behavior of 
the conductivity in a magnetic field than the P; 
series. Fig. 3a shows the P,''-type. None of the 
(111)-planes are brought by rotation into an 
extreme position with respect to the field, and 
therefore make no contribution to the change of 


1Complete polar-diagrams from 0° to 360° have been 
taken. 
2 Cf. Schubuikow and de Haas, Comm. Leid. 207. 
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resistance. In the P2*-case (Fig. 3b), however, 
one of the (111)-planes reaches a perpendicular 
position with respect to the field at g=18° 25’ 
(or at g= —18° 25’, depending on the starting 
position of the crystal in the field), and after a 
further rotation of 90° a parallel position. There- 
fore the maxima and minima of the curves 3b are 
shifted from their 0° and 90° positions. 

The P3; crystal shows a very similar behavior. 

The P» crystal is the first case observed for a 
nonferromagnetic material in which the magnetic 
change of resistance is greater when the current 
is parallel to the field than when it is perpendicu- 
lar. Such a phenomenon previously has neither 
been observed nor suspected but follows simply 
from the fact that the influence of one of the 
crystal planes is far stronger than that of the 
direction of the current. 

Finally, single crystals of the P\-type with im- 
purities have been investigated. These crystals 
show qualitatively the same picture as the pure 
Bi crystals, especially the appearance of the 
frequency A». Besides the general decrease of the 
resistance change with all kinds of impurities the 
frequency Ag is affected in a characteristic way by 
the foreign metal atoms. This is another reason 
for the supposition that the frequency A: is due 
to a secondary structure. The frequency Ag is 
diminished or increased with respect to the other 
frequencies A, according to whether the impuri- 
ties are to the left or to the right of the Bi 
column of the periodic system. 

It is shown in all the above measurements that 
most of the important crystal planes influence the 
resistance change in some characteristic fashion. 
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This is true of all the planes which are necessary 
for the description of the elementary cell of Bi 
lattice. Direct electric or magnetic measurements 
do not discriminate these planes. The magnetic 
change of resistance may therefore be regarded 


as a means of crystal analysis. This type of crys- 
tal structure analysis is not as generally useful as 
the optical methods but is a very sensitive one for 
detecting structure elements which have not yet 
been found optically. 
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Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
them to this depariment. Closing daies for this 
department are, for the first issue of the month, the 


Identification of the (2,2) Band of Neutral OH. Satellite 


Analysis of some plates which we have obtained with 
very intense exposures of the ultraviolet “water vapor” 
bands, due to neutral OH, permit us to make approxi- 
mately 150 new assignments of lines in the band sequence in 
the neighborhood of 43064. Among these are about 75 lines 
which form the six principal branches of the (2, 2) band 
with its principal head at about \3185. This brings the 
number of bands which have been identified in the OH 
system to a total of eight. The assignments are confirmed 
by application of the combination principle, both with re- 
spect to the upper level, by comparison with the (AF’)’s of 
the (2’, 0’) and (2’, 1”) bands, and with respect to the 
lower level by comparison with the (AF”’)’s of the (1’, 2’’) 
band, recently identified by Johnston, Dawson and Walker. 
Rotational and vibrational constants derived for the new 
(2, 2) band are in good agreement with the; values derived 
from the (1’, 2”) band. 

The remaining lines (approximately 75 in number ) in- 
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twentieth of the preceding month; for the second issue, 
the fifth of the month. The Board of Editors does 
not hold itself responsible for the opinions ex- 
pressed by the correspondents. 


Series of the (1,1) Band. Extensions of the (0,0) Band 


clude the ““R», and the Py, satellite series of 3122, 
which have, heretofore, been unobserved, as well as exten- 
sions of the principal R series of both \3122 and 3064 and 
of the satellite series of \3064. We also confirm Shaw's 
recent measurements? of 63 new assignments in the prin- 
cipal Q and P branches of \\3064 and 3122. 

A more complete description of the work will be sub- 
mitted to this journal. 


HErRrIcK L. JOHNSTON 
Davip H. DAwson 


Department of Chemistry, 
The Ohio State University, 
March 7, 1933. 


1 Johnston, Dawson and Walter, Phys. Rev. 43, 473 
(1933). 
2 R. M. Shaw, Astrophys. J. 76, 202 (1932). 


The Molecular Structure of Ozone 


The author had independently reached Badger and 
Bonner’s! conclusion that an obtuse-angled model for the 
ozone molecule is to be preferred to an acute-angled one in 
interpreting Gerhard’s* careful results on the infrared 
absorption bands of that substance. It still seems of interest 
to submit the present interpretation, which differs from 
Badger and Bonner’s in the frequency assumed for the 
missing fundamental v2, if only to emphasize the need for 
further measurements. Table I compares the various inter- 
pretations. . 


v2 =440 is thus seen to give somewhat worse agreement 
with the observed combination bands than does 760. 
Moreover, Wulf’s® assignment of the 440 cm~ interval 
between two progressions in the visible band system to v2 
of the normal state does not seem secure, as it is unsup- 
ported by measurements of temperature coefficient of 


1 Badger and Bonner, Phys. Rev. 43, 305 (1933). 
2 Gerhard, Phys. Rev. 42, 622 (1932). 
3 Wulf, Proc. Nat. Acad. Sci. 16, 407 (1930). 


TABLE I. Vartous interpretations. 











Observed Calculated 
Frequency (cm~) and type Gerhard Badger & Bonner Benedict 
— ; ve=528 (Z) ve=440 (D) v2=760 (D) 
879 (N2O; as impurity?) -— . 2v2=879 (D) —— 
. = 1033 = = 
1030-1060 (irregular D) {85 1055 ” v1 = 1050 (D) vi = 1046 (D) 
1355 (Z) Pins tb v3=1355 (Z) v3= 1357 (Z) 
ye Y3= — 2v2= 1520 (D) 
1515 (?) { 3v2= 1584 (Z) : 
2v3= 2066 (Z) 
2108 (wide, irregular) 2v2+v3= 2088 (D) 2v,;=2100 (D) { 2v, =2085 (D) 
4ve=2110 (Z) vot+v3=2117 (Z) 
2710 (?) 2v,;=2710 (Z) 2v3=2710 (D) 2v3=2710 (D) 
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absorption, and so may be a frequency difference of the 
upper state; and Wulf and Melvin‘ in other careful quan- 
titative measurements on the ultraviolet absorption, found 
no indication of such a low frequency. 

Further confirmation of the obtuse-angled structure, and 
a rough determination of the molecular dimensions, may 
be obtained from a consideration of the envelope of the 
irregular double absorption region 1030-1060 cm™. If 
yzone, with three atoms of like mass, has an apex half- 
angle >45°, the ratio of the least to the middle moment 
of inertia, A4/B=p<0.33, and the molecule is of the class 
of “slightly asymmetric” molecules whose infrared bands 
are sketched by Nielsen. Assignment of the absorption 
region to v; means that the band is of Nielsen’s type B, the 
electric moment vibrating along the middle axis of inertia. 
The general resemblance, particularly in the occurrence of 
a gap at the center of the band, is apparent; it must be 
remembered that many more lines will appear in the actual 
spectrum than on Nielsen's charts, transitions from levels 
as high as J=30 occurring with intensities of the same 
order of magnitude as lines near the center. The out- 
standing difference between the observed envelope and 
Nielsen’s sketch—the unsymmetrical appearance that 
led Gerhard to ascribe the absorption to two distinct 
vibrations—may be accounted for by interaction of vibra- 
tion with the rotation. The magnitude of such distortion 
is not so remarkable as it may seem on first sight: for in 
vy, the apex oxygen atom moves along the B axis with an 
amplitude twice the component of the amplitude of the 
end atoms along that axis. Thus on the crude mechanical 
picture of the vibration, a stretching of the bond between 
adjacent O’s sufficient to increase the effective moment of 
inertia around the B axis by AB/B=x, will increase the 
moment around the A axis by AA/A =3x. And since A is 
the most powerful constant in fixing the position of the 
modified Q-branches that contribute most to the absorption 
envelope, a small x, of the order of magnitude of the dia- 
tomic-molecule convergence factor a, will cause a very 
considerable distortion of a type B band. This mechanical 
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picture of the convergence is of course not exact, and 
further distortions are to be expected due to rotational 
centrifugal effects at the higher J values; however, en- 
velopes were constructed, using the above model and 
energy levels of the asymmetric rotator derived by an 
approximative method similar to that of Ray, and a fairly 
good duplication of the observed envelope was obtained 
with A=6.4X10-™°, B=71X10- g cm’, and x=0.015. 
These correspond to an apex half-angle of 61° and an 0-0 
distance of 1.29A. These dimensions are very rough, for the 
reasons given above and because the envelope of an incom- 
pletely resolved band may present misleading features, as 
Professor Dennison has kindly pointed out, but they are 
quite reasonable and possibly accurate to +10 percent. 

The other bands whose envelopes are known are in 
agreement with this model. »; at 1357 cm shows the 
characteristic zero branch, shifted slightly to lower fre- 
quencies by convergence, and widened by asymmetry so 
that the maxima of the P and R branches are no longer 
resolvable. In the region of 200 cm two bands, both 
presumably greatly distorted by vibrational convergence, 
overlap to give the observed irregular envelope. 

Substitution of the three fundamental frequencies 
ve = 1046, ve = 1357, 5=760 into Mecke’s’ formulas (modi- 
fied valence force system), gives real and reasonable values 
for the force constants only for apex half-angles ranging 
from 56 to 60 degrees, in satisfactory agreement with the 
angle obtained above. 

WituiaM S. BENEDICT 
Contribution No. 303 from the Research 
Laboratory of Physical Chemistry, 
Massachusetts Institute of Technology, 
Cambridge, Mass., 
March 1, 1933. 


* Wulf and Melvin, Phys. Rev. 38, 330 (1931). 

5 Nielsen, Phys. Rev. 38, 1432 (1931). 

6 Ray, Zeits. f. Physik 78, 74 (1932). 

7 Mecke, Zeits. f. physik. Chemie B16, 421 (1932). 


Vibrational Isotope Effects in Polyatomic Molecules. II 


From considerations of symmetry, Langseth’ has 
obtained expressions for the effects of isotopes of X atoms 
on the normal vibrations of pentatomic molecules YX,. 
Following the method already used for symmetrical tri- 
atomic and tetratomic molecules,? we have developed the 
Dennison equations for the normal vibrations of molecules 
YX,, obtaining relations for isotope effects due to isotopes 
of Y atoms and of X atoms. In the case of X atom isotopy, 
our results agree with those of Langseth for the inactive 
frequencies, but differ, for reasons to be pointed out in our 
complete paper, for the isotope shifts of the active fre- 
quencies. We shall merely state some of our results here, 
, in particular those independent of the force constants. 

Dennison*® has shown that the frequencies of the sym- 
metrical molecule YX, are: frequency w; (single and 
inactive), w2(double and inactive), ws; and «, (both triple 
and inactive). We shall refer to the frequency «; of the 
molecule with isotopes of either Y or X atoms as «;’. 


Let M be the mass of Y atom, m the mass of X atom, 
M-+AM the mass of isotope of Y, m+ Am the mass of isotope 
of X, and ux= M(4m+M)-. 

Consider, first, molecules 


yy, ™ y(@+4™) x (im) 


and 


Then there will be no isotope shifts of the inactive fre- 
quencies, since the motion does not involve the Y atoms. 
The isotope shift of each active frequency depends upon 
the force constants, but the following relation is inde- 
pendent of the constants: 


ws’! /wsa4 =u (1+4m/M+AM)). (1) 
1 A. Langseth, Zeits. f. Physik 72, 350 (1931). 
2 E. O, Salant and J. E. Rosenthal, Phys. Rev. 42, 812 
(1932), 
3D. M. Dennison, Rev. Mod. Phys. 3, 303 (1931). 
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For molecules YX, and YOOX, x (mtdm) 
we have different relationships, naturally. For the inactive 
frequencies: 


(2) 


The triple frequencies w; and w,of YX," are now replaced 
by two double frequencies, say w;’ and w,’ and two single 
frequencies, say ws’ and ws’, each depending on the force 
constants (this is the result essentially different from 
Langseth’s), but related by 


wy! /w, = we! /we = (1+Am/4m)~}. 


(3) 


w3'wa! /w3w4= 1 —(Am/8m)(u+ 3) 
and 


(4) 


The degeneracy of the active vibrations is, as Langseth 
pointed out, completely destroyed in the molecule 
YO X,0/ X,(™+4™) | but again we find that the frequency 


ws'we’ /w3w4= 1 —(Am/8m)(u4+2). 
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shift of each active vibration is not calculable from the 
masses alone, only shifts of certain pairs of frequencies 
being thus simply determined. This different result is, of 
course important in the use of the isotope effects for the 
assignment of the spectroscopic frequencies to the proper 
normal vibrations. 

We are greatly indebted to Dr. Louis S. Kassel of the 
United States Bureau of Mines for pointing out to us that 
the relative concentrations of the different isotopic species, 
as stated on page 433,? are in error and should read 


1:2a/b: (a/b)? and 1: 3a/b: 3(a/b)? : (a/b)*. 


The relative abundances of the C Cl, molecules with 
different proportions of isotopes of Cl were noted by 
Langseth. In general, if the relative concentration of 
X™ to X(™+4™ is b/a, then the molecules are present in 
the following proportions: 


VX,™ ; VX,(~MX(mtdm) » px Cm xX, (mt dm) » Vx(m) X (mt Am) » VX mt 4m) =1 : 4a/b : 6(a/b)? : 4(a/b)? : (a/b)*, 


Physics Department, 
New York University, 
March 8, 1933. 


E. O. SALANT 
Jenny E. ROSENTHAL 


The Nuclear Spin of Sodium 


We have made a direct measurement of the nuclear spin 
of sodium by means of the magnetic deflection of a beam 
of neutral sodium atoms.' Sufficient resolution is attained 
by the use of a magnetic velocity selector consisting 
essentially of an inhomogeneous magnetic field which 
spread the beam into a velocity spectrum, the field being 
of sufficient strength to cause a complete decoupling of the 
nuclear and electronic spins. By means of a movable 
selector slit we select a portion of this beam homogeneous 
in velocity to about 10 percent and corresponding to about 
half the velocity of the Maxwell distribution.” 

This beam’ is then permitted to pass through a weak 
inhomogeneous field of the Stern-Gerlach type. This field 
is so chosen that there is still strong coupling between the 
nuclear and electron spins. Each magnetic level is thus 
under the influence of a different deflecting force and the 
beam splits up into a number of components. The beam is 
then further allowed to traverse a strong inhomogeneous 
field of the same type, supplied by a third magnet, so 
arranged as to bring the whole deflection pattern into the 
neighborhood of the original position of the undeflected 
beam. This serves as a focussing device to sharpen up the 
lines. The length of path of the beam is 61 cm. The method 
of detection is the Langmuir-Taylor surface ionization 
detector used ballistically.? 

The beam splits up into four distinct components as 
shown in Fig. 1 which represents the experimental points 
of one particular run directly. We have checked this result 
under varied conditions. The peak marked A does not 
belong to the deflection pattern for the selected atoms. It is 
due to fast atoms from the small tails of the original beam. 
This we have proved abundantly. 

The total number of magnetic levels of the atom is 
(27+1)(2J+1) where J is the nuclear spin, in units of 
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h/2x and J is the angular momentum of the external ’ 
electron configuration. For sodium J=}, and the number 


1 Breit and Rabi, Phys. Rev. 38, 2082 (1931). 
? Rabi and Cohen, Phys. Rev. 43, 377 (1933). 
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of levels is 47+2. Our method of velocity selection in a 
strong field involves the use of only one-half the total 
number of levels. We therefore obtain 27+1 components 
in the deflection pattern. Since we have four components 
we accordingly find the spin of the sodium nucleus to be 
3h/4r. 
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A complete description of the method and results will be 
published in this journal. 
I. I. Ras 
V. W. Coen 
Columbia University, 
March 11, 1933. 


Angular Distribution of Low Energy Cosmic Radiation and Interpretation of Angular Distribution Curves 


One of us has made an experimental determination of the 
angular distribution of the cosmic radiation at two dif- 
ferent elevations,' but an interpretation of the data on the 
assumption of straight line paths through the atmosphere 
was not advocated since it was uncertain to what extent the 
distribution might have been affected by divergences 
between the paths of primary and secondary rays and by 
scattering in their passage through the atmosphere. Since 
it is known that the soft rays suffer the greater deviations* 
we have measured the angular distribution of the radiation 
which is stopped by 3.8 cm of lead in order to determine 
the extent to which such scattering affects the angular 
distribution. At every 15 degree interval between the 
vertical and horizontal measurements of the coincidence 
counting rates of three G.M. counters* were made alter- 
nately with and without a block of lead inserted just above 
the lower counter. 

Each counting rate was the average of fifteen or more 
independent data, an analysis of which showed that the 
fluctuations agreed with those to be expected statistically 
and for which the standard deviation is the quotient of the 
square root of the total number of counts by the total time. 
Two such complete runs taken several weeks apart were 
in good agreement. The differences in counting rates with 
their probable errors as calculated by combining both sets 
of data are shown by the experimental points in Fig. 1. 
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For comparison the continuous line represents the distri- 
bution of the rays which penetrated the lead plotted to a 
reduced scale. Although there is an indication that a small 
fraction of the soft rays, equal to about fifty percent of the 
horizontal rays, has been scattered through large angles, 
the majority are fully as concentrated about the vertical 
as are the penetrating rays, and we may conclude that the 


distribution of the total radiation is not appreciably 
broadened by scattering. We may, therefore, interpret the 
angular distribution of the total radiation on the basis of 
straight line paths through the atmosphere. 
Assuming a simple exponential law the theoretical dis- 
tribution is: 
j(0) =joe~** sec 6 


In Fig. 2 are plotted three sets of experimental data show- 
ing the distributions on Mt. Washington, New Hampshire, 
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’ Fic. 2. Distribution of total radiation. Points are experi- 
mental; solid curves are computed for 1=0.18 per meter 
of water for elevations sea level and 6280 ft. 


at a base station near Mt. Washington, and on the roof 
of the Bartol Laboratory. The theoretical distributions 
for sea level and the top of Mt. Washington for » =0.18 per 
meter of water are represented by the continuous curves. 
The agreement between the observed and calculated dis- 


1T. H. Johnson, Phys. Rev. 43, 307 (1933). 

?C. D. Anderson, Phys. Rev. 43, 381 (1933). 

3 T. H. Johnson and J. C. Street, J. Frank. Inst. 215, 239 
(1933). 
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tributions at sea level is remarkably good and the absorp- 
tion coefficient is in good agreement with the value 0.16 
deduced by Millikan and Cameron‘ from the portion of 
their ionization-depth curve adjoining sea level. For the 
Mt. Washington elevation there would seem to be a dis- 
crepancy between the observed intensities and those cal- 
culated from 4 =0.18. There is, however, some uncertainty 
in the relative intensities as measured at the two different 
elevations arising from a change in the efficiency of the 
counters with counting rate. Although the change in 
efficiency was measured! the uncertainty of the measure- 
ment is sufficient to account for this apparent discrepancy. 
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However, if the absorption coefficient increases with eleva- 
tion as rapidly as indicated by the data of Millikan and 
Cameron, unaccountable differences appear between cal- 
culated and observed intensities both as regards angular 
distribution and relative values. 
Tuomas H. JoHNson 
E. C. STEVENSON 
Bartol Research Foundation of 
The Franklin Institute, Swarthmore, 
March 13, 1933. 


4 Millikan and Cameron, Phys. Rev. 37, 235 (1931). 


A Neutron of High Velocity, and Energy Relations for Nuclear Disintegration by Non-Capture 


1. High velocity neutron 

One of the photographs which we have obtained of the 
disintegration of N“, by capture of a neutron from beryl- 
lium, into B" and Het‘ is remarkable with respect to the 
length of the track of the boron nucleus. The energy of the 
neutron which caused the disintegration is found by cal- 
culation to have been 16 X 10° electron-volts, which is about 
twice the value 7.8 X 10° usually! given for ‘fast’ neutrons. 
The velocity of this neutron is found as 5.610° cm per 
sec. From the mass* 9.0155 of the beryllium atom, it is 
found that 6.9 < 10° electron-volts is the amount of energy 
given by the change of mass in the reaction with the inci- 
dent a-particle, provided C” is formed, or considerably 
less if the product is three a-particles. The a-particles used 
in the disintegration of the beryllium were those from Th C’, 
which for the ordinary higher velocity have an energy 
of 8.8 10°® electron-volts. Thus, if no y-rays are emitted 
in the process of emission of the neutron, the total energy 
available is 15.7 < 10® electron-volts, or, making allowance 
for the kinetic energy imparted to the C” nucleus, the 
maximum energy of the neutron should be 14.7X10° 
electron-volts. The fact that the energy found experi- 
mentally is greater by 1.3 X 10° electron-volts might be due 
to the action on the beryllium of a particle of higher 
velocity than that assumed, but is much more probably 
due to the fact that present knowledge of the range- 





velocity relations of B" is not sufficiently exact to give a 
better agreement. The high energy found seems to indicate 
that the amount of energy which is transformed into y-rays 
when the neutron is emitted is zero or very small. 


2. Disintegration by non-capture 
If the neutron is not captured, then the reaction may be: 


N'"—>B!°+ Het, Am=0.0077 mass units (1) 
N4+C8+H}, Am =0.0038 (2) 
N4—+C!?+ H?, Am =0.0091 (3) 


However, (2) and (3) are excluded, at least in many of 
these events, by the fact that the track of the light atom 
has too great a line density of ions to be due to hydrogen. 

The mass increase in (1) is 0.0077 as given by Aston’s 
data. This corresponds to a gain of energy of 7.1X10° 
electron-volts. Thus the neutron must lose (7.1+4.4) x 10° 
electron-volts. The error of 4.4 X 10° electron-volts included 
here is the arithmetical sum of all of Aston’s estimated 
errors, which are not likely to be all in one direction, so 
that the probable error is less than this. 

The energy which must be supplied lies between 6 and 


11. Curie and F. Joliot, Nature 130, 57 (1932). 
? Bainbridge, Phys. Rev. 43, 367 (1933). 
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8X 10° electron-volts if Aston’s data are as accurate as is 
assumed by Chadwick.* In disintegrations of this type the 
kinetic energies of B'® and Het commonly amount to 1 to 
3X 10° electron-volts, so that in all about 9X 10° electron- 
volts should be obtained from the loss of kinetic energy of 
the neutron. There is thus at present no good explanation 
from the energy standpoint as to how such disintegrations 
by non-capture are able to occur, particularly in the work 
of Feather, who used the comparatively slow a-particles 
from polonium, and thus obtained neutrons of energy 
8X 10° electron-volts. 

Fig. 1 shows two views of a disintegration of a nitrogen 


A New Kind 


The complexity of the process by which neutrons are 
produced and their energy measured invites a careful 
analysis of the data obtained. Such an analysis has bearing 
on several other problems but the present note will deal 
only with the application to the Be® nucleus. A more 
complete account will soon be published. 

It appears that the experiments so far used to determine 
the upper limit to the range of protons produced by 
neutrons lead to an underestimate. The error may well be 
greater than twenty percent of the maximum range ob- 
served. It follows that Be® has more energy than two 
alpha-particles and a neutron.' From general considerations 
and especially on quantum theoretical grounds we are led 
to the conclusion that Be® should decompose spontan- 
eously, giving off alpha-particles and possibly a neutron. 
If a He of mass five is temporarily stable no neutron would 
appear for awhile at least. 

This matter has been investigated and the prediction 
verified. After a little trouble some radium-free Be was 
obtained. The activity was low enough to explain the 
failure of older workers to notice the effect. The con- 
tamination by the radium family as measured by radon 
content was less than enough to produce one percent of 
the activity. The thoron content showed, however, that 
the thorium family might have contributed twenty-five 
percent of the ionization measured. It can then be argued 
that the actinium family was also harmless. However, the 
question was decided definitely by range determinations 
with Al foils. The range proved to be about one centimeter 
in air and therefore quite easily distinguishable from the 
ordinary alpha-disintegrations. 

The corresponding half-life turns out to be about 10" 
years. This is in striking accord with the helium content 
of certain Beryls which has puzzled geologists since 
Rayleigh’s* measurements. The explanation which has been 
offered for the abnormally high helium content of Beryls 
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nucleus, seemingly without capture of the neutron, since 
the plane of the event is such that the neutron could not 
have had a straight line path from the source if the event 
had occurred with capture. 
WituiAM D. HARKINS 
Davip M. Gans 
Henry W. NEwsON 
G. H. Jones Laboratory, 
University of Chicago, 
March 13, 1933. 


3 Chadwick, Proc. Roy. Soc. A136, 692 (1932). 


of Radioactivity 


in terms of the radioactivity of Be* is now unnecessary and 
in fact certain strong objections can be offered against it. 
That is to say, we can estimate the mass of Be* from that 
of Be® or from the energy of the short range group of 
particles emitted in the bombardment of Li by protons. 
From these data it seems that the half-life of Be* would be 
extremely short and that this isotope could play absolutely 
no part in the chemistry of beryllium, it would have disap- 
peared before any Beryl crystal could have been formed. 
The helium contained in Beryl may prove extremely inter- 
esting in connection with the existence of He*. This isotope 
would exist in vanishing proportions in He from ordinary 
sources but might be a large fraction of the helium which 
results from the decomposition of Be®. 

Further investigations with counters and especially with 
the expansion chamber will, it is hoped, lead to results 
which will help to determine accurately the masses of the 
neutron and all the other light nuclei. Such information is 
essential to the development of the theory of the nucleus. 
In particular, alpha-disintegration can now be studied in 
an almost ideally simple case and defects in the present 
theory should be brought to light. 

R. M. LANGER 
R. W. Raitt 
California Institute of Technology, 
Pasadena, California, 
March 14, 1933. 


1 In an article which has just come to hand, K. T. Bain- 
bridge (Phys. Rev. 43, 367 (1933)) announces a mass- 
spectroscopic determination of Be®. His value is rather 
higher than we predicted but it is conceivable that our 
present determination of range is off enough to account 
for the discrepancy. 

2 Cf. Bulletin of the National Research Council, Physics 
of the Earth IV, p. 405. 


The Equilibrium Theory of the Abundance of the Elements 


It is possible to consider by statistical mechanics an 
assembly containing radiation, atomic nuclei, electrons, 
and neutrons; when all possible transmutations of the 
nuclei occur without the “annihilation” of any ultimate 


particles. One can calculate the abundances of the nuclei 
of the various sorts in such an assembly, when it is in 
equilibrium, in terms of the atomic masses and packing 
fractions. The simplest procedure is possible if the nuclei 








586 LETTERS TO 
(and neutrons) are made of electrons and protons; in this 
case three independent parameters, of which one is the 
temperature, suffice to determine all the equilibrium 
properties of the assembly. If the nuclei are made, on the 
other hand, of protons and neutrons as the ultimate 
particles, very nearly the same treatment will suffice. In 
the simpler of these two cases, a parameter ¢ is associated 
with the extranuclear electrons in the assembly and a 
parameter » with the extranuclear protons. Corresponding 
to any nucleus of mass number M, and atomic number N, 
there exists a parameter yu, given by 


py = EM e-NegMs exp (—0,/kT) 


where Q, is the energy required to build the nucleus s out 
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of electrons and protons. The number of nuclei present (in 
the equilibrium state) of sort s is then 


Xe=us(9/dus)= log f[us exp («/kT)] 


and there are other familiar expressions for the other 
equilibrium properties. 

It may happen that these equilibria are established 
nearly perfectly in the stars. The writer hopes to publish 
the theory shortly. 

T. E. STeERNE* 

Harvard University, 

March 15, 1933. 


* National Research Fellow. 


Conductivity of Rare Gases Irradiated by Their Own Resonance Radiations 


At the Atlantic City meeting of the American Physical 
Society the writers' reported on an investigation of the 
conductivity of rare gases irradiated with their own 
resonance radiations and on the effect of the addition of 
small traces of impurities on this conductivity. The con- 
ductivity was determined by the volt-ampere character- 
istics of a pair of parallel plane disk collectors in the 
irradiated gas. The currents obtained do not reach a 
saturation value at any of the pressures or voltages used 
(1-12 mm; 0-150 v). This observation was also indicated 
by Penning? and has been recently confirmed by Kenty.’ 
The writers have attributed this lack of saturation to the 
collection of those ions formed in the outer volume of the 
gas as well as those formed in the region between the col- 
lectors. To test this hypothesis volt-ampere characteristics 
were taken for a pair of parallel plane disk collectors fitted 
with guard rings. The collectors which were 5 mm in 
radius and had guard rings 5 mm wide were separated 5 
mm. Collectors and leads were insulated with glass except 
on the adjacent faces of the collectors and rings. When 
they were mounted at right angles to the tube axis so that 
the space between them was shielded from the direct 
radiation, then, except at the lowest pressures, no measur- 
able current (<10~'° ampere) flowed between the guarded 
collectors. The total current (to collectors and guard-rings), 
however, was the same as that previously observed between 
unguarded collectors. When the collectors were turned so 
that the primary radiation illuminated the space between 
them, appreciable currents were obtained and the much 
larger total current was of the same order of magnitude as 
before. In no case were the collector surfaces exposed to 
the primary beam. 

These results lead to two conclusions: first, that the ions 
which give rise to the observed conductivity must be 
formed in the space between the collectors; and second, 
that any metastable atoms which may produce these ions 
must also be formed in this space. It is difficult to reconcile 
these results with the propagation of resonance radiation 
through the gas by diffusion. In previously reported experi- 
ments, where collectors were found to receive measurable 


currents when shielded from the direct primary radiation, 
the collectors were unguarded and the observed currents 
probably were due to the collection of ions formed in the 
space surrounding the collectors. These ions could reach 
the collector edges. Hence the difference between the 
currents received by such an unguarded collector when 
turned parallel and then perpendicular to the primary 
beam is not to be attributed solely to the effects of the 
direct and scattered radiations. In Ne which contains a 
trace of A, with the collectors floating with respect to the 
discharge, the currents measured to the guard rings may 
be ten to fifty times those measured to the inner disks. 
Moreover, the currents measured between the guarded 
collectors exhibit decided maxima for small collector 
voltages (<10 v). These were first observed by Penning* 
who attributed them to ions diffusing in from the surround- 
ing volume when the field between the collectors was so 
small that the ions were not trapped by the guard rings. 
This explanation the writers believe to be correct. Pen- 
ning’s use of guarded cylindrical collectors did not involve 
a change in their orientation with respect to the primary 
radiation. 

If resonance radiation is not propagated to any con- 
siderable extent by scattering or absorption and re- 
radiation in cases such as these where impurities are 
present which dissipate the energy of the scattered radi- 
ation, one must look for a new explanation for the con- 
ductivity observed in a side tube at right angles to the 
main discharge tube. It is possible that some of the primary 
radiation may be reflected into the side tube by the tube 
walls. In this case the degree of reflection by the tube walls 
must be greater than has generally been supposed so that 
less of the primary radiation is absorbed by the walls than 
has been supposed. Kenty? has already suspected that this 
may be true. If the observed conductivity is due to the 
action of the direct primary beam rather than to that of 


1 Duffendack and Smith, Phys. Rev. 43, 374 (1933). 
2 Penning, Zeits. f. Physik 78, 7 and 8, 454 (1932). 
3 Kenty, Phys. Rev. 43, 181 (1933). 
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the scattered secondary radiation it is difficult to explain 
the relatively large conductivity observed in rigorously 
purified gas under conditions such that the photoelectric 
action of the primary beam is very small. However, if the 
scattered radiation plays any important rédle one would 
expect at least measurable currents for any orientation of 
the collectors. More complete and accurate data are 
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necessary before any final conclusions can be drawn. The 
experiments are being continued. 


O. S. DuFFENDACK 


The Torque on Faraday’s Magnet 


Much of the recent discussion! of the torque on Faraday’s 
revolving magnet concerns merely the choice between 
several simple ways of making the calculation, all of which 
yield correct results. On the other hand there is also a 
natural desire to understand this torque as arising from 
the same simple laws that are definitely known to hold 
for forces upon currents alone and upon magnetized 
matter alone; and in regard to this point divergent views 
have developed owing to the fact that here the torques 
upon current and magnetized matter cannot be separated 
in experiment and the distinction between them must of 
necessity rest in part upon speculative hypothesis. 

By considering simple examples it can be seen that in the 
present state of knowledge no general answer to the 
question can be given. Two different possible distributions 
of the current among the molecules of the magnet will be 
described which would be indistinguishable to an observing 
physicist and yet which are of such a nature that in one 
case the torque certainly acts upon the magnetized matter 
while in the other it just as certainly acts upon the current. 
These two examples thus confirm the conclusion reached 
by various authors? on the basis of electron theory that in 
such cases the distribution of the true forces as between 
current and magnetic matter must depend upon the exact 
path of the conducting electrons through the material. 

For simplicity let us as usual give to the magnet the 
form of a cylindrical bar uniformly magnetized longi- 
tudinally and mounted so as to revolve freely about its 
axis, and let the current enter with symmetrical distribu- 
tion about the middle of the bar and leave at a point on the 
axis at one end. We shall also suppose the external circuit 
to be so arranged that the current exerts no torque upon 
that part of itself which lies within the magnet. The total 
torque on the bar is then known to be i@/27 where 7 is the 
total current and & the flux of induction through the 
cross section at which the current enters. 


First distribution of the current 


Let the magnetized matter consist of very thin wedges 
meeting by their edges along the axis of the bar, and let the 
current flow entirely through similar wedge-shaped spaces 
left between these wedges of magnetic matter (Fig. 1a). 

In this case the observed torque arises entirely from 
forces exerted upon the magnetic matter by the partial 
field H; due to the current. In the spaces through which 
the current flows we have, besides H;, only the relatively 
weak ‘“‘demagnetizing” field H,, due to the magnet, which 
at most points in the bar has a direction opposite to that 
of the principal flux; the effect of this field, acting upon 
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the current, is only to produce a small back-torque in 
opposition to the principal one. 

At first sight one is tempted to think that, on the con- 
trary, the field of the current cannot possibly cause the 
magnet to revolve because of the symmetry, for if the 
magnet turns through an angle d@ the work done upon one 
of the Amperian circuits in it is measured by the change 
in the flux through that circuit, and there seems to be no 
such change. The displacement d@ must, however, be taken 
indefinitely small; and just as the Amperian circuit begins 
its displacement it approaches the wedges of current on 
one side and recedes from those on the other (the current 
being of course assumed to stand still). The flux through 
the circuit will therefore vary a little, and as a matter of 
fact a calculation of this variation leads at once to the 
actual torque. Displacing a wedge through d@ is equivalent 
to transferring a slice of thickness dé from one side to the 
other, say from ab to cd. If now we integrate along a path 
abcdefa, where e and f lie at the other end of the bar, 
SH;-ds=0; but f2/=0, and, since be forms part of a 
circle enclosing the entire current, JSv'Hi- ds = (d0/2x)4xi 
=2id@. Thus /H;-ds differs along the two sides of the 
wedge by 2id@ and the work done upon the displaced slice, 
if its cross-sectional area is dA, will be 2iJd@dA =id6d®,,/2x 
(J =magnetization, ®,,=42rAJ). The torque on the mag- 
netized matter is thus (i/27) fd®,, and subtracting the 


1 Zeleny and Page, Phys. Rev. 40, 299 (1932); 24, 544 
(1926); Kimball, Phys. Rev. 28, 1302 (1926). 

2Cf., eg., Kennard and Wang, Phys. Rev. 27, 460 
(1926). 
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small back torque on the current itself mentioned above 
we have, as usual, 16/27. 


Second distribution of the current 


Going now to the other extreme, let the magnetized 
matter have the form of very flat narrow cylindrical rings 
placed coaxial with the bar, and let the current flow 
alternately radially and longitudinally through similar 
cylindrical slots between the rings, as suggested schemat- 
ically in Fig. 1b; let us suppose that the current otherwise 
ignores the molecular structure and distributes itself with 
exact symmetry about the axis. (The current must also be 
prevented from slipping sidewise by a suitable constraint.) 

In this case the torque arises from the field of the magnet 
acting upon the current. For in this case Hj; is strictly 
symmetrical and the argument just given shows correctly 
that the torque on the magnetic material vanishes; all of 
the lines of B are cut by the current and the usual calcu- 
lation gives for the torque exerted on the current alone 
JS riB dr =ib/2x where = { B2ardr. 


The real case 


Now so far as one can tell by inspection a magnet might 
actually be constructed according to either of the two 
models just described. If the conducting electrons were to 
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move with entire disregard of the atoms, passing freely 
through them, then the mean force on the electrons would 
be determined by the mean value of the true magnetic 
intensity or B, and we should have an equivalent of our 
second distribution. It seems improbable, however, that 
the electrons do so move; they should tend to follow the 
spaces between the atoms (or to do something equivalent 
to this in quantum mechanics), and the result will then 
depend upon local conditions. If the interatomic space 
through which an electron slips resembles more nearly a 
longitudinal crack, we have an approach to the case of our 
first example (mean torque on magnet); if a transverse 
crack, an approach to the second (mean torque on the 
electron). It seems thus to the writer most probable that 
the torque in actual cases arises partly in one way and 
partly in the other, the amount arising from each cause 
being at present unknown. 

If this be granted, then we cannot say that the force on 
a current is definitely either jx B or 7X H. We can only 
assert that either of these expressions, if paired off with 
the proper assumption as to forces on the magnetic 
material, will give correctly the total force-action upon a 
rigid body (or an incompressible liquid). 

E. H. KENNARD 
Cornell University, 
March 20, 1933. 


Improved Calculation of Ground State of H» 


The work of Hylleraas' on He shows that a great im- 

provement in the treatment of a two-electron problem can 
_be obtained by using wave functions in which the elec- 

tronic separation enters explicitly. We have extended this 
method to the Hz molecule, and wish to communicate our 
preliminary results. 

The dissociation heat as measured experimentally is 
4.46+0.04 volts according to Richardson and Davidson,’ 
4.44 volts as given by Mulliken. Adding the zero-point 
energy, 0.27 volt, gives —4.73+0,.04 as the energy of He 
at the equilibrium nuclear separation, 1.4 Bohr radii 
(0.74A), referred to zero for separate atoms. Our computed 
values for three nuclear distances, R, are as follows, com- 
pared with the values read from a Morse curve constructed 
from spectroscopic data and passing through the minimum 
just described: 


R 1.3 1.4 1.5 
E(calc.) —4.63 —4.69 —4.59 
E(Morse) —4.67 —4.73 — 4.68 


These results were obtained with a wave function in the 
form of a series: 


¥=ZmnikpCn njape ©O+9) [Ar™A2"urius*p? +A1"A2™ U1 "ua! p? J. 
Here Ax, wi and Ae, we are the ordinary elliptical coordinates 


of the two electrons respectively, and p is the ratio of their 
separation to R/2. The summation extends in principle 


over all not negative values of the indices such that j+k 
is even, but we find such rapid convergence that no index 
need be carried higher than 2. The screening constant 6 
and the coefficients C are parameters to be varied in the 
usual way, in order to minimize the energy, and have 
different values for each R. So far, however, we have com- 
puted only with =0.75. The terms used, and their coef- 
ficients for R=1.4, are as follows: 


mnuijgkp C mnijgkp C 
00000 1 101 1 0 +40,0841 
fo 0 0 2 0 +0.5160 200 0 0 +40.0300 
0 0 11 0 —0.2134 000 0 1 +40.1550 
f1 0 0 0 0 —0.3852 001 1 1 —0,0306 
f1 0 0 2 0 —0.0333 000 0 2 —0,0115 
102 0 0 —0,0635 


A few other terms were tried and rejected, as they led 
to negligible decrease in the energy. In fact, the converg- 
ence is so rapid upon introducing the terms in succession, 
that we can state with some confidence that the inclusion 
of five or six more would bring the energy to within 0.01 
volt of the convergence limit, and that this limit is 
—4.73+0.02 volts, in complete agreement with experiment. 


1E. A. Hylleraas, Zeits. f. Physik 65, 209 (1930). 

2.0. W. Richardson and P. M. Davidson, Proc. Roy. Soc. 
A123, 466 (1929). 

3R.S. Mulliken, Rev. Mod. Phys. 4, 78 (1932). 
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Analysis of the problem indicates that the chosen value 
of 6 is very nearly the optimum for R=1.3 and 1.4, but 
that for larger R a correspondingly larger 5 would give 
slightly better results. We are now investigating this. 

It may be of interest to give the best energy, —4.2 volts, 
obtainable from a combination of terms not including the 
electronic separation explicitly. This is the same value 
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found by Hylleraas, by a long computation which he does 
not describe. 
HuBERT M. JAMES 
ALBERT SPRAGUE COOLIDGE 
Harvard University, 
Cambridge, Massachusetts, 
March 17, 1933. 


Influence of Inner Shells on Atomic Interactions 


In applying wave mechanics to the study of the inter- 
action of atoms more complex than He it has been custom- 
ary to omit the inner shells from consideration. That the 
apparent success of computations of this sort does not 
justify this omission may be seen from a study of the 
normal state of Lig. 

The first attack on this problem was made by Delbriick.! 
In this work the inner electrons were included, but so 
many approximations were made in the numerical cal- 
culation that the results were seriously affected, as can be 
seen by comparing them with those of the writer, in 
obtaining which no ‘such approximations were made. 
Bartlett and Furry? seem to misinterpret the work of 
Delbriick as showing the inner shells to be negligible and 
omit them entirely in a work which was otherwise more 
careful than Delbriick’s, obtaining a result within 5 percent 
of the experimental value for the energy of dissociation. 
This close agreement must be viewed with some suspicion, 
however, for the corresponding computation in the case 
of hydrogen, that of Heitler and London, completed by 
Sigiura, gives a result in error by 34 percent. 

There are three fundamental approximations made in 
work such as that of Bartlett and Furry: (1) an approxi- 
mation in simplifying the operator used; (2) the omission 
of the inner shells; (3) the use of atomic wave functions to 
describe the outer electrons (the Heitler-London method). 
Such work does not seem to warrant a conclusion as to the 
importance of the inner shells in the interaction. The 
writer has, therefore, carried through without computa- 
tional approximation a similar calculation including the 
inner shells and using the complete rigorous Hamiltonian 
operator. The wave functions for the outer electrons were 
the same as those used by Bartlett and Furry. The results 
of this computation are: 


Internuclear distance (A) 2.78 2.98 3.18 
Energy of interaction (volts) —0.259 —0.294 —0.309 


A Morse curve fitted to these points has its minimum very 
close to 3.18A, and the dissociation energy is thus com- 
puted as about 30 percent of the experimental value. By 
making the approximation of using an “interaction 
operator,” as Bartlett and Furry do, the agreement with 
experiment was somewhat improved (dissociation energy 
0.328 volt, equilibrium distance 3.07A). A comparison of 
these results with those of Bartlett and Furry indicates 


that the inner shells in Liz are responsible for a repulsion 
between the atoms which is decidedly important in com- 
parison with the total energy of binding. 

An attempt to better the results by adding functions of 
ionic type failed to give a perceptible improvement. By 
changing the shielding constant for the outer electrons the 
computed interaction energy can be increased by about 
10 percent. It seems, however, that no rigorous com- 
putation designed along these lines can yield for the dis- 
sociation energy of Liz a value greater than 40 percent of 
that observed. 

In view of these results it seems to be necessary to 
assume that the approximate agreement of the results of 
computations on the alkali metal molecules with experi- 
ment is due to a cancelling of errors, the neglect of the 
repulsion of the inner shells being counteracted by an 
underestimate of the attractive forces due to the outer 
shell. The interesting results of Rosen and Ikehara* must 
then be considered as of doubtful significance until some 
reason for expecting a systematic cancellation of these two 
types of error can be discovered. 

The writer has for some time been of the opinion recently 
expressed by Furry,‘ that the best way to treat diatomic 
molecules is by using one-center functions for the electrons 
of inner shells, two-center functions for the others. A 
treatment of Liz by this method, using for the outer elec- 
trons the functions applied to the treatment of H, by Dr. 
A. S. Coolidge and the writer,’ is now well under way. A 
more complete statement and analysis of the results of the 
research on which this note is a preliminary report, as well 
as a description of the methods used, will be published later 
in connection with a report on the more promising method 
mentioned above. 


Huspert M. JAMES 


Harvard University, 
Cambridge, Massachusetts, 
March 17, 1933. 


1M. Delbriick, Ann. d. Physik 5, 36 (1930). 

2 J. H. Bartlett, Jr. and W. H. Furry, Phys. Rev. 38, 1615 
(1931). 

3 N. Rosen and S. Ikehara, Phys. Rev. 43, 5 (1933). 

‘W.H. Furry, Phys. Rev. 43, 361 (1933). 

5 Preceding letter. 
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Modified Scattering by Crystalline Hydrogen Halides 


Using the Hilger constant deviation spectrograph we 
have measured the modified lines due to scattering of 
mercury lines \\4047 and 4078 by crystals of HCl and 
HBr. 

Detailed comparison with the infrared absorption band 
of crystalline HCl, reported by Hettner' at 3.66u, and with 
the liquids of HCl and HBr? will be made in our complete 
report. We may remark here, however, that the Raman 
lines appear about 40 cm™ broad, their centers being, for 
HCl 2760 cm™, for HBr 2465 cm, showing a (0, 1) vibra- 
tional energy difference of about 10 cm™ less than those 
of the liquids, and that these displacements are much too 


great to be attributed to Lorentz-Lorenz forces. As for the 
liquids,? HBr was the stronger scatter. 


E. O. SALANT 
Dixon CALLIHAN 
Department of Physics, 
Washington Square College, 
New York University, 
March 18, 1933. 


1G. Hettner, Zeits. f. Physik 78, 141 (1932). 
2 Salant and Sandow, Phys. Rev. 37, 373 (1931). 
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